#shp重投影
import geopandas as gpd
from osgeo import gdal
import pyproj
# Read in the shapefile containing the point
point_gdf = gpd.read_file(r'F:\Solar_Radiance\BSRN_SSR\Station_Point\54102.shp')
# Open the image file
image_ds = gdal.Open('F:\Solar_Radiance\BSRN_SSR\Factor\A20223040520Factor.tif')
# Get the projection information from the image
image_proj_wkt = image_ds.GetProjection()
# Create a pyproj.Proj object from the WKT string获取影像的投影,并将其改写成下面的albers的形式即可
image_proj = pyproj.Proj(image_proj_wkt)
print(image_proj,'image_proj')
# Reproject the point geometry to the same projection as the image
albers = "+proj=aea +lat_1=25 +lat_2=47 +lat_0=0 +lon_0=105 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs"
reproj_gdf = point_gdf.to_crs(albers)
# Write the reprojected point to a new shapefile
reproj_gdf.to_file(r'F:\Solar_Radiance\BSRN_SSR\Station_Point\54102new.shp')
基于Python依据遥感影像的投影方式对矢量进行重投影
最新推荐文章于 2024-04-19 03:48:37 发布