- 博客(30)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 读论文《Learning to Measure Changes: Fully Convolutional Siamese Metric Networks for Scene Change Detec》
CosimNet：直接量测变化的场景变化检测模型
2022-07-07 16:29:45
1140
原创 读论文《SNUNet-CD: A Densely Connected Siamese Network for Change Detection of VHR Images》
本文提出了一个用于变化监测的稠密连接孪生网络，命名为SNUNet-CD（联合孪生网络与Nested-UNet）。SNUNet-CD通过编码器与解码器之间紧密的信息传输减少了深层网络训练中位置信息的丢失，同时还提出了一种集成通道注意力机制（ECAM）进行深层监督。通过ECAM可将不同语义级别的最具代表性的特征进行细化，并用于最终的分类。实验结果表明，本文所用方法极大地改进了许多评估标准，与其他先进的变化检测算法相比，本文所用方法在精度和计算量上有更好的权衡。...
2022-06-14 16:57:54
3547
9
原创 Python计算图像纹理-灰度共生矩阵
基于Python探究灰度共生矩阵(GLCM)那点事儿 - 知乎一、什么是灰度共生矩阵？ 灰度共生矩阵（Gray-level co-occurrence matrix；GLCM）和相关的纹理特征计算是图像的一种分析技术。给定一个图像，该图像由各自具有一定强度（特定灰度级）的像素组成，GLCM矩阵在图像…https://zhuanlan.zhihu.com/p/367213524GitHub - 1044197988/Python-Image-feature-extraction: Python实现提取图.
2022-05-26 15:51:11
6714
1
原创 python ASF API下载Sentine-1数据
说明：后面再补说明import os,sysfrom selenium import webdriverimport timeimport base64import pandas as pdimport geopandas as gpdfrom urllib.request import build_opener, Requestfrom urllib.request import HTTPHandler, HTTPSHandler, HTTPCookieProcessorfrom htt
2021-10-19 11:09:45
700
原创 基于RTX30显卡的ArcGIS Pro2.8深度学习环境配置
问题描述基于不明原因，ArcGIS Pro2.8并不支持RTX30显卡，在深度学习模型训练中会出现模型精度极低，loss等值为nan的情况，具体可看ESRI讨论贴：https://github.com/Esri/deep-learning-frameworks/issues/17帖中总结出的问题具体如下：解决思路进入正题，配置argis pro2.8深度学习环境以支持RTX30显卡的思路是：1、安装CUDA11和对应cudnn；2、arcgis pro2.8深度学习环境库包替换。1、
2021-10-14 19:10:14
5759
30
原创 python LAADS+Selenium下载MODIS数据
说明：1、NASA数据的存放地址参照下图，各类数据存放在不同的网站，MODIS的L1级数据存放在LAADS网站上，此次从此网站进行数据查询和下载，需要在Earthdata网站进行注册并登陆，登陆之后在profile拿到token。图片链接：https://earthdata.nasa.gov/eosdis/daacs2、LAADS数据查询需要在网页上选择各项参数，然后显示数据列表，数据列表可下载为csv或者json文件。具体查询可以直接到网站进行操作体验：https://ladsweb.modaps
2021-06-07 11:27:31
576
2
原创 python landsatxplore下载Landsat8和Sentinel-2数据
说明：1、通过landsatxplore进行query和download，需要在earthexplore网站注册2、查询需要landsat8和sentinel-2的轨道文件，轨道文件与查询范围相交查询影像，landsat8的轨道文件网上较多，sentinel-2的轨道文件在另一篇博客中有提供下载地址（下载为kml，可转为shapefile文件）3、参数：数据集名称（按规定填）、开始时间+结束时间（年-月-日格式）、云量、查询范围（shapefile文件，只查询第一个面要素）、存储文件夹（自定义，会在文
2021-06-07 10:52:24
1604
8
原创 哨兵二号轨道数据下载
哨兵二号轨道数据下载https://sentinels.copernicus.eu/web/sentinel/user-guides/sentinel-2-msi/naming-convention
2021-05-25 16:26:42
586
原创 ArcGIS Pro发布服务
ArcGIS Pro发布服务向联合服务器发布服务，enterprise环境即server与portal进行联合托管。在arcgis pro版arcpy中，将原来arcmap版arcpy的arcpy.mapping模块迁移至arcpy.mp。1、使用arcpy.mp.arcgisproject()函数引用 ArcGIS pro工程 (.aprx) 。2、根据需要是否要进行注册数据，若要注册可以使用AddDataStoreItem()将文件夹或数据库注册到 GIS Server 站点中，在注册之前也可以
2021-03-19 16:18:48
2422
原创 python读取Himawari8数据
中间把def read_hsd(inputfile): resolution = int(inputfile[-12:-10]) if resolution == 10: rows = 1100 cols = 11000 elif resolution == 20: rows = 550 cols = 5500 else: rows = 2200 cols = 22000
2020-05-29 15:26:27
1905
6
原创 IDL MODIS 生成查找表
author:PQPRO MODIS_LUT;,month,day,iwave,idatm,iaer,lutpath,lutname igeom=0;自定义几何条件 phi0=0;卫星方位角 month=10;月份 day=27;；日期 idatm=2;：大气模式中纬度夏季 iaer=1;：气溶胶模式大陆型 v=0;；能见度 xps=-0.4;；目标物高度 xpp=-1000;；星测 iwave=44;：自定义1输入波段范围和反射相函数42为modis的red
2020-05-29 15:05:48
1712
23
原创 批量裁剪栅格 python gdal
arcpy裁剪的缺点：1.arcpy不支持中文2.arcpy裁剪会对影像重采样，裁剪前的原图像元不对应，有偏移.代码功能：以一个shp裁剪多个栅格输入：1.shp全路径和为栅格命名的依据字段2.栅格文件所在文件夹和栅格名称的后缀（.tif/.dat）3.裁剪后的输出路径（文件夹）输出：在输出文件夹下以原栅格名称+‘_’+命名字段内容+‘_clip.tif’保存裁剪后的图像import ogrimport osimport numpyfrom tqdm import tqdmim
2020-05-29 14:44:57
2529
9
原创 6S参数
PRO temp COMPILE_OPT idl2 ;几何条件 Parameterxml='***********.xml' SOZ=90-float(getGF1PMS_Parameters(filename=Parameterxml,nodeName2='SolarZenith')) SOA=float(getGF1PMS_Parameters(filename=Parameterxml,nodeName2='SolarAzimuth')) SAZ=90-float(getGF1
2020-05-26 13:40:11
783
1
原创 IDL调用6S大气校正
sixs为6s源码编译后的exe全路径sixs='*********\sixsV2.1.exe'代码传入的参数中有paraID，表示了影像是否为纯水域（完全在海上）Function sixs_atmos,sixs,temp_path,Parameter COMPILE_OPT idl2 ;Parameter=[SOZ,SOA,SAZ,SAA,MONTH,DAY,paraID,xps,iwave,idatm,rapp] SOZ=Parameter[0];太阳天顶角 SOA=Parame
2020-05-22 15:57:20
2496
11
原创 python投影Himawari-8数据
IDL是可以构建GLT对影像进行校正的，在python里面应该是用gdal.warp()实现的，据传可以编写VRT文件写入经纬度文件的路径就可以实现和GLT校正同样的功能，可惜我太菜了不会搞，就弄个笨办法来用了。输入文件为HSD数据提取的影像或者NC文件中得到的影像数据，具体长这样：import gdalinputfile='*\\HS_H08_20190916_0300_FLDK_B3_R20.dat'outtif='*********.tif'memDs = gdal.Open(inputf
2020-05-13 11:02:20
2351
10
原创 Python 下载哨兵2影像
调用sentinelAPI,输入.shp文件，查询覆盖shp的影像并下载from sentinelsat.sentinel import SentinelAPIimport ogr#两个路径都能用# api = SentinelAPI('*', '*', 'https://scihub.copernicus.eu/apihub/')api = SentinelAPI('*', '*', ...
2020-04-24 15:51:07
1260
2
原创 IDL 读取葵花8（Himawari-8）HSD数据
可以实现直接输入.DAT文件，读取获得定标完成的影像、时间、太阳位置的三个数组inputfile为.DAT文件，例如HS_H08_20170623_0250_B01_FLDK_R10_S0110.DATresolution为HS_H08_20170623_0250_B01_FLDK_R后的数字，代表了分辨率代码会自动对数据进行对应的定标（可见光波段定标为表观反射率，红外定标为亮温）时间提前...
2020-04-01 18:39:17
2545
6
原创 IDL处理葵花8Himawari-8标准HSD数据——制作大气校正数据集（太阳角度数据集）
IDL制作Himawari-8太阳角度数据集（太阳方位角、天顶角）计算卫星方位角和方位角与计算太阳方位角的原理相同，根据《卫星轨道坐标系及坐标转换》（向汝建）https://max.book118.com/html/2015/0312/13226794.shtm一文中的总结，计算方位的步骤为以下几步:计算太阳角度，首先计算儒略日JD，需要从HSD数据中读入observation start ...
2020-03-07 11:57:43
3350
5
原创 IDL 解析葵花8Himawari-8标准数据（HSD），辐射定标、重投影、裁剪
代码解析全云盘FLDK全波段（16个）的数据并进行定标、重投影和裁剪，重投影用的GLT生成方法参考之前的博客PRO Himawari8 COMPILE_OPT idl2 start=systime(1) inputfolder='***\HS_H08_20190916_0300_FLDK' ;获取影像时间 basename=file_basename(inputfolder)...
2019-11-15 19:25:01
5917
8
原创 IDL处理葵花8Himawari-8标准HSD数据——制作大气校正数据集（卫星角度数据）
IDL制作卫星天顶角和方位角数据集（对应NC文件中的SAA、SAZ）GMTED2010.jp2为ENVI自带DEM数据loc.dat为之前生成的经纬度文件，为HSD数据中每个像元对应的经纬度生成的SAA与SAZ需要进行GLT校正，校正文件同样用之前博客中的地理校正文件函数翻译自C++ Himawari预处理代码PRO temp COMPILE_OPT idl2 e=envi(/h...
2019-11-12 11:36:26
3551
7
原创 Himawari-8葵花八HSD数据处理——几何校正
生成像元对应坐标填补边角无效值生成GLT文件cols,rows,LOFF,COFF,LAFC,CAFC参数根据分辨率不同有所区别给出的是1KM分辨率的代码，2km的GLT可以用经纬度文件重采样再生成计算lon,lat的代码参考另一个Himawari8预处理C++资源中的函数PRO Himawari8_LOC COMPILE_OPT idl2 e=envi(/headless)...
2019-11-07 16:34:07
9859
39
原创 python读取Himawari-8葵花8标准数据（HSD）
搬运自https://github.com/smft/Read_Himawari_binary_data，有改编，只测试了全云盘FLDK的数据，还有很多字段没弄明白存放的什么数据###########################************************#"""author : qizhang NJU"""#************************###...
2019-10-23 18:30:25
10449
20
原创 IDL合成landsat8OLI光谱响应函数
Function GetFliterfunction COMPILE_OPT idl2 e=envi(/headless) root_Dir = e.root_dir filterfolder=root_Dir+'resource\filterfuncs\' file1=file_search(filterfolder,'landsat8_oli.sli',/test_regu...
2019-10-22 10:09:23
2107
原创 ENVI对GF-5高光谱数据进行FLAASH大气校正
参照博客：https://www.ixxin.cn/2019/05/10/gf5_ahsi_ac/1.辐射定标ENVI添加中国卫星支持插件，打开GF-5高光谱.xml文件：Metadata中有Gain和Offset（可见光部分的辐射定标参数应该是有问题的）打开辐射定标工具：Radiometrig Calibration，使用FLAASH设置：2.FLAASH大气校正输入文件选择辐...
2019-07-22 12:03:34
5791
20
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人