![](https://img-blog.csdnimg.cn/20201014180756919.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Earth Observation
文章平均质量分 54
there2belief
这个作者很懒,什么都没留下…
展开
-
下载带坐标的高清历史影像(Google、World Imagery Wayback)
需要下些带坐标的高清遥感历史影像,但通过GEE、PIE、地理空间数据云、USGS等途径都找不到符合需求的数据。然后,在网上找了许多教程,发现目前可以通过ersi的World Imagery Wayback和谷歌地图能看到高分辨率的历史影像;但是要下载这两个来源的历史影像还是比较麻烦的,需要其他软件、科学上网等途径作为辅助。接下来我就把这两天学到的下载World Imagery Wayback和谷歌地图带坐标的历史影像分享一下,希望对在为下载高清历史遥感影像烦恼的同学有所帮助。转载 2024-01-17 10:15:20 · 6296 阅读 · 0 评论 -
rpb/rpc文件说明与matlab读取
rpb文件是用来存储用于遥感数据几何校正的RPC(Rational Polynomial Coefficients )模型的文件。类似的还有RPC文件,rpb与rpc文件只是格式不同,但包含的信息一致。其用于从图像坐标转换到地理坐标,L1级别的遥感数据是没有经过地理编码及几何校正的,因此可以通过数据包内的RPC文件进行数据几何校正(地理编码)。原创 2024-01-12 14:13:59 · 2222 阅读 · 0 评论 -
卫星RPC有理多项式模型原理与正反算实现(python)
后续也将基于这个简单的小库,做一些应用的演示。RPC可以将其理解为一个带畸变参数的相机模型,其描述的是从三维的地理坐标到二维的卫星影像坐标之间的转换关系,一般称之为从物方到像方,我们可以理解为3D到2D,这也被称之为正投影,具体为已知lon,lat,h,求得像素坐标s,l。而反向投影的则是从影像坐标系到地理坐标系,即从像方到物方,这里要注意的是,这仍然是一个3D到2D的关系,因为2D是无法升维到3D的,具体为已经s,l,h,求得lon与lat。分别为归一化的行列号,使用归一化值是为了减少计算过程中的误差。转载 2024-01-11 17:27:34 · 1327 阅读 · 0 评论 -
ESA卫星任务查询网站
欧空局的卫星任务查询网站,卫星任务十分丰富,介绍也很详尽。原创 2023-06-27 18:15:36 · 145 阅读 · 0 评论 -
GDAL重采样与裁剪图像示例
GDAL重采样,可以通过写文件时改变图像尺寸和geo_transformes的分辨率信息实现。GDAL裁切先计算索引,通过索引获取裁切后的图像,同时重新计算geo_transformes里的左上角点坐标。原创 2023-03-09 17:52:39 · 539 阅读 · 0 评论 -
机载高分辨率SAR数据(~0.1米)
美国桑迪亚(sandia)国家实验室提供一系列机载SAR数据,包括MiniSAR、FARAR等。数据分辨率4英寸,约0.1米。原始数据下载地址,数据是复数据,以不同格式存储,数据读取可用网站上提供的matlab代码,例如.gff可以用GFF_READER。原创 2023-02-16 21:43:57 · 2128 阅读 · 14 评论 -
亚米级土耳其地震影像数据下载
亚米级土耳其地震影像数据下载转载 2023-02-09 09:19:40 · 560 阅读 · 3 评论 -
微波遥感(三、SAR图像特征)
SAR是主动式侧视雷达系统,且成像几何属于斜距投影类型。因此SAR图像与光学图像在成像机理、几何特征、辐射特征等方面都有较大的区别。在进行SAR图像处理和应用前,需要了解SAR图像的基本特征。转载 2023-02-03 10:37:50 · 7967 阅读 · 0 评论 -
世界港口数据获取
目前获取世界港口信息的网站包括searates、marinevesseltraffic。searates提供全球海港数据的地图查询与基本信息查看。marinevesseltraffic 提供港口查询,航行信息、live照片和video等信息查看。原创 2022-08-25 11:41:23 · 4515 阅读 · 0 评论 -
全球机场数据获取(基本信息与影像地图)
目前比较好的数据源是OurAirports。网站目前(2022.8.25)提供全球73,106个机场的信息,功能主要包括查询、全球机场概览、机场列表、下载、评论信息和地图浏览等。原创 2022-08-25 10:56:48 · 1782 阅读 · 0 评论 -
干旱预警系统建设的国际经验与借鉴
干旱是一种比较特殊的自然灾害,持续时间长、波及范围广、治理难度大是其主要的特点。干旱作为一种全球普遍存在的自然灾害,有120多个国家和地区每年不同程度地遭受干旱灾害的威胁。我国是饱受干旱灾害之苦的国家,根据南京大学符淙斌院士提供的数据,自上世纪70年代以来每年因干旱所造成的经济损失高达上千亿元人民币。云南省从2009年秋季开始一直延续至今的60年一遇的严重干旱,已经造成超过100亿元的经济损失。世界著名的咨询公司美国麦肯锡公司于2009年11月发布的《从粮仓到旱地?》的报告指出,除非中国加大抗击干旱的投入,转载 2022-01-19 09:34:57 · 591 阅读 · 0 评论 -
遥感应用发展观察(更新)
1 综合报告金融公司报告中金 | 卫星遥感应用:新生的蓝海,崛起中的龙头(2021)太平洋证券 | 遥感智能:来自天空的“海康威视”(2021) 完整版原创 2021-12-14 10:10:43 · 127 阅读 · 0 评论 -
遥感影像数据集更新
遥感影像数据集汇总 - 飞桨AI Studio - 人工智能学习实训社区 (baidu.com)原创 2021-12-06 20:51:58 · 291 阅读 · 0 评论 -
BigTIFF is not supported by this version of GDAL and libtiff.
gdal读取大于4GB的bigtiff影像时import gdalgdal.Open("*.tif")出现错误ERROR 4: This is a BigTIFF file. BigTIFF is not supported by this version of GDAL and libtiff.解决方案:使用高版本gdal 2.2.4 +libtiff 4.0.9 +但是某些版本还可能出现问题,使用 -c conda-forge安装gdal而不是其它源,更加稳定..原创 2021-12-03 12:24:05 · 928 阅读 · 3 评论 -
自然资源部卫星遥感应用报告(更新至2020)
自然资源部每年7月份发布上一年度的自然资源部卫星遥感应用报告,是对自然资源领域典型应用的总结,对于遥感应用探索有一定参考。自然资源部卫星遥感应用报告-2020版自然资源部卫星遥感应用报告-2019版 摘录2020年版应用简介: 全国应用服务2020 年,面向自然资源调查监测、督察、执法、海洋防灾减灾、海洋生态 保护等需求,持续开展了全国自然资源常态化季度遥感监测、水资源遥感监测、 红树林遥感监测、冰川遥感监测、青藏高原湖泊变化遥...原创 2021-11-19 12:55:57 · 17699 阅读 · 0 评论 -
标准规范查询与下载
1 标准规范查询国标、行标、地方规范等可以在全国标准信息公共服务平台查询,例如国标在:国家标准目录查询- 全国标准信息公共服务平台 (samr.gov.cn)2 标准规范下载如果标准规范右侧显示全文链接,则可下载。点击全文点击在线预览打开链接后,使用另存为无法保存pdf,可以使用打印,选择输出pdf可以保存出pdf文件。...原创 2021-09-24 09:50:55 · 5074 阅读 · 1 评论 -
Landsat Fractional Snow Covered Area Product(Landsat雪比例产品,含云掩膜)
Landsat Fractional Snow Covered Area Product是一套30米的雪比例产品(含云掩膜) 、该产品包括2部分。一是雪比例产品fSCA,产品文件说明LXSS_US_HHHVVV_YYYYMMDD_yyyymmdd_CCC_VVV_PACKAGE.tar(e.g., LE07_CU_003011_200008285_20170918_C01_V01_SNOW.tar)L LandsatX Sensor (“C” = OLI/TIRS, “E” = E原创 2021-09-01 15:04:19 · 246 阅读 · 1 评论 -
土地覆盖/土地利用简介及数据集
1 简介土地覆盖:地球表面当前所具有的自然和人为影响所形成的覆盖物,是地球表面的自然状态,如森林、草场、农田、土壤、冰川、湖泊、沼泽湿地及道路等。土地利用:是人类在生产活动中为达到一定的经济效益、社会效益和生态效益,对土地资源的开发、经营、使用方式的总称。两者的区别:• 土地利用表示与土地相结合的人类活动而产生的不同利用方式,反映土地的社会和经济属性。• 土地覆盖表示地球表面存在的不同类型的覆盖特征,强调的是土地的表面形状,反映土地的自然属性。土地利用/土地覆盖分类系统LULC分转载 2021-08-23 10:00:50 · 10823 阅读 · 1 评论 -
GDAL driver查找
矢量driver:https://gdal.org/drivers/vector/index.html栅格driver:https://gdal.org/drivers/raster/index.html以矢量driver为例,代码中使用的driver是Short name,例如ESRI Shapefile、GeoJSON原创 2021-08-19 10:53:20 · 629 阅读 · 0 评论 -
核心期刊查询
中文核心期刊查询可以在知网的期刊库中查询,地址如下:出版来源导航 (cnki.net)查询时可以按类别导航选择,或直接搜索关键词,通过“核心期刊”筛选出核心期刊原创 2021-08-03 09:32:36 · 1712 阅读 · 0 评论 -
泥石流与滑坡的遥感解译特征
泥石流的解译特征 标准的泥石流由三部分组成:形成区、流通区和堆积区,如图 3-1。其中,堆积区是识别泥石流的重要标志,同时也是评估灾害危害程度的重要指标。 泥石流的形成区作为泥石流固体物源的集聚区,植被覆盖度较差,土壤裸露程度较高,表面岩体风化侵蚀严重,多为三面高地环绕的盆地低谷地带,且多群发小规模崩塌或滑坡等不良地质现象。在遥感影像上,形状多呈漏斗状、瓢形或桃叶形;色调相比周围较浅,且面积较大;由于表面岩体破碎严重,纹理粗糙,且伴有繁杂的雨水冲刷细沟。...转载 2021-01-05 19:50:15 · 6490 阅读 · 3 评论 -
gdal 栅格转矢量线
This recipe converts raster pixels with a specified value to vector lines. For example the blue pixels (value = 0) are converted to vector lines.import ogr, gdal, osr, osimport numpy as npimport itertoolsfrom math import sqrt,ceildef pixelOffset2co转载 2020-12-30 14:02:50 · 1531 阅读 · 2 评论 -
pyshp写投影信息
pyshp库对shp文件的操作本身不包含投影信息,如果需要增加投影信息的话需要单独写.prj文件1、使用pyshp写shp文件(包含shp, dbf, shx)"""using pyshp create shpfile(1) type:You can reference shape types by the numbers or by constants defined by PyShp:shapefile.NULL = 0 shapefile.POINT = 1 shapefile.POL原创 2020-12-27 19:19:30 · 689 阅读 · 7 评论 -
欧洲主要城市变化数据集(Urban Atlas LCLU Chang)
欧洲环境局提供了时序10米分辨率的主要城市地表覆盖变化数据集(Urban Atlas LCLU Chang)。The European Urban Atlas provides reliable, inter-comparable, high-resolution land use and land cover data for 785 Functional Urban Area (FUA) for the 2018 reference year in EEA39 countries. The spat原创 2020-12-12 17:22:20 · 798 阅读 · 0 评论 -
Arcpy批量合并多个gdb中相同的要素类
Arcpy批量合并多个gdb中相同的要素类,流程为:1、将所有的gdb放在一个目录下,即代码中的inpath (path_allgdb_except_one)2、剪切一个gdb作为合并后的输出gdb,即代码中的outgdb (one.gdb)3、执行下面脚本 (arcgis10.3成功)# -*- coding: cp936 -*- #import osimport sysimport arcpyinpath = r"D:\path_allgdb_except_one"outgd原创 2020-08-13 23:09:41 · 3011 阅读 · 4 评论 -
Landsat Fractional Snow Covered Area Product(Landsat雪比例产品,含云掩膜)
Landsat Fractional Snow Covered Area Product是一套30米的雪比例产品(含云掩膜)、该产品包括2部分。一是雪比例产品fSCA,产品文件说明LXSS_US_HHHVVV_YYYYMMDD_yyyymmdd_CCC_VVV_PACKAGE.tar(e.g., LE07_CU_003011_200008285_20170918_C01_V01_SNOW.tar)L LandsatX Sensor (“C” = OLI/TIRS, “E” = ETM+,原创 2020-07-03 20:00:15 · 448 阅读 · 0 评论 -
GDAL 矢量属性数据修改(python)
Case:需要给一个现有的shp数据创建一个字段,并将属性表中原有的一个文本类型的属性转换为整型后填入新创建的字段。Problem:新字段创建成功,但是赋值操作无效,即无法成功给字段写入值。solution:对字段进行赋值后需要,重新写入Feature,否则赋值无效,即layer.SetFeature(feature)。特别注意:在对数据进行读写操作,一定要以读写的方式打开,即Open(filePath,1),该方法的原型为Open(pszName,int bUpdate = false),并且转载 2020-05-13 08:56:54 · 1322 阅读 · 0 评论 -
gdal 空间参考(投影)操作 - 栅格
获取已有栅格文件的空间信息(geo_transform, projection)import gdal, osrrasterfn = 'rasterfn.tif'newRasterfn = 'newRasterfn.tif'raster = gdal.Open(rasterfn)geotransform = raster.GetGeoTransform() # geotransformprojection= raster.GetProjectionRef() # projection原创 2020-05-10 10:56:24 · 1807 阅读 · 0 评论 -
gdal 空间参考(投影)操作 - 矢量
Create Projectionfrom osgeo import osrspatialRef = osr.SpatialReference()spatialRef.ImportFromEPSG(2927) # from EPSGReproject a Geometryfrom osgeo import ogrfrom osgeo import osrsource = osr.SpatialReference()source.ImportFromEPSG(292转载 2020-05-10 10:44:35 · 2773 阅读 · 2 评论 -
使用GeoServer+QGIS发布WMTS服务 | Publishing WMTS Service Using GeoServer+QGIS
Web GIS系列:1.搭建简易Web GIS网站:使用GeoServer+PostgreSQL+PostGIS+OpenLayers32.使用GeoServer+QGIS发布WMTS服务3.使用GeoServer+OpenLayers发布和调用WMTS、Vector Tile矢量切片服务.4.Leaflet入门:添加点线面并导入GeoJSON数据使用GeoServer发布服务请...转载 2020-05-01 14:50:01 · 1616 阅读 · 0 评论 -
Python矢量操作的几种方式-以矢量相交为例
矢量操作可以使用 ogr这类基础方法进行,但同时还有一些高级库支持更便捷的矢量操作,例如Fiona(shapely)、pyshp(shapefile)、geopandas等。教程资料导航ogr教程资料:1.GDAL/OGR例子合集2.GDAL官网API3.网上的学习笔记,教程来自于犹他州立大学4.中国OSGeo教程shapely教程资料:shapely空间分析...原创 2020-05-01 14:36:53 · 5183 阅读 · 0 评论 -
栅格切片与WMTS服务发布
WMTS – OGC Web Map Tile Servicehttps://gdal.org/drivers/raster/wmts.html切瓦片https://www.cnblogs.com/wang985850293/p/5381684.html服务发布QGIShttps://blog.csdn.net/sf2gis2/article/details/52854059...原创 2020-04-29 09:01:26 · 1523 阅读 · 0 评论 -
国家2000坐标系CGCS2000(EPSG4490和EPSG4479)参数
EPSG:4490China Geodetic Coordinate System 2000 (CGCS2000)经纬度坐标系Unit:degree (supplier to define representation)Geodetic CRS:China Geodetic Coordinate System 2000Datum:China 2000Ellipsoid:...原创 2020-04-27 14:18:02 · 37815 阅读 · 6 评论 -
spacenet 卫星数据集亚马逊云链接更新
语义分割类:SN1_buildings、SN2_buildings、SN3_roads、SN4_buildings、SN5_roads、SN6_buildings通过以下方式查看aws s3 lss3://spacenet-dataset/spacenet/SN1_buildings/aws s3 lss3://spacenet-dataset/spacenet/SN2_build...原创 2020-03-26 17:27:02 · 2271 阅读 · 2 评论 -
python gdal,ogr实现geojson转shpfile
本文修改自:【python】geoJson转换Shapefilegdal,ogr实现geojson转shpfileimport jsonimport gdalfrom osgeo import ogrdef create_polygon(coords): ring = ogr.Geometry(ogr.wkbLinearRing) for coord in co...原创 2020-03-25 14:37:07 · 2446 阅读 · 0 评论 -
python gdal ReadAsArray 按块读取栅格
gdal ReadAsArray支持按块读取影像,有两种方式:(1)针对datasetdataset = gdal.Open(file, gdal.GA_ReadOnly)data = dataset.ReadAsArray(xoff=0, yoff=0, xsize=512, ysize=512)(2)针对banddataset = gdal.Open(file, gdal...原创 2020-03-20 11:54:38 · 14753 阅读 · 5 评论 -
减灾云服务平台 IN-CORE
IN-CORE:Run your scientific analyses that model the impact of natural hazards on a community and the resilience of those communities.TheNational Institute of Standards and Technology (NIST)funded ...原创 2020-03-05 18:31:13 · 356 阅读 · 0 评论 -
Building damage dataset
1.xView2-xBDhttps://xview2.org/register to downloadxBD: Annotated high-resolution satellite imagery for building damage assessment6 disaster types,850,736 buildings, 45,361 sq km2.xVi...原创 2020-03-04 22:52:33 · 536 阅读 · 0 评论 -
Remote sensing for natural disasters: Facts and figures(遥感在自然灾害中的应用)
By:Sian LewisSian Lewisexplains how remote sensing can be used to help manage natural disasters and highlights ongoing efforts and obstacles.Severe geo-physical or climatic events, including ear...转载 2020-03-03 21:33:49 · 1665 阅读 · 0 评论 -
arcpy.da.SearchCursor RuntimeError: cannot open '.shp'
arcpy在使用arcpy.da.SearchCursor读取shp时报错RuntimeError: cannot open '.shp',代码如下:#coding:utf8import arcpypath=r'D:\work\data\temp\test.shp'with arcpy.da.SearchCursor(path,'FID') as rows: print ro...转载 2019-12-23 13:17:35 · 3854 阅读 · 13 评论