基于Python的Grib数据可视化           利用Python语言实现Grib数据可视化主要依靠三个库——pygrib.numpy和matplotlib.pygrib是欧洲中期天气预报中心(ECMWF)的GRIG API C库的Python接口,通过这个库可以将Grib数据读取出来:numpy是Python的一种开源的数值计算扩展,这种工具可用来存储和处理大型矩阵:matplotlib是python著名的绘图库,它提供了一整套和matlab相似的命令API,十分适合交互式地进行制图:

MATLAB地图工具箱学习总结(一)从地图投影说起 前言 本学期地图投影课上,李连营老师建议我们使用MATLAB完成每周的作业.从大二上学期开始接触MATLAB学习数学运算和地理数据处理的我,自然不会放过这次机会,每周找了点时间好好研究了一下,把作业比较轻松地就搞定了.不过由于网上相关的教程.资料比较少,在夹杂了多种专业词汇的同时,又没有中文翻译,也是让人感到头疼.在此想总结一下相关的函数,以帮助其他需要使用工具箱的人.作为一个学生,我所接触的仅仅是皮毛,也必然会有一些错误,希望看到的人能够指出

地图投影(Map Projection) http://baike.baidu.com/view/94066.htm 概念: 地图投影是把地球表面的任意点,利用一定数学法则,转换到地图平面上的理论和方法. 由于地球是一个赤道略宽两极略扁的不规则的梨形球体,故其表面是一个不可展平的曲面,所以运用任何数学方法进行这种转换都会产生误差和变形,为按照不同的需求缩小误差,就产生了各种投影方法. 方法: 1.几何透视法:几何透视法是利用透视的关系,将地球体面上的点投影到投影面上的一种投影方法.如假设地球按比

一.简介 目前Web地图已经是一个非常普遍的应用,百度地图,高德地图等等极大的方便了我们的生活和学习.本项目主要是在Web完成一个简单的中国地图的绘制,实现Lambert投影. 二.制图源数据信息及来源说明 本项目所用数据来自老师所提供的1:4M 中国地图,坐标系为地理坐标系,北京54坐标系. 三.数据收集 在QGIS中用老师所提供的Python程序将源数据中的中国边界图层(线图层)和中国省会图层(点图层)导出为文本文件. 四.程序设计 1.编写程序分别读取线文件和点文件,绘制中国地图 2.编制

PROJ.4学习——地图投影(坐标系投影) 前言 PROJ是由大量的基础投影库构成.这里主要讨论学习PROJ库的相关参数. 这里大部分是讲如何将3D坐标系投影到2D平面上.投影时,涉及到基准线,单位,比例英子,偏移量,轴转换等. 很多实际操作中却涉及到3D坐标系转换为其他3D坐标系(下一篇更新). 参数列表 下面是可以应用于大多数坐标系定义的PROJ参数列表.此表不尝试描述特定于特定投影类型的参数.   参数 描述  +a  椭球体长半轴长度  +axis  轴防线  +b  椭球体短半轴长度

地球椭球体(Ellipsoid) 众所周知我们的地球表面是一个凸凹不平的表面,而对于地球测量而言,地表是一个无法用数学公式表达的曲面,这样的曲面不能作为测量和制图的基准面.假想一个扁率极小的椭圆,绕大地球体短轴旋转所形成的规则椭球体称之为地球椭球体.地球椭球体表面是一个规则的数学表面,可以用数学公式表达,所以在测量和制图中就用它替代地球的自然表面.因此就有了地球椭球体的概念. 地球椭球体有长半径和短半径之分,长半径(a)即赤道半径,短半径(b)即极半径.f=(a-b)/a为椭球体的扁率,表示椭球

此处直接给出原文链接: 1.地图投影的N种姿势 2.GIS理论(墨卡托投影.地理坐标系.地面分辨率.地图比例尺.Bing Maps Tile System)

[前言]最近在编写一个气象应用程序,用来显示某一时刻某一地区的气温等值线和风场,程序主要用到了第三方库matplotlib及Basemap.在编写的过程中发现,如果不进行擦除操作直接绘制新的等值线或风场,新的等值线(风场)会与原来的等值线(风场)叠加在一起,而绘制的等值线及风场没有单独的remove方法,所以如果想要擦除已经绘制的等值线就要将地图重新投影一遍,如果地图投影精度高一点,整个投影过程就会特别漫长.通过对等值线及风场的返回结果进行研究,我找到了一个不必重新投影地图就可将等值线及风场擦除

1. 通常GIS项目涉及到的坐标系 (1)面向局部区域的大比例尺二维平面:高斯投影(横轴墨卡托) 说明:在市一级的小范围区域的GIS系统,比如规划局.国土局.建设局的系统,大都使用高斯投影,以便与地方地形图测绘.工程报建一直采用的坐标系相一致.高斯投影的特点是投影分带,适合小范围局部,不适合应用于大省.全国等大范围应用,若是强制按某带投影,则远离中央经线的区域的角度.距离.面积全部变形严重. (2)省级.国家级范围:兰伯特投影(圆锥投影) (3)国际级.世界地图:经纬度等间隔直投 说明:“经纬度

转载地址:http://blog.sina.com.cn/s/blog_517eed9f0102w4rm.html 一篇题为<我们看到的地图一直都错得离谱……>的文章在朋友圈里莫名流行起来,里面列举了横麦氏地图的种种"谬误".小编知道,其实地图册里最常见的世界地图并非横麦氏地图,而是这样的: Fsdfafd 但有人指出:"等差分纬线多圆锥投影是一种任意投影,所以也是没法儿保障等积的……" 地球是一个3维的球体(并不准确但暂且管它叫球体吧),而地图是一个2

特别感谢:1.[张天旭]的D3API汉化说明.已被引用到官方站点: 2.[馒头华华]提供的ourd3js.com上提供的学习系列教程,让我们这些新人起码有了一个方向. 不得不说,学习国外的新技术真的是一个非常艰苦的过程. 在学习D3绘制地图的过程中.有朋友建议看一下当中投影的说明比較好.于是,凭借我这半吊子不到的英文水平,大致给翻译了下来,仅供參考: 原文链接:https://github.com/mbostock/d3/wiki/Geo-Projections#albers D3中一共提供了1

我国的基本比例尺地形图[1:5千.1:1万.1:2.5万.1:5万.1:10万.1:25万.1:50万.1:100万]中,大于等于50万的均采用高斯-克吕格投影[Gauss-Kruger]:小于50万的地形图采用正轴等角割园锥投影,又叫兰勃特投影[Lambert Conformal Conic]:海上小于50万的地形图多用正轴等角园柱投影,又叫墨卡托投影(Mercator),我国的GIS系统中应该采用与我国基本比例尺地形图系列一致的地图投影系统. 相应高斯-克吕格投影.兰勃特投影.墨卡托投影需要

在用axesm函数创建地图坐标系的时候可以指定地图投影(设置projinfo参数),地图投影可以通过projinfo函数来创建,里面的参数依据proj4投影字符串,可以参考此网页:http://remotesensing.org/geotiff/proj_list.投影后的地图坐标系显示范围的设置可以用axism函数,参数是一个列表(list),列表里的4个元素分别是起始.结束经度和纬度.脚本程序: #Set data folders basedir = 'D:/MyProgram/Distri

update20090601:EPSG对该投影的编号设定为EPSG:3857,对应的WKT也发生了变化,下文不再修改,相对来说格式都是那样,可以到http://www.epsg-registry.org 网站输入SRID进行查询.         Google Maps和Virtual Earth等的流行程度不用多讲,然而他们所使用的Web Mercator或Spherical Mercator在很长一段时间内并没有被EPSG的投影数据库所接纳.EPSG认为它不能算作科学意义上的投影,所以只是给

背景 近期开发以MongoDB为基础的分布式地理数据管理平台系统,被要求做一个简单的demo给客户进行演示.于是笔者便打算向数据库中存储一部分瓦片数据,写一个简单的存取服务器,使用Openlayers客户端加载瓦片进行展示.本来是个挺简单的功能,可是谁曾想,"打了一辈子鹰,却被鹰啄瞎了眼".在使用Openlayers客户端进行瓦片加载时,出现了加载的瓦片偏移.行列号不正确的问题. 问题描述 如图所示,地图出现了上下偏移的问题,并且,这些瓦片的实际行列号与当前Openlayers计算的行

背景 近期开发以MongoDB为基础的分布式地理数据管理平台系统,被要求做一个简单的demo给客户进行演示.于是笔者便打算向数据库中存储一部分瓦片数据,写一个简单的存取服务器,使用Openlayers客户端加载瓦片进行展示.本来是个挺简单的功能,可是谁曾想,“打了一辈子鹰,却被鹰啄瞎了眼”.在使用Openlayers客户端进行瓦片加载时,出现了加载的瓦片偏移.行列号不正确的问题. 问题描述 如图所示,地图出现了上下偏移的问题,并且,这些瓦片的实际行列号与当前Openlayers计算的行列号并不一

SpatialReference mSR4326 = SpatialReference.create(4326); SpatialReference mSR3857 = SpatialReference.create(3857); 1. 屏幕坐标转换成投影坐标 @Override public boolean onSingleTap(MotionEvent point) { SpatialReference sr = map.getSpatialReference(); Point dp = m

VPB和OSGGIS安装 转自:http://blog.sina.com.cn/s/blog_668aae780101k6pr.html 第一部分VPB安装 VirtualPlanetBuilder是一种地形数据库创造工具,是能够阅读各种地理图像和高程数据并建立小面积地形数据库,乃至大规模如整个地球的庞大数据库.这些数据库可以上传到互联网,并能够提供像在线googleEarth一样的风格漫游整个地球的数据库,或保持对本地磁盘高速接入等所需的专业飞行模拟器.前面所讲的VirtualPlanetBu

PROJ4初探(转并整理格式) Proj4是一个免费的GIS工具,软件还称不上. 它专注于地图投影的表达,以及转换.采用一种非常简单明了的投影表达--PROJ4,比其它的投影定义简单,但很明显.很容易就能看到各种地理坐标系和地 图投影的参数,同时它强大的投影转换功能,也是非常吸引人的.许多的 GIS软件中也将其集成在内.Proj可以在 window的命令下有可运行的 EXE文件,其实它更主要的是一个库!可以用来编一些批处理.在 Linux下除了可以直接运行外,还可以作为库来进行更高功能的开发.

MATLAB地图工具箱学习总结(四)自定义投影 这是本系列的最后一篇文章,准备给大家讲讲自定义投影怎么做.在做这项作业的时候,自己也是花了不少时间,将所有地图投影源文件都看了一遍,简单分析了一下源代码,就开始着手修改了.虽然也曾画出了一些奇形怪状的“艺术品”,但最终还是找到了画图的一些诀窍,使得自定义出来的投影即使会有bug,但大体上还算能看得过去. 在这里呢,我就想以最简单的一个地图投影源文件开始介绍怎么修改. 首先,让我们找到自己MATLAB安装目录,依次点击toolbox->map->m

MATLAB地图工具箱学习总结(三)地图工具箱的基本知识 今天想要介绍的是一些比较基础的函数.了解了这些函数,地图投影的基本概念才能真正明白.而要想继续研究MATLAB中有关地图投影的函数,尤其是未来我要提到的投影文件源代码,知晓这些函数的功能必不可少.本篇文章将会罗列三个案例,并在后面一一进行讲解. 1                    作业案例:地图投影作业1 这次的案例从作业1开始.作业1是要求计算出地球椭球体的一些基本参数,包括子午圈曲率半径.卯酉圈曲率半径.平均曲率半径和纬圈半径等