微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何使用 GDAL/OGR 栅格化在 jupyter notebook 中,EPSG 中矢量文件的 crs:4326

如何解决如何使用 GDAL/OGR 栅格化在 jupyter notebook 中,EPSG 中矢量文件的 crs:4326

我试图在 python 中光栅化矢量文件

我使用了这个代码

CREATE TABLE #TempVOUCHER_LINE (BUSInesS_UNIT varchar(10),VOUCHER_ID varchar(16),VOUCHER_LINE_NUM varchar (4),MERCHANdisE_AMT decimal (8,3)) --DROP TABLE #TempdisTRIB_LINE CREATE TABLE #TempdisTRIB_LINE (BUSInesS_UNIT varchar(10),VOUCHER_LINE_NUM varchar (4)) --DROP TABLE #TempVOUCHER CREATE TABLE #TempVOUCHER (BUSInesS_UNIT varchar(10),FREIGHT_AMT decimal (7,3),SALETX_AMT decimal (7,INVOICE_DT datetime,ACCOUNTING_DT datetime,CLOSE_STATUS varchar(4) ) --DROP TABLE #TempPYMNT_VCHR_XREF CREATE TABLE #TempPYMNT_VCHR_XREF (BUSInesS_UNIT varchar(10),PYMNT_ID varchar (10)) --DROP TABLE #TempVOUCHER_LINE INSERT INTO #TempVOUCHER_LINE (BUSInesS_UNIT,VOUCHER_ID,VOUCHER_LINE_NUM,MERCHANdisE_AMT) VALUES ('11000','00720667','1',17480.820),('11000','00721196','5',4680.200),'00721021',82.64),'2',542.12),'3',126.89),'4',1022.03) --DROP TABLE #TempdisTRIB_LINE INSERT INTO #TempdisTRIB_LINE (BUSInesS_UNIT,VOUCHER_LINE_NUM) VALUES ('11000','1'),'5'),'2'),'3'),'4') --SELECT * FROM #TempVOUCHER --DROP TABLE #TempVOUCHER INSERT INTO #TempVOUCHER (BUSInesS_UNIT,FREIGHT_AMT,SALETX_AMT,INVOICE_DT,ACCOUNTING_DT,CLOSE_STATUS) VALUES ('11000',59.50,0.00,'05-03-2019','05-14-2014','A'),55.50,512.890,'08-10-2020','08-12-2020',129.01,'12-13-2019','12-16-2019','A') --DROP TABLE #TempPYMNT_VCHR_XREF INSERT INTO #TempPYMNT_VCHR_XREF (BUSInesS_UNIT,PYMNT_ID) VALUES ('11000',''),'') ; --sql QUERY: WITH CTE AS (SELECT CONCAT(A.BUSInesS_UNIT,A.VOUCHER_ID) AS INVOICE_ID,A.VOUCHER_LINE_NUM,A.MERCHANdisE_AMT,C.FREIGHT_AMT,C.SALETX_AMT FROM #TempVOUCHER_LINE A LEFT OUTER JOIN #TempdisTRIB_LINE B ON B.BUSInesS_UNIT = A.BUSInesS_UNIT AND B.VOUCHER_ID = A.VOUCHER_ID AND B.VOUCHER_LINE_NUM = A.VOUCHER_LINE_NUM LEFT OUTER JOIN #TempVOUCHER C ON C.BUSInesS_UNIT = A.BUSInesS_UNIT AND C.VOUCHER_ID = A.VOUCHER_ID LEFT OUTER JOIN #TempPYMNT_VCHR_XREF I ON I.BUSInesS_UNIT = A.BUSInesS_UNIT AND I.VOUCHER_ID = A.VOUCHER_ID WHERE C.INVOICE_DT > '01-03-2019' AND C.ACCOUNTING_DT < '06-01-2021' AND I.PYMNT_ID = '' AND C.CLOSE_STATUS <> 'C' AND A.VOUCHER_ID IN ('00720667','00721171','00326705')) SELECT C.INVOICE_ID,C.VOUCHER_LINE_NUM + ROW_NUMBER() OVER (PARTITION BY C.INVOICE_ID,VOUCHER_LINE_NUM ORDER BY VOUCHER_LINE_NUM) - 1 AS VOUCHER_LINE_NUM,V.LineType,V.[Value] FROM CTE C CROSS APPLY (VALUES('Merchandise',C.MERCHANdisE_AMT),('Freight',C.FREIGHT_AMT),('Tax',C.SALETX_AMT))V(LineType,Value) WHERE V.[Value] <> 0;

(矢量数据的crs在epsg中:4326)

  1. 在 jupyter 中绘制矢量文件时,它显示了坐标 X轴:74~78(经度值) Y轴:35到38(纬度值)

使用此代码光栅化后 X 轴:0 到 250 Y轴:250到0

即它的 crs 值丢失了(我猜)

  1. 它创建了一个光栅文件,但像素分辨率(0.015625 度,0.01171875 度)已经足够大了。

如何在此处调节像素大小。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。