如何解决如何使用 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)
使用此代码光栅化后 X 轴:0 到 250 Y轴:250到0
即它的 crs 值丢失了(我猜)
如何在此处调节像素大小。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。