GIS
Theropod
这个作者很懒,什么都没留下…
展开
-
IDL调用MRT批处理MODIS
IDL调用MRT批处理MODIS产品我能说为了在IDL里批量拼接MODIS我查和试了有20个小时吗? 总之就是envi自己的函数都不怎么样。。。调用envi的seamless mosaic看起来不错,帮助里说现在都用envi::openraster和ENVIMosaicRaster。 然而前者打不开hdf4,只有在envi里手动选;后者又要求输入必须是object reference,但用老的E原创 2016-10-09 03:26:54 · 3857 阅读 · 1 评论 -
将几种北京建成区数据下载裁剪重分类后以供比较
Beijing-Impervious-Compare4 geotiff representing the impervious surface of Beijing (30m)1是建成区,其他是0urban_beijing_classifiedBeijing City Lab自己做的2015年30m城市制图 用GEEhttps://www.beijingcitylab.com/...原创 2018-12-15 01:25:59 · 1255 阅读 · 2 评论 -
【已经看不懂了】曾经记录的在c++上用gdal的一些经验
当时随便记录了一下没有组织起来,现在已经看不懂了。。。可见及时整理的重要性编译好的gdal平台版本(x86,64)需要和项目属性里面以及和调试的时候选择的一样。。。否则是链接错误 在x64/debug下面要放上所有的dll 否则报错比如找不到gdal202.dll(202是版本号)gdal allregister环境变量里面要有GDAL_DATA输出路径不能已经存在文件数组指针,函...原创 2018-03-23 11:06:36 · 1291 阅读 · 0 评论 -
【待细看】关于GDAL的VRT格式
最近才知道Gdal有这个功能,简单来说就是用一个xml描述对一个或一组栅格的处理(投影、mosiac、各种计算),生成一个.vrt文件,这样并没有真的进行计算,也没有存储文件,只有当要使用vrt的时候这些操作才会进行。比较省心。 Lazy raster processing with GDAL VRTs 官方文档 也可以导入到postgis里面...转载 2018-03-13 14:04:30 · 4559 阅读 · 0 评论 -
postgis不同的栅格存储方式和查询的速度以及区域统计
很多应用都需要综合处理栅格和矢量数据进行统计,如果两者在同一个数据库里就很方便。 但是有些影像太大了,而且还经常得更新,所以我这次想用文件式的管理,只把索引入库,一查postgis是支持这种的。raster2pgsql -s 3857 -R -C -I -M -F "raster.tif" -t 50x50 public.tablename |psql -U postgres -d da...原创 2018-03-16 17:19:14 · 3549 阅读 · 0 评论 -
关于openlaeyrs获取谷歌卫星地图的无偏瓦片
谷歌地图的卫星影像在国内是偏移的,比如国内能访问的那个http://www.google.cn/maps,还有Google Map API上专门提供给中国的http://maps.google.cn/maps/api/的这两个都是偏移过的,而https://www.google.com/maps/,以及https://maps.googleapis.com/maps/api就是没有偏移的。 一个验原创 2017-05-26 23:02:13 · 8176 阅读 · 3 评论 -
Geoserver的CORS问题解决
这已经是第四次了,每次都不太能记住Geoserver默认CORS是关的,所以读发布的地图中的数据和提交给geoserver数据的时候会有CORS错误。做两个事情:到geoserver目录下的webapps/geoserver/WEB-INF/web.xml里把cors的filter打开(有提示,uncomment to enable CORS)这时重启geoserver会报错找不到 org.ec原创 2017-09-13 23:00:01 · 951 阅读 · 0 评论 -
ArcGIS依照谷歌地球空间校正
从谷歌地球上能得到经纬度,但是空间校正工具的链接表里都是投影后的坐标格式的。所以不要试图从经纬度计算链接表里的坐标,而是把得到的经纬度做成图层,工具里手选链接快一些原创 2017-06-23 16:35:47 · 1429 阅读 · 0 评论 -
论如何让一个人心态爆炸-尝试改变奥维地图浏览器下载的地图源
小学期实习需要做一个土地利用调查,在根据影像绘制好地块的面之后,将其转成kmz导入到一个奥维互动地图的app(因为据说野外用这个手机上编辑属性方便一些?),发现图斑有偏移。。。 有可能是给我们的影像没做好几何校正,矢量化的时候描歪了,也有可能是火星坐标系的问题,因为看了一下奥维下载的地图清一色都是国内的源。 由于这个奥维的app可以使用其PC客户端离线下载的地图,我就尝试在PC端设置成khm3.原创 2017-06-21 22:39:14 · 4933 阅读 · 0 评论 -
使用Openlayers v4时感觉需要记一下的东西
ol.map有一个forEachFeatureAtPixel的函数,可以用来响应鼠标事件,得到feature的属性。但是对于TileWMS的source,或者说用这个source的Layer,这个函数没起作用。没有看源码,估计是因为WMS图层通过http请求tile来生成地图,而这个函数不包含这种方法。官方的示例用getGetFeatureInfoUrl代替,得到的是服务器返回的一个html元素。原创 2017-03-06 20:21:50 · 1711 阅读 · 0 评论 -
Shapefile 乱码与字符截断问题
Shapefile 乱码与字符截断问题把shapefile转成geojson的时候中文字段名会出现错误,是编码方式的问题。看这个博客转载 2017-02-23 12:56:19 · 795 阅读 · 0 评论 -
MySQL关于空间数据的一些入门
http://m.blog.csdn.net/article/details?id=23099407转载 2016-12-15 08:50:11 · 1816 阅读 · 0 评论 -
操作记录-在leaflet中加速大量geojson/topojson多边形的显示
场景需要展示覆盖北京的25万+个规则的网格,进行常规的点击显示数据、绘制勾选、着色等等操作。但是希望能直接使用geojson/topojson(原来构想的是随机森林和pagerank跑出结果,匹配到多边形,实时输出到geojson),不转格式,同时想尝试后端越少越好。。。实现首先想到的自然是类似于geojson-vt的实现,将json转换成切片之后显示。实际用了leaflet.vector...原创 2019-01-06 20:28:12 · 3888 阅读 · 1 评论