地理坐标系一般是指由经度、纬度和相对高度组成的坐标系，能够标示地球上的任何一个位置。经度和纬度常合称为经纬度，把球面上的经纬度显示在平面地图上需要采用某种地图投影
地理坐标系由于地面高低的不同和地球形状的不正规，天体测量所得的信息不足用于明确地计算地理位置。一般采用某个标准的大地测量系统所规定的地理坐标系统，目前最常用的标准是WGS84，它被卫星导航系统如美国GPS、中国的北斗导航系统使用。
投影坐标系
投影坐标系最出名的墨卡托投影，这里不过多讨论。
需求现在有一份栅格影像tif文件，需要把这份栅格文件去除投影坐标，变为地理坐标。准确地说是，让这景tif的坐标系从GCS2000 变为WGS84 。
代码#!/usr/bin/env python
# -*- coding: utf-8 -*-
# @Time : 2023/12/5 14:40
# @File : reprojection.py
# 栅格数据可以重投影，但比矢量数据投影更复杂。
# 栅格数据需要处理栅格数据中像元会弯曲和移动的事实，一对一的映射并不存在
# 通常用最近邻域插值法、双线性插值和三次卷积插值法进行 ...
需求现在有一份SHP文件，需要把这份矢量文件制作为深度学习的标签。这时候第一个步骤是：把矢量栅格化。
代码#!/usr/bin/env python
# -*- coding: utf-8 -*-
# @Time : 2023/11/27 10:18
# @File : shp_clip_img.py
from osgeo import gdal
import os
# 根据shp 裁剪影像
def clip_image_by_shp(input_file, output_file, clip_shape, nodata=0):
'''
@todo 根据shp 裁剪影像
@param shp_file: shp file
@param img_file: image file
@param outpath: output path
@return:
'''
outpath = os.path.dirname(output_file)
if os.path.exists(outpath):
...
python
未读为什么要使用多个版本的python？前几天想研究insar,发现了有个第三方库pygmtsar。在pip insatll pygmtsar 完成后，测试一下import pygmtsar 语句，竟然报错。
搜索一遍后，发现是 distributed这个库不支持python3.9.1。
简单的解决办法是把python的版本升级。但是升级是不可能升级的，升级了其它库是不是得重新装。
多版本python的windows10 环境于是不如多下载一个python版本。
安装好后，把python.exe复制粘贴，把复制的文件改名为python311.exe。
把D:\py311 添加到系统环境变量。
把D:\py311\Lib\site-packages 添加到系统环境变量。
在终端/cmd命令行输入
python311
测试完成。
进入D:\py311\Scripts，把pip.exe改名为pip311.exe
为什么会用虚拟环境？作为一个使用python开发过大型项目的同学，一定会碰到python库的版本依赖。
例如 直接依赖：
A项目：使用 django 2.2
B项目：使用 ...
GADM: 全球行政区划数据库GADM主页：https://gadm.org/
GADM，全称Database of Global Administrative Areas，是一个高精度的全球行政区划数据库。其包含了全球所有国家和地区的国界、省界、市界、区界等多个级别的行政区划边界数据。
警告GADM提供的中国国界数据不符合中国的领土主张，省界、市界、区界等数据也不一定是最新的版本。在正式刊物中发表使用此类数据的图件时需格外谨慎。
数据下载GADM提供了两种下载方式：
下载全球所有国家和地区的所有数据 https://gadm.org/download_world.html
按国家下载 https://gadm.org/download_country.html
由于全球数据量巨大，建议根据需要按照国家下载数据。
需要说明的是，GADM 中对country 的定义为“any entity with an ISO country code”。因而如果想要下载完整的中国数据，实际上需要下载China、Hong Kong、Macao和Taiwan 四个数据。
由于GADM提供的中国国界数 ...
1问题
java的版本、编译的环境
解决办法
Maven >= v3.2.5, recommended 3.9.2
Java JDK = v11 (Any of the available OpenJDK distributions should work).
编译的环境是window 10
2问题
netcdfAll-5.3.3.jar 太大了，mvn下载时失败。
解决办法
单独从web下载netcdfAll-5.3.3.jar
然后在命令行输入
mvn install:install-file -Dfile=文件的绝对路径
该解决办法失败
3问题
mvn clean install
通不过测试TEST。
使用Maven打包的时候，单元测试有可能导致打包失败，而且单元测试比较费时间，会减慢打包速度。一般需要在打包时跳过单元测试。mvn clean install -Dmaven.test.skip=true完成对snap-engine的编译解决办法
</table
kml格式是一种在GOOGLE地图浏览器中用于交换地理数据的文件。它基于XML，并定义了一些TAG用于规定地理数据的显示方式。KML可用于定义的地理特征包括地点、描述、叠层、路径和多边形等。
实际工作中我们拿到的矢量数据格式大部分是shp格式的文件，如果需要在谷歌软件或者其他第三方软件中使用的时候，有时候需要转换为kml格式，便于使用。在《QGIS中如何将kml转换为shp文件》一文中介绍了如果将kml文件转化为shp文件，参考该文章也很容易实现shp转kml格式，当然借助ARCGIS强大的工具箱也可以实现该方法。
本文介绍基于ogr2ogr工具的转换方法。在进行各类GIS文件格式转换的时候，使用ogr2ogr工具可以实现dxf、shp、kml等很多种格式之间的转换。
一、ogr2ogr工具安装与启动。工具的安装与启动可以参考文章《借助ogr2ogr工具实现SHP文件编码转换》、《GDAL的安装与使用》。
二、shp文件转kml格式。在命令行界面输入“ogr2ogr -f “KML” +转换后kml格式文件 +需要转换的shp文件。
为什么要转为KML？
arcgis
未读arcgis闪退再一次出现了正在画矢量时，arcgis软件闪退，导致shp文件不能打开，这次详细的记录下来，解决方法。
普通的修复方法在ENVI软件打开对应的shp文件，把shp拖拽到envi中。
点击Load Data.
在左侧的Later Manager中右击对应的文件，选择Matadata Viewer,查看shp有多少个矢量。
上图显示，该文件有180个矢量。
在 dbf viewer 2000 软件中打开 shp同一个文件夹下对应的dbf文件
在 dbf viewer 2000 发现dbf文件中的记录只有170个。与之前在envi中看到的180个不相等。这时候只要添加10个空记录（180-170），就能修复好shp文件。
成功的在arcgis打开修复好的shp文件。
最后还有一步就是，别忘了重新给shp赋属性。
总结解决办法就是，令shp文件的矢量个数与dbf文件的记录个数相等。
Landsat 9从USGS下载Landsat9影像，Landsat9影像和国产影像不一样的地方在于：
1.Landsat9影像已做好了几何校正，无rpc文件。
2.Landsat9影像以单波段的形式，将每个波段保存为单独的文件。
Landsat9 的真彩色RGBA众所周知，Landsat9的第2波段是蓝色，Blue,简写B
Landsat9的第3波段是绿色，Green,简写G
Landsat9的第4波段是红色，Red,简写R
接下来，要写一个函数，实现这个功能：
输入Landsat9解压后的文件夹，得到Landsat9的第2、3、4波段的文件名
这个函数的作用是，有利于自动化，每景影像的名字是唯一的，如果每次都要自己手动指定Landsat9的第2、3、4波段的文件名是费时费力的。
注意分析Landsat9的tif名字的格式，提炼出共同的命名规则。
LC09_L2SR_018113_20230209_20230310_02_T2_SR_B1.TIF 的特点是
B1代表band 1 ，第一波段。
该文件的后缀是大写字母 TIF。
根据这两点可以自动获取第2、3、4波段的名字。
细看以 ...
RGBA格式RGBA是一种色彩空间的模型，由RGB色彩空间和Alpha通道组成。RGBA代表红、绿、蓝和Alpha通道。 RGB是RGBA采用的颜色，它可以属于任何一种RGB色彩空间。而艾德文·卡特姆和匠白光在1971至1972年间提出的Alpha参数使得Alpha渲染和Alpha合成变得可能。
RGBA，更进一步的展开简单的来说，RGBA格式比RGB多了一个Alpha通道，并且通过这个通道来控制整张图片的透明度，alpha通道一般用作不透明度参数。以我们的公众号的头像为例，如下图。
观察上图，可以发现，RGBA图片和RGB的最大区别在于，RGBA的背景值是透明的。RGB格式下，我们的公众号头像如下：
RGBA格式的好处在于，可以利用背景值为透明的特点，与其他图片进行相加。通俗易讲的来说，RGBA格式是在‘抠图’的必不可少的一步。换个说法，RGBA格式下的图片的背景值，是虚的，而RGB格式的图片的背景值，是实的。
RGBA格式在遥感中有什么应用呢？
最为常见的一个应用是在webgis中：缩略图的可视化。比如你在欧空局、USGS下载数据时，点击数据浏览，看看数据长什么样子，这时候，你 ...
HDF5、NetCDF默认HDF5、NetCDF这两类文件是相似，读取文件的步骤也是相似的。
HDF5的文件后缀一般是.h5，NetCDF的文件后缀一般是.nc
以下用nc文件代指HDF、NetCDF文件。
对程序员来说，nc文件和zip、jpeg、bmp文件格式类似，都是一种文件格式的标准。netcdf文件开始的目的是用于存储气象科学中的数据，已经成为许多数据采集软件的生成文件的格式。利用NetCDF可以对网格数据进行高效地存储、管理、获取和分发等操作。
由于其灵活性，能够传输海量的面向阵列（array-oriented）数据，广泛应用于大气科学、水文、海洋学、环境模拟、地球物理等诸多领域，例如，NCEP(美国国家环境预报中心)发布的再分析资料，NOAA的CDC(气候数据中心)发布的海洋与大气综合数据集(COADS)均采用NetCDF作为标准。
从数学上来说，netcdf存储的数据就是一个多自变量的单值函数。用公式来说就是f(x,y,z,…)=value,函数的自变量x,y,z等在netcdf中叫做维(dimension)或坐标轴(axis),函数值value在netcdf中叫做变量 ...
