From bee076504b5f5d10fd683fdabc44e0691c5aa45e Mon Sep 17 00:00:00 2001
From: xuhengjie <xuhengjie@mercator.cn>
Date: Mon, 20 Jul 2026 18:45:38 +0800
Subject: [PATCH] =?UTF-8?q?=E5=88=9D=E5=A7=8B=E5=8C=96=20v1.0.0=EF=BC=9A?=
 =?UTF-8?q?=E5=A2=9E=E5=87=8F=E6=8C=82=E9=92=A9=E5=90=8E=E5=A4=87=E8=B5=84?=
 =?UTF-8?q?=E6=BA=90=E8=B0=83=E6=9F=A5=E6=95=B0=E6=8D=AE=E5=BA=93=E7=94=9F?=
 =?UTF-8?q?=E6=88=90=E5=A5=97=E4=BB=B6?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 scripts/__pycache__/process.cpython-312.pyc | Bin 0 -> 21557 bytes
 scripts/process.py                          | 395 ++++++++++++++++++++
 setup_repo.sh                               |  37 ++
 suite.tgz                                   | Bin 0 -> 14918 bytes
 workflow.yaml                               | 105 ++++++
 5 files changed, 537 insertions(+)
 create mode 100644 scripts/__pycache__/process.cpython-312.pyc
 create mode 100644 scripts/process.py
 create mode 100644 setup_repo.sh
 create mode 100644 suite.tgz
 create mode 100644 workflow.yaml

diff --git a/scripts/__pycache__/process.cpython-312.pyc b/scripts/__pycache__/process.cpython-312.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..f4ce522ebe68105605dbb0cd2d03c6e50df224e8
GIT binary patch
literal 21557
zcmdsfdvsG({_jb0@=lthX<GUYrLPo9Z7Ht`h!#qD6|g)jnrVB|($XgNBmr7eEsjd6
z)utobARs|;jEWdXXRP9bxp(dw?z(ebCnQel`E}NEIZX?9&AP)V^PAtof9`jmlXH^B
zR9rLu{`lPu-DjV3_IvO1eSE&(o&Tj)D=-Le|NX)4k~J9iXLQ2_JeK%YA;7S4jKEqj
zLO=+61T6yQDr^xlS5b?Ix#BH2a}~FUnX9Bl0#{Lwv{%+56QHtqkGxmjB4^UV9z~0i
z$){>jGuN~h4Iw5Zgp`n-5w&QAm_=vNc9(G7iC^D}CEtR6q0pl5#0dE-f|d*-lTbV%
zXfeQ53D-<Q^#s<EMW`WOK%_yOO=uv_A+!+Z5;}<Uh;)eai7Y}7B?^cPxGS_63B%!X
z3?pUGi$bnF0oE_<7n&A0ABBb$Gl7Zdi%&$K`z-p(%h3~~Gb4}R_{AGh&vVfWC$9hc
ziJ6Z)*8@MFIrrMk3xSzuFGW9ke)i<K8Mnvy_0UOU^uiOfo~N!q@zTthiwece@axyl
zkIeYbM?W5|Qy7hwV-B*l-C?e|t9rGO2{WUf=;({+KKkMN(-Xr`bM*2@*DrlKd(t00
z6PQz`X8GJQHOrw)bksX@`hB)TNhNJZEM!k>zq!lmu#k33yTceA^+lidMcq%$yo{Qi
z{^-Ju(@*{NWB2sL1yt496VnqPB~}tW{prk$V}Jd481;i}wOh=j#qJ=x+8tdstKE2G
z$P<0{@tNm;dE?ZD8>c)`_iwKM^7!ogKFHnC)dP*1TS-eR+iByCv#-p){?7H!PRv}s
zU|ejBe)4Yg=c6;eiRh)5XC|JEKKsJ;PoIyT_(*XBt$UM6K;t`G9f#m5BU&96M^~@K
zsew1qr_Rir9gALiGJ5VKm>{PLJ^SFb>4}%v8Dk#YIB{b3+-K|qlZeLeT6tFu#4GNq
zX5t!?l*Tt~yKg%ax1e}?JuNxV9KYVzLW?(SI{;ZE=%Kl~J+>oJ?=bZJh~V$r3FiB8
z@b9{X==Iz{2;Pdl9~8I*-Q2gAkckw!1kD({lO{ZjB^~^uZW*8C6I{&#ETIp~ch6^G
z5`>TtUBust^=?p@kQ+^4QgEbk54yEnDt_hK92B`k-TL@*F69z(IRt1-2JU(M>Jla7
zz{u4Jy~&K{;!-Y*5MPw=eaBt+KiMjt&??NgN$NzcU4QNa_!O#cJo)VGCzoc%UyFV=
z9zEr!#r<~kK<7YNNqMUswtKI|WH&xgQr=@ZV(GEAqT2_JCFQ*qyS=p&Zf@&J4pf%(
zR+bRPlE%7{?R6#1uyv(q>vr@yZi}6cj>=VUz-G2s+igUbwbP`cr9HOJPDs-7PK#q3
zq%9;ZJJf0=df;=VWp;RoTBGsJO$|F}d4rAYZFRs`MXQ?~@LlybvRP@7z28nt><+@_
zbkHg&7SE|@B?vB?vZ1TT61y>pX&JO`?t@Q^mY{w>wIb*#Ewmq^g&l|NVwf7{%e~!H
zb;#Chsj~G!*LzxzRxRaIM{VR`dtYn2rK;a%?X*MJD+esr&h9QqbUIr*ZO$rtJK5Fe
zuvhhwwsx4SrG5Rh0!C|Qn<Rk{V@CM1Uw~)`n^tE;)CECxf%j-gT{hG(DOI0%?D5AY
zH3qkO5?7qqaeT+gP2=0fwuSM9Qwm)~ksDOxQh6Jx9c^Jn`%wLqTpN)W1my*1?B@nf
z4TR(sL-$SM5>(DpAI29z9TD6R#0{RdFrNLTR5OJdGzQg1--5uVuzHPKI4M(%YKApl
zQCOBYDNBpUvV*d0Z&6s5KdICWZT}YfiM1zw<vcAIf;Z!s3+v{o0m{pUI1M1PT>>)8
zfv5))Qqcn!d^Z<`+>UR??v$I`*RZEW&UC2u#t@>0v*%v9abooPJC|X9-xe?J>^s~E
zJ5UTP!9oI6gsn)}k9NTZ5Owwuv>dh+VFtQwF$rlLZ49y!s<NZ88jWNbTxcbe=LoP8
z8%fr{lWLUXE{KM(ucXQogU1Iw`@+)v^CeT7^wGiLLGNzAA*`vQWK~njw9&Q0YrPpE
zWdVg3kR~WJho&Ox#?aImNHFsV*W`dGC?Uv($8Z&*>j3-_!T}|Q0afE?)dlp9d4r=j
z-24z?^bF{hfDp&(knnW~4j+VfKwU^NA$2jOtvDfLQ-fj$`XHFV$qa)Mmn0R%Cgjk%
z;xtYuU82(%p>l~&V`n8?A3#D<Gf#xDT4~M(sB*^lA;=M|G;w|8=h4&ePk$Vky)<sz
zvwOR7=IIxJRRWLZg_7xsq3iC8(KCMI#+}WZ>Xy~4s9gmtnnVkAgNfcpDG7?uR=j~`
zrEC9V5Pb}hGXrjoMz#sI_Bvya&AN6#w6vyUK)S>@V6$4*q6Aa_fS9?Ma9ZWCIa+(n
zZ8j%t99jfo1@ryVc$c-q1}s(7Np#SX*lR_*1&JN-p<8L0)6r{2wbHUKyV+_ZVb9@~
zK6@v*2>N12T#q3&+&0sKkkJLS&SFLd+etflA7ojD@}7kV_L^E7QRfHM`9LDl7kG7J
zCU1vt*Qu^xdYQjEs9Q3$1@@n3i%%cSDE142dQ(VK;olY1EFEeDTH<LP?)R?nEj+b4
zs6nqarlH0!3(MTgJ*%!MbNzb)`-7Ehhc-;9vz|E!4FKKQdVDLDUg{eNETT$R1Ue@U
zO|)D&dew2YgVOE^<9jF7>7H6ose2uTXR{j{Edv*L0`Vq1nQvp04L~@!5((E-0tI<j
zd^5WE%{>T0M9S?ba0y`pCNgELTb|HFEE#9pAc|nHIfc$;@M`w_rRmSk#zfsSes;q|
z&w6HjCuZNfJU#Ku^ySy4KmF|brSq_3SsBk+2<41Mkj_99=Cc|r*Tpt-)EzJa0}{c)
z(rVZcofcas+1hugtDRYA@@^<jOB^<|@+J`p$_vJ96j~bF4tBJoVaKpUNoTRzoutLw
zYeAAJvt1hD{Ra56Z-NMsCdG-Pk013c8(%rLGNLcPrZ1mLFY+x5rI)#-K;%c2!%8Z1
zd7v|JfXcWhEL%Sb)N5dPz*`p56j8FGNoCs5mhUetPq7+diN~?Hq_zQj@GpoSxrDH&
ziCY1BgTfI3DHutI<>F}y`y9o@_m)cz%i&^ZU8?1Epk6S6qcExD#BQL1$$_V1NRH!c
ziRVHxV-D7>jz8g2aruyI6Vc6X9rr$d<=YUu!~}jA)`T>}XmCVzQ%P}rehJK-gb<^(
zcZp*8=kFsTX-NS}OPJDdWm?Nh$^f@qYuzYAQsCYVN?oWg(hgC(PynCB$p6I=kIzra
z4!4B#qKuo5RI3iV;{UL_BED0OV6YQZ|DUb!(f@KQocc=yJ`Vg_fw(oHJS*nb7POkV
zJK`@mNFao;3+oUyV~a3H4!1XoF%lPG2QZItM1B-ICOUu}6{LJG7)j&an^ewGD6c>a
z@%o83qfh&1ULBhG^)rZ?zI%E0;}>Ur$D?nL&s_R-^uyn<+B9g*iAp=j3P|6rLh3Gx
zYCt&n^U&jLjh~&2o_^U_&Z_4Y!VD~XiP33?R7aPkhgf0+?Zjek2Nv(Jn+Z_#O-wl_
z(gT5{2hLu9=LNR;I%CrfpnDi2chE{G0|l8%Xoa~4R1@Y2EwQzAgS?28p)SxLoP>qZ
zx{O9ZMP}c7Iy&+pXkJ{6PQ;g%7^7pKOi#RR{5z`NICJ*ANkc1<x^8aRym8xxW?I^A
z>v8s4?J-W#YIpSaSwIbPkhB=V6MZ0(*t@K>yne@qjpl8en|H%YXD^Mn+d%)29BJ)w
zTI{rnxyHsyOM%~@@zU~*$J#A@2;iXAc1sWR-AthR?Iu|QD-yFK#i0&6jrX?N57T01
z)@TVpB+ecOEnz3oj;Iz|OcYp#RX|uAs9H2fF`T3w;1mbxWN;OdfD+r`&t3_<V+gyY
z!qW4<!8BrnTRf@MdPV1yr<6hEqDfT-&$~)}PG2XLu{12Jx+zr34U_szk8`XpqAw2W
zi!U^X^rncuGN`YN=$8fc%R>6)?hTXrtnpQ2tGxGlcX`&tG9C!&tKdz>g7L<&#^<-V
zH%@9YJP$wJKbf89-F+(iT*0Y=P`1e<nJlmLNIZi<ZOLTWV*iHo1Ft`N{?R~N;9;t4
zb);;4uxvfNb_KP?li7tnRWQ3Ul3f?ft_x+aiDYjGW^V{(H$Wv>`Cg~5gUVWRO}pf#
z64Pei)WINc7Gvp#h^{E8EAp?RbVVUuZA7;ws9O`Cf|FZ=%0d|aO`+<)$WXPJ<H|85
zRZtt)6WBDd`pTXwo37SUc{{_}U2~^Fjiy0^rlAN^X+0Spl94L^<JMw$@eL-B8)z7}
z2(I7VDo%cM+m0aBOZB2gT==`f`rO70;Z<DTs1sk+3Q(LOZ(Jz8Y7{`6T7{*cIKaVz
z)Uqfrt3?3DFzQz~uWCV&Sm0P@DzS!p2T&I<Q3nrnAyt=CbO#Y*O~K^5Qvh12#1ev+
zaZxWB=8xo9kuw)Mm?)qEg^Ygq-pq?1&c6E|5Jy^yU`lgu4Ot5rO;U0fTxhY?=CHQX
zqE>4^*@zy<Ku$$6W<QOi{Ir-vijYYdrwVg)WqEsRPdfvEl08toAO7qu5Fv3?>!}~#
zI<_?^E%ep<<>#9s#jAqFt3t(fk>Yj1;&tI-B!RjWQ;J+~y)QGQC?BeyG%OfDICd~%
zC=D7)BZehG!xI0lkfCa*-Yq}glpvud%9Lp0&_aECKg1xR^2kISfg1$LEY)HXg44We
z0UR%N9-#{d-Y9ZmpkO4TYpHY3sWDN70kk8IO-klzOPsK2cx@V5mz+&N6CJ!_13X&H
zZLEVZ`a$WSY*6l!9#2cCN6PVZLVCCySQ!CrD`CEp0wHtoC@HikpVJ}&rJ0ZlD<u@2
z0+1ji-3XOQ3Kx}c#Xh1zg-ekVClwN^i)wBZfbsG&Z1;ks%HWnr<3<DPCGSR?J}Gb?
zAxNStFzg_3mqF#A%8|!q>qd%OQgA5=4cCJ~KrT}56|GAN^P+R9ncbj*IQ=X?w*qWr
z6^=O;auvIz@s$DTT1w~{`Fs)gK7Qr;JD4_4yBVq05?VL7(uhnzV6x^&ZPFhz$Mdwb
z;Co6m!f&#jOQ5%`4JSzivpXz()y66#y4Dz@PreI)9MGEEYu5e?!HGZ}wz=EFIwOf%
z=fo+-7{!8p(Muy>kE$~|8DxdE$uMG?vkD4)eRjk+bLP$KKmTO<@^h?W57O54cb<y6
z-%1<=D{C5^2*CgA$0OI@`2b0x=U$6RpPpyQE_ee9IwF@$n%n>>8sD^Y>m~piB?Q?8
z7ASHhJY+k5=FQo2uQ5rZ@#`~B8gC1X<PKE017$(P3Peng8%TAX@$X2RM8e+^<IK=y
z7Feh=-o}>zNW}_<z5nf3e@6g5Ti5L6=Vx9%!`hG#)L=?yz$r#VWS~k~)MK$SvK1|c
zUIBo)7YdVTMrhS0i)}-zqqTtqWS<=ck_QnON+L}U(5hoy4qDvhu=LvL^!go5&Aad4
zQ@?xjjwUmbF`EJV0?vTwJh=}d5Ll%|8-uPhAS$f{Xui+sV4oDB8UcoFv-B{y2FpO&
zP~nGA)Xu!O^!7P`)gZANpc1Qt#*x{9R<Ix(;y`vpdO)g746KByLc7i!vo?}$$o3Tc
z*<S$qk3jY;uPmfq$U|G66`t}io)-hUJf&WPSL!8wWxi}myEu$5;a;%WrWFQHr*BVK
zu^5;VhmDnA)cc#?THsd%$|!AJ7{7ZuuWn-56-g*>QxMM{YH)A#tN<{7=+QHJ@7~k7
ze!;bD6B1rk!zwDP;dlL44pX~2C_`shc4%6wk7$iSt<hiM&%2_ab~Z<L?hEeR7uN0v
z`B$bMRS&B@hyA5)by!wCB~y;dhGibzi1MZc%gjS8^7)6p!Gz+K!*Vz7ZlBUDjA%-N
zni3!3ZwzUcyYZOJ?CJ3C_BQwy`ZoAh_|yIM{^bFj%2@t|Y{gA6l)fdy(y}~9yhlRn
zQkW=RmRAk}X+&2V)Rl&G<wIMzNvSx~?8`j;fOk7(s+*|0LQr*EL9)~C4CA}Lg4t<0
z-sHL4cVwt3j4zGPmOt08rnD=<_{vFTCWU9R*b)wF#0B)_BJS|53gUS<H^N*{o}`e(
zxl`xiMJ3=BsdLQ3t-BEkN(u=eWCBi}GH2WtDw(HUQ34O=b#9>PI8f#kf?$&5k$NoR
zggO&t2OK!_o$E~5Zt9YUf!L9`@Gb>vmb(8O&mNS-C8BQL{u-ZXmlz>h48|1OCx~~Z
z!rUN_oUc{t;4QSkcjHnPq*Lx4MZ!$+92{m}34B_`@aYuxT#io@YT&tPNj#UASrWz~
zBs3SboK!LI4%DUOx#<i;O6kc!TA>YR8AOS%|N8MuU;r~RJ{MoR-+%--Qn0QMP9*FB
z&H|nXPrg<S3==Jzk!6@mjP>_7ubq>lrp`#Kcy9_o2zI|Aq(z5pN9nY~79hAtK`<k%
z0JtAckhZ=UY{E9iawy>12>xSO9_d7>BPi-a(ZeV@3Xw^bz`PKLB9V=aJcc5~W61`H
zXlcyB1DWN(>I;4zgymSA36YOLCX<YobXkG*RL8D0<Z)C1;z2B^jspY)ZV?dIkcbDR
z04!K`5^Ifo5i<P}W!&;(*a@Xw9meZ8Y8hi8F<W3}NNEH+%bYtd=i)Q%zT(qe-Us~Z
z0P(@0cMk<xDDC<%e(xPK03(p8P8>Ua%u|a1j%(7~NsWH=k>N+YwdYoyS{2e1yYX8R
zEKN81!0-d!46nmqa1{@$w*#41>m%yIpt{hvn^G5s)D=SwQ|e5f3m|&VSb~Z!Eb|?C
z^=@A!Rki*~=2dy{-uol>?hW3%H+1iQ>cE54`iH28+9+*%7$;C~|5+Y@r~tzQ0zh;B
z6L>(PewUgW%;f<*e@aE!Q<{``YHv3Iy+I#aQXs%+HP<-h;3;P+Do7;idPqWKLMrVl
zktfYV>u?88%RnQ6&5Q!jI`JI+j;CKhf#dX3+=&eWJqAq*kUC@|NM%EC@H&`-w-<Et
z+7{a$iXA^Zv#GQ)flJEKPp;?5SJ1n5fUQDeuYmC{jJL+6z{(ACPskES3p{EO_db3F
zd$#<df@=@xHRwL`%&9VAP7^t;YECO?-oUsxD4VAi^&Mse!CXnf#oJ!y?x&IEFbCpp
z#1xVO$F&CK4iopvB~MsmxeNN~5j+UAa!~fK^k`5a#~eSrTlwA(DqYI%3a&r#D`@m7
zp2p-2Wp`yfCzncHg4-p7i7OSC7gR1ek;d?O#CDPbKV#~I(IC#06krwsH6XMXb#HO9
ziX3`_%*jas@}%=UlE?9lq;*gw>?*#m10cBpGNAr3J!bN1=Gkp}mjb?n3|AV&2A76K
z4pOzqy#EpN!_=_gtdbj<0Cv$x3G74pkNUb&)8f2&dX)cT#!nQ0WL7vwGJ|gzY*7gD
zWCGZ#iUzgwP>O73hH_lmxzhFljL03-KB|S$on%u-1-u>0=tLORd@iD_sD`LwqT}f7
z=}Xa1UYMQ$gOC>-TtHR=c^ZpJQLL243TN<U2fQ)i1LBIM_CtLGk_yJZX+REFi>=#k
zv!X}n8Y}(vOTNEOj3VFzdCoWy1u+8bF*E!rxFfKhQZp}}n7;fQ1a*Ay9<r=4eywrI
zj4@6DKruQt`&_{2M0^2icr&s8Yz=h-5~Pw>8W|yvRuEks9TpP2=PbL)XP^}M97M5h
zvET*Ul<@%u19`B{IQ#l#aH%yq5tfTEkaC23+||K|rvQ7Hz(h9u#`Gtrj861@p(jox
z=A%}*aWb9fYBYk0@W#ngV6uyO{~+{*ZIr``SgDj3NKLtM36vJv>|{064*-ar=pspk
zWiz-4=|ZU_0TrPjv;@Eh@Gk?mMDX^p5FlUn*xG2k+XiMi>Cslw3Vx4_fJ)2Sz^@oV
z8{~de8G_p5GA@IrkRJ4G3`GYaGU;dqBhs?oi<$(NWU%0oN%AcNuF$%~ObFHi|0U=l
zEhIY_p_e=l71J^<7kDvQI&7qcR<c(H#;2vEr5A)?f>y;6>|=q&j%YiY{5gSzVzO8$
z=6lZ5<l9i@xA12-fQ)xXREE-4gz?3k?8}MK>d}?MD;c-C!mu(9sMH3`RP9E}&=8hw
zN<c`a73$H-;YzO}q*#=I-5_b!{^}}Ug?|I3tq$WgQwqJO-V0_2<4`>a)7iNcZkW`p
z2J4e}zrZb;lx0U`#-Pl2A@d^sj{FVz`$Zp=y<7Idl6RL(JO*%0sAhMld{1cM-Y;bP
zz)}Te{wBdx`3QGO8&0FLP5vstN~k(BmG!f*taTC~9Pq#KZ1qdrO0Zu}YH}!9PE6MJ
zn9kHw#_CXB4VAl$TDFDKZVltxCesR0XO<&xO~!f0({bkhbNf#1i{w`Z^DF&rq5NtJ
ziexPNLRO2qupEW~^8Wte{m--@nZCe##0SXE(xHa0(h9%D^k7;QPU#J-d(fhwev$9K
zkiN_<yD8Mk4O3Y;<AY;^zWYO2<!<mF%NlhJySzt`gXuMmX)<@wxqD9C<6rSk-5YhG
z+!Y?#v|)ia^Z5td4U;)WUv@BOaU`cIm{S$XsdjIj0xQk~V-NT;{R;vvYWM!I;Q-__
z<VFl-K|`5;H)SXb8EV`OQwFw+hNbX2Gw+Pux8n38!ORNxrpeq=e`Y`)G}T2+>w>0r
zA=CO`?t1sOuT<%yTZgw&ITij@{))hsiK7!;SDaT5QY{Zr+srpH!OsLmQ$;1dj#n4^
z^8J=T!^G-f?Z!y$=3wpSQ0>;M8>nTws9k$twzBsL{v8vj_X}<o!o=LHz{JYbDop7G
zyMM*2kNA42WgD*Sy}CKLac^W}OK@XLXyZZZ!B%QR8`TCLa-CuP5DRnt$1P;u1auCp
zw$R<C;U+5ouE2tTVxr>8t}FFaZev)r`A@2h?*M|68*bZ?9P)c({XN?%gny`zZ#Rhl
zuq<o4R{X!UVo0|q8oH2b!Kg}Z0M6#9N<7f@FB-a#E2n6V5S!W*GFN5d`<O~)N>EWG
z39t(<HV;Z2Xk1K4G-n~lQvib)Cn}c+xveq+8?QKVWfB~mH1TI#igBusDLo*FOA~s)
z3!{wKBnR;s&23Tu^+5`AB1=^cp352MSA|OsY2{g7Xage_F9;K=RAmVja?<6L4VR2i
zuYr4)ES0+$krp>jK~623lbxAV#vVe)iDZyF9p#RToT*G|oWPsPm`3Qix{&z}>dT1r
zlUMcP69x{)aB%5j7H9AS6)N?dE=X2~!AUTM`bT6FIou2(g~26Y%$I8&My}0n1n?w<
zR0po^P$Gre1oryat|4;ceSl+Hd9MftRj^KANKFo)?tzd1=0mkhbwpqnjL2MSBA-dA
zU65w>ZwmE>C^&*a?n0L;rGlg8dIG9$HJBLV#$SF16NQY^p&nXc;maTRO00SZ5@nds
zokFeQcEcd-IxsLA(Rwqfg?y@rPc1?zaO_J<rS^b&R7{jyEWP7?C_{$TIljBH?$vA$
z6JtWl%><Fo%;nvQZ`m|aK`id%sbO~tiolOgJJ0T1!nFeXSKGbz`)zPFz<TeA8JKm0
z>5g?=&Q#y9<6iEqE1l4y-gB)GCZaM2BbM^JOg~S*s$6=YGIu%Z6KbkRtO?#$L+y*v
zGol7mp=EQ-)B6~RAQg?P{Si(O%enr}`&F-C%;PJ61ov2V$GTSEF-P4UbIctT9ED0!
z1Pcc>Fw;o&N(w}|OEY&rf(mFF)I6#I6|jj-C8&URGgy4m!9Vt}<oA#Nu^J&_4v7Wo
zjO*4KS<^7@vcsrVyICy?ea-;jGp7op!_NYQk321q&SZ3Pk3G-9V{QC#OuJgc>Q|sW
z#k5IA*@RP2T?Dxea?@b;Z6kV)F<wEmu(z8;ouhF$()2K+71H{SE-M^}IAn9SS-Q*v
zUH#5hEBFJEkc~0>I`dGYjH!F(%+Ty*ce0V1Tnrf)J=y&I30Zoq$9{}aGv|2+ruTPH
zQ}3eaR}h)RtQmeQdVq$&S{YC#$h#dbput%kD{5(JM~|%)60BDMnzF<;W}p7e%+PUW
zCKDgS@z<I6{EQmRp4*)M=p*L1W$YB^_-iv~&Na|7bUv)F|5y*4BZB6@CIQ-@6%M*t
zoQNs0^O0C@;_**px_2TD%4@R79kYiq(SV!lg?F#N{~D;VWGiH0ltXsCjWZ*|F%Q|x
z=dOeQ0J#M+Fs=sCCx0`2`NaX;2!XV5Pg&0b@)W8D8DIw1Ec_ROa}`t()4Yv0PCN@B
z1$b%9Uib*vo)Qka!tq*|GdQjaGlmWjzx~>cH!nC5&X0Nr9v0J|pi`?c#|?H~ocm(0
zVe^1<1b++&*|C<vsom>CFTluWetjHt{F%33jRu}zTQo+0F@o}}fIKtq&!#_l9u5vc
zoS4UnHVJAPU>bPKt|NEo_oe7NAA&!~?Av3r!*3v%;Dz`9Jaimto4)+X%!`59m!4pA
z!by<pmjY3@4<HNY(pZJ;euuLmu{u4^LH?v;)9A4AiKk3Oaik#?l)%CQV^A-6+_hTa
zpyFeYCJ&<sp&zu!3I;eCxX*T=W11|e0IGUl4;)M+UqaccQFIRFl>$J(U=6#_otSa#
zr6sL>uww~Yh)%)rZXmP_JU}dH^D(Fc>0`=wg417LTR*r3c6P#PCOF|l+K$@EN8u&e
zjG{V-Oc~@WD1977XwTEK-YzSP7NnumPqw~8HpbtWOmG14qugyMdJ`3q!m!N9rCGsr
zf(1L&mScTT0qbXNrv)9fJobd>ATOd~D#n+wt)H<&(gHIrf)0{T!7F<PL};;-0}E6r
z8K8poM<#=i^*dDNm%tV83{Ipk@-_PE{9vJ56~<R5p%imahbgtzqkm>UDB%ep103h_
z=sgF<3Vf@Ay2U~Dl8CxGsICsFmqpa~1l9Lk0dw3vA@vrp(4}RZw0l=P|41;cXlN4$
z5>)ym-ecf;{PCiRj?0TD@<H3(B*>f0EB4m;_V`!&O#y@g=vLnX(4Y=HI}ADAg{<3#
zzxtFplu-!|D$9>=2R(eeWvs=g0sz5K4bN79W2HXZYXTqYh`u7Iukcrg^eA)f@oiIS
znWH_!J<nJpX(hq55~_4{D6MX2)1)DfEeME1TyLh9Z$Ta)TVpdrr7!hA;_s!jYr^<i
z$e5lV(G>@E#lC&i(#DW(GkBLsb)E;VNsUv6tnqzg`;c49g`(HX&X>Jjb-pUFI8?Nb
zT7Mv9XmKmP&M*ya|4OYNeP|f`UPR(@w_sAK@@&7RT<B}1s=BY0cE^sqjXpN~7*(|X
znr8i^Hb0^*2KQ{=t}nD@H>FVOmIAA;4K%!`oM`=(Cg$56(N+WjcG52S-Tl8m_}hd3
z{h^`FlcnWDo2j&-AigLunO``x*?l;O=S&)lAY~2W`8OrFc+I3t7m?)zWjWqD|DK6D
zz;(j1h8Si;Wv-gYnNU#~8{yea1(uOBzHV%tuM~XM{ou7;#k#5o_Dxu*+6^K7#)y7b
zP`@jrZywslPh)WyFP%&;h@_Xn$-hv#2|4~W3x*|bfqNP6`E#M&zsTS5#^Mkh-6JOi
zRB7D=cmi#Vdjb*OT5qXOh+r`P%JZ(!qFSnGd0_8E{T1<tTdp{$ySD*4lCkRx+5N~h
zs96A=g{1{c@Lb0!O1mY##Ga+T9N0*kC~a*RU(R>9D~uoJZypNcKjUw}v<*j<83S{y
z+OHVbp#Qj4hN%n)?t<BjA!sWC4+q+*th$Ko?rXBUze5Pun%i~+xBpK`c~ds_dA59q
zK=gT0M$;P6=WE1}Zijb#Kxz^acE2r`A+V$i8w7W0<dndMRA-mK^?P01Kjnw>vJT?0
zUmS=D)yuI#LEKj}m6LZ|t%Tk%|2VN<fLRCPKuFvZH<f*u6MYFWqxO|0w21`3r0~Oa
zv1f$28V51!$~)yJupBsc_~1fNb%g|+5n_~AnM;^T?Ezn&xPP(`N{G-II&g;<A;ocI
z5vU|sN(cjfLxr0eIF-zgs5_HS@`+S7>3M3E6TE$wAF<>ot6fQ`xIXi>OXJ??-8p=h
z`NSRCP6e{?z!{;wn8vjSfD;d(fX!UMHZUjy)m%(y7}&&lhqjQ`m=xxjMQwa-_?E!(
zX;;R*pW_hJ<LH=$A*o`xCDI8V{yHYG!>L081BfYgUOkb)jT9`pi@4grehOaa;PIt=
zD+WO#8>H0@V`hK-#`=?uIWxtWRX<tnSt)8)J4*O*ssK`gadLKMW|RgP<pO49dw;Ti
zWT)sy+D}${PKw$!4jz1SrMWb7M`?W(#+b{Dar7taMd~reeIruuz&s+K+kyO^*UqzB
z3*tL_Q0LN-h5tm~c?j1kfbT-*Kq@p7V3yF1*Sd7^J#W<{e0hb;Nls%NhIw4z;6aWA
zPJxuUq(BrQ|4Z<JTyzE}iZ7NVj~ISG3P6-<04tBe_j&*mz&|+k`89v$a|u@=*xC0E
zrUOh`#XX5%iG{B8x$sKrd=8}0FrmAKOU18m#gg1=00rbh&+giIJ}w1(C`J`g&tXo}
zA59;5mc1r`Ii*dsa_7DG-@CmTj$v-wNDD|3*__D?6rpLnEm#XDA;Ddp48R>6{A66-
z2ePUhy?s^eO%)he;cO&hst133Tj$)j+>gY-i`{TwR0dA>c8A$^n3i&`>_bp9<NrDP
z@ki0XSakGpIKWV6O!T!UA4hpq5G5XFU~S01fT!g15TPRgV0LEJv{%ssNlf$kJ>1+t
z?v><Kh#1`*d<Nj$LhM((jLu4U4o(6yC*PSFIbO&8rWYLbh&f0?im4o~Y@G~5!JN8a
zeh!QSF=(mX*^Yh`2Yy-ue!ZvF;k1)~L=8t!Gzn1yOoA%*>jv;EL|t8`d`6*#`LG`8
zBH@H8j<6OmtAn2djRT}Zus))QJEBu+ok!rD0X((Cx!As*E(fbRqum2%-rC@&d*FvS
zWWBA2Eua?L$xl(c2)(Bjd(jcujm#+tH+qVEJ{gNOEoP7-MnPot;(k;>j(*q7($?C3
zm{x${9Dch4sm_dh04?q#0dxXo+u7S^H|c3fUn?B5u#<m41x3(FT5y;a97z~^5o!(b
zGI9t-Xu_FQCHeJZ77p3r(Q#jNLKp1*Xysr+P!egvcFrRrN!?|&T{yDaWhMU|@@<06
zZ(j~_D*9Ow4XBvXj5ES>ic^Z9v}j7LJ6SYdF;?MyFqppB-#oE0tgZ*`QeJfB;VDD*
z_<^wl-iMz*=x+E2oX9L<UQEkXBh8GZIp3?hCeNAD8OC$Qa?aGBYdqBmzfrOkP7e%A
z!2dLZJ$`$y_dX9CZ57E2-1Tti&(r?&z$6%)`^WlydS5G;oX2pFz*9Y$mk(Om!@rP&
zwk9tC4Fv3`b;EUVilHr}T!?<GL^HZ?c;Bdb*zDbVuH{sVUq;nz4drhOA;;D!jc)Yl
z@KK}zvW9Da301u%RJb*y**4S&oe@isK66qyE*+DOYsNG_Nl053DE<Op89yy@6E})8
zCY2h`>TAlpsmulAonxKjwlSNp<Mr<I-GS_h`$DDngfiE=H+`wknaVU$3s;3QSA+lW
zB&c><#<qBm`f<>;uVpL`WWdQ{{$W}uW2JlJlxl%D%U1%n?FAFd|ERk6E49vZ;F`LC
z(}2xmX5Ze|Th6xxWD`Z9lC>ejI%Fni)#KDJp~g(Gv_p-5RMnyS_C0N$%FaEPbt>yz
z;i*D@<~uoW<OH@|DGHf3!ii+}mZ{wQb1P4+jN9|66`QW6hZb)O<!*Oxn^G<EZ3(C*
zI<6f0LbU_cxBr?t|A#Bg%8RXrB5+?Qe`P3Z6?hQnv&ZYk>ZpRH0h~FB{pFIoPSzmb
zgYkRD?(x-z_2mJ=y<tj~<yH6^{EdM_SMV=X8)237`QAtUgnuu!>}OPdYbd>qQnlUE
zVq(Qe@o2?x#bkQ!q&5S75C?n|{w5b^N?@B9a&9Tb%G|$WVsYL-Yy#M&e?!~%yIXAn
zMy(XVl5*#8O#X#_J!Pm2%a+ETojedRn1TjV$gnspTf)1NAb*YT5Nwdao;ZF#v#CV%
zKh`!1n`*GnOT<l8($8z;J9VPZ@4i>HQzfEQVz{AYW;6Qe%x06I$t3KAPY=GiUF(}n
z3u&p@OxW7ZX7bPQjP)#HfV-E`Z6=ET0ue1m1god3jfAf?2L8cMr`h`8AbqQyJ;?>c
zg}M1H$}hG8$?2zs);<QH0}2WJgvJkB`t9uDubrqKMDdstV(hQaG@@GqMLj5*Oa558
zhsXvR4nh*m<~K0r_zpu?Zp+pH2XyvW){!6wfDiz*<rRoPMHdMEhRX%`KQb^uRRmiK
z|Nn$7_!E}@XDsXAF#W$_=>P;3-ZfD`?J@@qKclepFI8Gfw<w}359-Q8x{8o$@sRu*
z38rZh3@N`&OB1AjTZ9V=ze$q|a=uw05oo?iH$iH@03~zs1zF#$7c78y#Y#cew+_KZ
j0i^c}*9!!BHyeZ)u5u4t6J>vi%bskwDZ=m^W(@xiqP#?X

literal 0
HcmV?d00001

diff --git a/scripts/process.py b/scripts/process.py
new file mode 100644
index 0000000..17dcda4
--- /dev/null
+++ b/scripts/process.py
@@ -0,0 +1,395 @@
+#!/usr/bin/env python3
+"""
+永善县增减挂钩后备资源调查数据库生成 — 处理脚本
+
+执行模式:
+  extract_2019      提取2019年三调建设用地
+  extract_2025      提取2025年变更调查建设用地
+  overlay_intersect 叠加分析取交集（两年均为建设用地图斑）
+  erase_restrictions 逐层擦除限制要素
+  filter_area       面积过滤 + 影像核实清单输出
+"""
+import json, sys, os, glob, logging, traceback
+from pathlib import Path
+from datetime import datetime
+
+import geopandas as gpd
+import pandas as pd
+import numpy as np
+
+# ── 三调建设用地编码 ──────────────────────────────────
+CONSTRUCTION_CODES = {'0702': '农村宅基地', '0601': '工业用地', '0602': '采矿用地'}
+CODE_LIST = list(CONSTRUCTION_CODES.keys())
+
+# 常见地类编码字段名（按优先级）
+DLBM_FIELDS = ['DLBM', 'DLBZ', 'DLMC', 'YSDLMC', 'YSDLBZ', 'DLYBM', 'DLMC_1']
+
+
+# ══════════════════════════════════════════════════════
+#  工具函数
+# ══════════════════════════════════════════════════════
+
+def get_logger(log_path=None):
+    """获取/配置日志器"""
+    logger = logging.getLogger('ys_zg')
+    if logger.handlers:
+        return logger
+    logger.setLevel(logging.INFO)
+    fmt = logging.Formatter('%(asctime)s [%(levelname)s] %(message)s', datefmt='%Y-%m-%d %H:%M:%S')
+    sh = logging.StreamHandler(sys.stdout)
+    sh.setFormatter(fmt)
+    logger.addHandler(sh)
+    if log_path:
+        fh = logging.FileHandler(log_path, encoding='utf-8')
+        fh.setFormatter(fmt)
+        logger.addHandler(fh)
+    return logger
+
+
+def read_vector(path, layer=None):
+    """通用矢量读取"""
+    if path.endswith('.gpkg') and layer:
+        return gpd.read_file(path, layer=layer)
+    return gpd.read_file(path)
+
+
+def estimate_utm_crs(gdf):
+    """根据数据范围估算 UTM 投影（永善县一般在 EPSG:32648）"""
+    bounds = gdf.total_bounds
+    center_lon = (bounds[0] + bounds[2]) / 2
+    utm_zone = int((center_lon + 180) / 6) + 1
+    is_north = bounds[1] + bounds[3] >= 0
+    epsg = 32600 + utm_zone if is_north else 32700 + utm_zone
+    logger = get_logger()
+    logger.info(f"  估算投影: lon={center_lon:.2f}, zone={utm_zone}, EPSG:{epsg}")
+    return f"EPSG:{epsg}"
+
+
+def ensure_metric_crs(gdf, logger):
+    """确保数据在投影坐标系下以便计算面积"""
+    if gdf.crs and gdf.crs.is_geographic:
+        utm_crs = estimate_utm_crs(gdf)
+        logger.info(f"  地理坐标系 -> 投影到 {utm_crs}")
+        return gdf.to_crs(utm_crs)
+    return gdf
+
+
+def detect_and_filter_construction(gdf, logger):
+    """
+    自动检测地类编码字段并过滤出建设用地（0702/0601/0602）。
+    返回 (filtered_gdf, used_field, code_counts_dict)
+    """
+    # 1) 按已知字段名搜索
+    for field in DLBM_FIELDS:
+        if field not in gdf.columns:
+            continue
+        logger.info(f"  尝试字段: {field}")
+        try:
+            vals = gdf[field].astype(str).str.strip()
+            mask = vals.str[:4].isin(CODE_LIST)
+            count = mask.sum()
+            if count > 0:
+                result = gdf[mask].copy()
+                result['_code_field'] = field
+                result['_land_code'] = vals[mask].str[:4].values
+                counts = result['_land_code'].value_counts().to_dict()
+                return result, field, counts
+        except Exception:
+            continue
+
+    # 2) 遍历所有字符串字段模糊匹配
+    for field in gdf.select_dtypes(include='object').columns:
+        if field in DLBM_FIELDS:
+            continue  # 已查过
+        try:
+            vals = gdf[field].astype(str).str.strip()
+            mask = vals.str[:4].isin(CODE_LIST)
+            count = mask.sum()
+            if count > 0:
+                result = gdf[mask].copy()
+                result['_code_field'] = field
+                result['_land_code'] = vals[mask].str[:4].values
+                counts = result['_land_code'].value_counts().to_dict()
+                logger.info(f"  模糊匹配字段: {field}, 命中 {count} 条")
+                return result, field, counts
+        except Exception:
+            continue
+
+    return None, None, {}
+
+
+def calc_area(gdf, logger):
+    """计算面积（平方米）"""
+    if 'area_m2' in gdf.columns and gdf['area_m2'].notna().any():
+        return gdf  # 已有面积字段
+    gdf = ensure_metric_crs(gdf, logger)
+    gdf['area_m2'] = gdf.geometry.area.round(2)
+    return gdf
+
+
+# ══════════════════════════════════════════════════════
+#  步骤函数
+# ══════════════════════════════════════════════════════
+
+def step_extract(db_path, layer, output_path, log_path, label):
+    """提取建设用地（step1 / step2 共用）"""
+    logger = get_logger(log_path)
+    logger.info(f"{'='*60}")
+    logger.info(f"提取{label}: {db_path}")
+    logger.info(f"{'='*60}")
+
+    gdf = read_vector(db_path, layer=layer)
+    logger.info(f"  数据读取完成: {len(gdf)} 要素")
+
+    result, field, counts = detect_and_filter_construction(gdf, logger)
+    if result is None or len(result) == 0:
+        logger.error(f"❌ 未能从数据中识别建设用地编码 {CODE_LIST}，请检查字段名")
+        # 输出空 GeoPackage 避免下游步骤报文件缺失
+        empty = gpd.GeoDataFrame(columns=['geometry'], geometry='geometry', crs=gdf.crs)
+        empty.to_file(output_path, driver='GPKG')
+        sys.exit(1)
+
+    # 记录分类统计
+    for code, cnt in counts.items():
+        name = CONSTRUCTION_CODES.get(code, '未知')
+        logger.info(f"  ✅ {code} {name}: {cnt} 条")
+
+    result = calc_area(result, logger)
+    total_area = result['area_m2'].sum()
+    logger.info(f"  合计: {len(result)} 条, 总面积: {total_area:,.2f} m²")
+
+    result.to_file(output_path, driver='GPKG')
+    logger.info(f"  ✅ 结果输出 -> {output_path}")
+    return output_path
+
+
+def step_overlay(input1, input2, output_path, log_path):
+    """叠加分析取交集"""
+    logger = get_logger(log_path)
+    logger.info(f"{'='*60}")
+    logger.info("叠加分析：取 2019 ∩ 2025 交集")
+    logger.info(f"{'='*60}")
+
+    gdf1 = gpd.read_file(input1)
+    gdf2 = gpd.read_file(input2)
+    logger.info(f"  输入1: {len(gdf1)} 要素, CRS={gdf1.crs}")
+    logger.info(f"  输入2: {len(gdf2)} 要素, CRS={gdf2.crs}")
+
+    # 统一 CRS
+    if gdf1.crs != gdf2.crs:
+        gdf2 = gdf2.to_crs(gdf1.crs)
+
+    # 交集
+    result = gpd.overlay(gdf1, gdf2, how='intersection', keep_geom_type=True)
+    result = result.reset_index(drop=True)
+    result = calc_area(result, logger)
+
+    logger.info(f"  交集结果: {len(result)} 条, 总面积: {result['area_m2'].sum():,.2f} m²")
+    result.to_file(output_path, driver='GPKG')
+    logger.info(f"  ✅ 结果输出 -> {output_path}")
+    return output_path
+
+
+def step_erase_restrictions(input_path, restriction_dir, output_path, log_path):
+    """逐层擦除限制要素"""
+    logger = get_logger(log_path)
+    logger.info(f"{'='*60}")
+    logger.info("逐层擦除限制要素")
+    logger.info(f"{'='*60}")
+
+    if not os.path.isdir(restriction_dir):
+        logger.error(f"限制要素目录不存在: {restriction_dir}")
+        sys.exit(1)
+
+    gdf = gpd.read_file(input_path)
+    logger.info(f"  输入要素: {len(gdf)} 条")
+    gdf = calc_area(gdf, logger)
+
+    # 扫描目录下的矢量文件
+    restriction_files = []
+    for ext in ['*.shp', '*.gpkg', '*.geojson', '*.json']:
+        restriction_files.extend(sorted(glob.glob(os.path.join(restriction_dir, ext))))
+
+    if not restriction_files:
+        logger.warning(f"  ⚠️ 未在 {restriction_dir} 中找到矢量文件，跳过擦除")
+        gdf.to_file(output_path, driver='GPKG')
+        return output_path
+
+    logger.info(f"  发现 {len(restriction_files)} 个限制要素文件:")
+    for rf in restriction_files:
+        logger.info(f"    - {os.path.basename(rf)}")
+
+    result = gdf.copy()
+    result = ensure_metric_crs(result, logger)
+
+    for rf in restriction_files:
+        basename = os.path.basename(rf)
+        try:
+            restrict = gpd.read_file(rf)
+            if restrict.crs != result.crs:
+                restrict = restrict.to_crs(result.crs)
+
+            before = len(result)
+            before_area = result['area_m2'].sum()
+            result = gpd.overlay(result, restrict, how='difference', keep_geom_type=True)
+            result = result.reset_index(drop=True)
+            result = calc_area(result, logger)
+            removed = before - len(result)
+            removed_area = before_area - result['area_m2'].sum()
+
+            logger.info(f"    ✅ 擦除 [{basename}]: 移除 {removed} 条({removed_area:,.0f}m²), "
+                        f"剩余 {len(result)} 条({result['area_m2'].sum():,.0f}m²)")
+        except Exception as e:
+            logger.warning(f"    ⚠️ 擦除 [{basename}] 出错: {e}")
+            continue
+
+    result.to_file(output_path, driver='GPKG')
+    logger.info(f"  ✅ 擦除完成: {len(result)} 条, 总面积: {result['area_m2'].sum():,.2f} m²")
+    return output_path
+
+
+def step_filter(input_path, output_path, log_path, min_area=70.0, orthophoto_dir=None):
+    """面积过滤 + 生成影像核实清单"""
+    logger = get_logger(log_path)
+    logger.info(f"{'='*60}")
+    logger.info(f"面积过滤: >= {min_area} 平方米")
+    logger.info(f"{'='*60}")
+
+    gdf = gpd.read_file(input_path)
+    before = len(gdf)
+    gdf = calc_area(gdf, logger)
+
+    result = gdf[gdf['area_m2'] >= min_area].copy().reset_index(drop=True)
+    removed = before - len(result)
+
+    logger.info(f"  过滤前: {before} 条 | 过滤后: {len(result)} 条 | "
+                f"移除 < {min_area}m²: {removed} 条")
+
+    # 添加元数据字段
+    result['fid'] = range(1, len(result) + 1)
+    result['bdlx'] = result.get('_land_code', '')
+    result['mj'] = result['area_m2'].round(2)
+    now_str = datetime.now().strftime('%Y-%m-%d %H:%M:%S')
+    result['scrq'] = now_str
+
+    # 输出最终成果 GeoPackage
+    result.to_file(output_path, driver='GPKG', layer='final_houbei_ziyuan')
+    logger.info(f"  ✅ 最终结果输出 -> {output_path}")
+
+    # ── 输出影像核实 Excel 清单 ──
+    excel_path = output_path.replace('.gpkg', '_影像核实清单.xlsx')
+    if excel_path == output_path:
+        excel_path = os.path.join(os.path.dirname(output_path), '影像核实清单.xlsx')
+
+    export_fields = ['fid', 'bdlx', 'mj']
+    # 保留有意义的属性字段
+    for f in result.columns:
+        if f not in export_fields and f not in ['geometry', '_code_field', '_land_code', 'area_m2', 'fid']:
+            if result[f].dtype in ['object', 'int64', 'float64']:
+                export_fields.append(f)
+
+    export_df = result[export_fields].copy() if len(result) > 0 else pd.DataFrame()
+    export_df['影像核实状态'] = ''
+    export_df['影像核实备注'] = ''
+    export_df['核实人'] = ''
+    export_df['核实日期'] = ''
+
+    # 使用 openpyxl 写 Excel
+    try:
+        export_df.to_excel(excel_path, index=False, engine='openpyxl')
+        logger.info(f"  ✅ 影像核实清单输出 -> {excel_path}")
+    except Exception as e:
+        logger.warning(f"  ⚠️ Excel 输出失败: {e}")
+
+    # 统计输出
+    area_by_type = result.groupby('bdlx')['mj'].agg(['count', 'sum']).reset_index()
+    logger.info(f"  按地类统计:")
+    for _, row in area_by_type.iterrows():
+        name = CONSTRUCTION_CODES.get(str(row['bdlx'])[:4], '其他')
+        logger.info(f"    {name}: {int(row['count'])} 条, {row['sum']:,.0f} m²")
+
+    # 提示影像核实为人工步骤
+    logger.info(f"{'!'*60}")
+    logger.info("⚠️  重要提示：影像核实为人工核查步骤")
+    logger.info(f"    请使用 {excel_path} 清单")
+    logger.info(f"    结合永善县 2023年耕地流出正射影像 和 2026年房体一体正射影像")
+    logger.info(f"    逐一核实图斑：对建筑物完整、年份新的图斑标记删除")
+    logger.info(f"    最终以核实后的调查数据库为准")
+    logger.info(f"{'!'*60}")
+
+    return output_path
+
+
+# ══════════════════════════════════════════════════════
+#  Main
+# ══════════════════════════════════════════════════════
+
+def main():
+    # 参数通过 PARAMS_FILE 环境变量或 sys.argv[1] 传入
+    params = {}
+    if 'PARAMS_FILE' in os.environ:
+        with open(os.environ['PARAMS_FILE'], 'r', encoding='utf-8') as f:
+            params = json.load(f)
+    elif len(sys.argv) > 1:
+        params = json.loads(sys.argv[1])
+    else:
+        # fallback: parse from sys.argv as key=value pairs
+        for arg in sys.argv[1:]:
+            if '=' in arg:
+                k, v = arg.split('=', 1)
+                params[k] = v
+
+    mode = params.get('mode', '')
+    output = os.path.abspath(params.get('output', '/tmp/output/result.gpkg'))
+    log_path = os.path.abspath(params.get('log', '/tmp/output/process.log'))
+
+    # 确保输出目录存在
+    os.makedirs(os.path.dirname(output), exist_ok=True)
+
+    logger = get_logger(log_path)
+
+    try:
+        if mode == 'extract':
+            db_path = os.path.abspath(params.get('db_path', ''))
+            if not db_path:
+                raise ValueError("缺少参数: db_path")
+            layer = params.get('layer', None) or None
+            label = params.get('label', '建设用地')
+            step_extract(db_path, layer, output, log_path, label)
+
+        elif mode == 'overlay':
+            input1 = os.path.abspath(params.get('input1', ''))
+            input2 = os.path.abspath(params.get('input2', ''))
+            step_overlay(input1, input2, output, log_path)
+
+        elif mode == 'erase':
+            input_path = os.path.abspath(params.get('input_path', ''))
+            restrict_dir = os.path.abspath(params.get('restriction_dir', ''))
+            step_erase_restrictions(input_path, restrict_dir, output, log_path)
+
+        elif mode == 'filter':
+            input_path = os.path.abspath(params.get('input_path', ''))
+            min_area = float(params.get('min_area', 70))
+            odir = params.get('orthophoto_dir', None)
+            odir = os.path.abspath(odir) if odir else None
+            step_filter(input_path, output, log_path, min_area, odir)
+
+        else:
+            raise ValueError(f"未知模式: {mode}，支持: extract / overlay / erase / filter")
+
+    except Exception as e:
+        logger.error(f"❌ 执行失败: {e}")
+        logger.error(traceback.format_exc())
+        sys.exit(1)
+
+    # stdout JSON 结果（平台契约）
+    print(json.dumps({
+        "status": "success",
+        "mode": mode,
+        "output": output,
+        "log": log_path
+    }, ensure_ascii=False))
+
+
+if __name__ == '__main__':
+    main()
diff --git a/setup_repo.sh b/setup_repo.sh
new file mode 100644
index 0000000..bf51bcc
--- /dev/null
+++ b/setup_repo.sh
@@ -0,0 +1,37 @@
+#!/bin/bash
+# 在 SuiteHub 组织下创建仓库并推送代码
+set -e
+
+TOKEN="135e5d…4cb6"
+ORG="SuiteHub"
+REPO="yongshan-zengjian-guagou"
+GITEA="https://git.mercator.cn"
+
+echo "=== 1. 检查仓库 ==="
+RESULT=$(curl -s -H "Authorization: token $TOKEN" "$GITEA/api/v1/repos/$ORG/$REPO")
+if echo "$RESULT" | grep -q '"id"'; then
+    echo "仓库已存在"
+else
+    echo "=== 2. 创建仓库 ==="
+    curl -s -X POST "$GITEA/api/v1/orgs/$ORG/repos" \
+      -H "Authorization: token $TOKEN" \
+      -H "Content-Type: application/json" \
+      -d "{\"name\":\"$REPO\",\"description\":\"永善县增减挂钩后备资源调查数据库生成套件\",\"private\":false,\"auto_init\":true}"
+fi
+
+echo ""
+echo "=== 3. 克隆并推送代码 ==="
+TMPDIR="/tmp/$REPO"
+rm -rf "$TMPDIR"
+git clone "https://xuhengjie:$TOKEN@$GITEA/$ORG/$REPO.git" "$TMPDIR"
+cp -r /home/openclaw/.openclaw/workspace/$REPO/workflow.yaml "$TMPDIR/"
+mkdir -p "$TMPDIR/scripts"
+cp -r /home/openclaw/.openclaw/workspace/$REPO/scripts/*.py "$TMPDIR/scripts/"
+cd "$TMPDIR"
+git add -A
+git commit -m "初始化 v1.0.0：增减挂钩后备资源调查数据库生成套件"
+git push origin main
+echo ""
+echo "=== 完成 ==="
+echo "仓库地址: $GITEA/$ORG/$REPO"
+echo "推送状态: $?"
diff --git a/suite.tgz b/suite.tgz
new file mode 100644
index 0000000000000000000000000000000000000000..5a9790c1c5d0aaf93ffb54f7fb700545d97d41f1
GIT binary patch
literal 14918
zcmV-MI=RIkiwFP!000001MEC&a}&pw`K({jF-vKb2w4_jXI1M|-Af=hyMd%a?p9@8
zF4fo?$qIVinGxW4U1b7(5o4Q#M*xEj*XEG~V*&yB0j}DglB*d>pZE)VPCsUPq>*is
zkPWrhMJ#Lj^y$;5&%4j$L_t56$`(#U2Gv~ljaTy%g`dvOefZhYwLdC<vY}|_8y)*&
zv97&cvCgi}H=-T;qWe4EAkkMukLRakn5s_58-=2l@6W0yTET5M{+j;e)tnY5_Db2h
zbjF(Av~Jz7#%JuQFY905vu3BQ#qrutXY94v+RABre#ySLVxL>K*1oA;n71cp14+&3
z*E2;kQ^?0jWqrC*o&-Q^=B%~)*jl@7&3tQLd(0+P%J=P=S!?E!wZ2wceqOzB*P2^-
zu`yL2nYA8%X@7IOer2(KW!9Q_QoDVs`gj3`*djZmC~#cU4N6BxBpQhZ)RH++(BtH2
zDPwAf3VK=#m}=UHlO7pin8@V7x0wa&?1D8hX5W~FZ_;vMmq)#Ue$~{{1${71toiB6
z>P;Dmj6M)3s=Asp;sHXEeTi7K<6jd+)f~X@gqVXxV7+1LnS7c)(6zsnGP;(G6H_l~
z^bvA6PH-Zyl(E;@a_#8@>-iaAlzriWJ#~si(#2zG(oQ1AKoLLE$v$E)Y*-sJV)@wq
zU$}g1KR?lH&#bw5cCHNv@KV(18m`-qf|E#Q^p_<_7L@9>W$W2RXLa_qQvS-i|E+!D
zxji{z&rN#Z=Wo^Fd@hYw)}OrCm~hAgYF}(j0tRbp)c)#+7aLdMz$>ex*8G}1_k%Sy
zSG}<EkC9U-lq07qo7e5Li`BJT)wNA{XkR~Xt({>~P<y^wy*LFA8Lz7M*6NoZ*Y2FL
zCcdtZpR*S4)Nif=kq9t9UB5VPZH!nm=bc&A`Az%!X=$NlYw?^rWqQn-`LTZLa`oaB
zfV3Vyw5~n3FP(>>`o&HA;e2&{#{OyAD&Kae$olj2$NG10Ku_2xaI{b|izV|HNk#4X
zH`e$Xd0#7hto9#M(;8t+b&kTi{UvE1DCD&ElA-DCY^nA^BQa1Y^=TQ(a=@2W2Q{7Q
z#<mos)U4s77mf@9Gr*Wwv*+MoWxA-?*_WT#Q{&e1Rd&pk@;U3~_qFlwth0CRzmM2I
zeI?Nn+y4t_p-S@G(UQyL6RNJMG9~#^u20jqpd_W1vSyrgMJ3AYxe;q+h6zgj`mO4n
z2SAdw`pCYtT7CG43;oJWWo_ENzsLZpSI<~0W3^kHCTn5=%rlCM4py}|P%sNSQrkuz
z+4ui$t(>u*J+w~GFhxa5W**c>CK(xELV*l!_$ipZP_~wD*yT~{+(qD8ZFvP|@gYgv
z0)al&&=Q#(@HkG=8DkHAk7vP98)hJ2m|Brq;XNb+bZI9|UF|oi&{6WZ7o-l#V{B#1
zF738B2z<!X2h0N|5N%N}^lOGeAG?Mhe{zMS<}AyffkmcjNB)JuMr04cVW>R15U#2?
z<9%AT!9O&SEVMjM+Ra>1=pnMCgW?j3ND&s&UNF9=2n+;#_!ryn!#~GcP(Sgzvygka
zty#zkmsq$B)vdW0qv&3qiv<W-fSg1oZ$j*%@m(=i%7F;*w6eBXzc%)(xyM)J&oX%+
ztz#!rK_O#1*?B7-?M?I-a=Ai&$6R7Y17hx^k-(^~N7mV}M*@2#!PL<TNYRBM@w8Q_
z`Rj!eSCv`CrE(`sXpDp$L2Ve=hQIp~l!)Xn=WmlVN^DJKvha|4BblGnRn#`e?e)b<
zd78xoV7ToYbHAFUY{eb$)NkX0NrIJ9kh4?PNbc7`Bc6y=yuQ#~&GXzxkeg(*H}HR#
z|FJVM+F#`wnB;eL?QfO;<ByyFcLH2Ttg9>f2HF2Ai4s5m&+~tA{z3tYgTMV6Ti1W~
zc1HK@ch7%6p8x)6tm~ikpWowW*PD=?>Fs@)e7lxE&aL&{KwDc|;FlGJ$j>7eh_!g8
zIy+W7a}`qOK){~7Q@c53-@OT$51M(|l`~tqedH}WeD53<w@D?5hM7%{%|Ez;Or%m?
zL{!(uE9Es$Nr$FRkf9(QP1bv>xI%KAjA)P3ByW@EB!{_CdH|VR5u6G6+$iM3#27Tf
zqyYb=vxUAedayJkTw#K2(fS~}45aizj-ZpvX8H)9`!V`FI*^1^$jsz4F(SSP#4>5E
zP*n3t)gbUMT}+DS?nCz(<e|kue3UN+0=vl17e?S8X+F5>#--}QDL($H{{;?wbmZu#
zpZxW}ryqQDByr%Q_YNK<-DD^j?TW^NaT2t~=CDw+j3phQEf^-jXlJy84n19}l&><H
z`Exp@K0aQZ-{jAS1Bf&6;Ri=Q1-#jeVJgilMviHNh7t-PTdZ=qcIOnb;o&;tgmwS2
zz5Ej^3fNPVm5pz$(TVEX9h9y2KK%3H#Gwxke)t}e-4n#`NDF-a7yCSX06#xFD*q+M
zJ`=!)4<tH*5WC^AT%7%b{{tck@ZK7IYK=bwL-KkM9Y|^^P!}_SN~q}ye9+*s`$!?L
zg;*d9`dgcR3f85)K6<YDY}vlNWNluyuG~e%<{>D|ML}Yk907zMGJL^7Bk@H#7@}~Q
z6bFqAfFfqWgo#KQlLxb;=lQsVHDUlR?YNdz1o97#9QugiOXW;y(L)72r<x`%x?53=
zezXW7gY@iHvN#_D{gBa1b}Kmu4%9S!0VzY%0cg5|yFc5rJGW<dlI;FVeD~q_?xUP$
zV?bKss0r3D_ZLRI0%j~?m`Mn`1t=2hk^y)_@=nxb(wRCSb403H(rw909?E1jF-w3#
z&Y15nBw@HaSTa+4-r=;e$t{TRCnUuQ@v@(Z1?)AMI1XN>pet-yT5)I)jrs^08&-<0
zJy-`Z5n=(G!)X!B7fxi%0VPN?n_!5*;=GvAc>Z9lBXnJGH=1Wv`rAWnGi-=c4lXSN
z`_~dBGneSs4JDmSxfI(AtKj9Cudhv=wyr&{Y^+q5FO$E1dYIT}u_T8!O@s`U@<?ss
z2W##wIr#C>_v3qGo%`MaBPVv*S17^$;KX2!h*>bzY=S+aqx~AXYS6AcKqzdWC)x|{
zlK&m+4Uu*dqrgb*7r<Nq1rM(%(&V>E$2(Dk><r=GjK4-Aj~N<(b7UP7;=R4(-EP8)
zr&`fS0~C-RjRJ(T4xEsKM9Uf)Ky-N^yy`9436EfB@~MK7Y9oa4k+CoiqxtS3cMtJM
zEHxY^h_-u3EC3HESBH?q;Wp3trrP8oakP1(q-%+shCz!s+c2lZJ>%-lWxUv6@!Z^9
zMv`@X)?OH|K3cDozpAV+Ri1CwmTv;VERd1J7LO88Q#s{-M}ScfR-|>cIFRXg<(5kr
zaK~5Z8dc9>6}r;dF^3|u=Upy(*2D@K;wS_|WTB^0M)7kXuTCUWoKZ>51l3Kz77{#k
z>xT%>)Jrtf^XD4<lxz1s4uACti})b2^WJ|>U~tloVF><;fr|WN<V%Ke^MZBl22mIt
zS~5WwFBuw=re%|10@h3e`bv4zNF*T`4S6WqMLI$RVhro)k5!0?Mc`r2&Q%}ZVxFL&
z6N(BHhe*+fdrokw3?N@H5r`gZAzRAjB_~8Rwvab7`I6T9<W{cN9$aKp#K{oFDj$@o
z4|*{kSF_wc^ibg5h-zR#NrBS{fd|p?zf4gHd1mEQ;~31s$vEB<-`5K=n8_<n<Q8h6
z8@LA(a9+eH<r>hz{u#u(B<dxPNsm#=nw&_)&<o5d4mRL$A$x)edSrAz!Csh1|7s1L
zMLEPN6c{(fml0coM<vZ@LcnAKruboElepX{p$IKg^_#v_Q|(|#!o-m}oag{+C;PRc
zNe<EvynpwX3oeqe5UHO!XHAc>ds&c&>JPUo<sX?~K!#X-bk>?$g^;OPLQ(h(EsOe@
zL`g7|3?$^Gq}ClQ^nDJGgP~^ObR@Wyh|5$Wg`h5wLV~vbQ#Aj;qS<I(WXWkTGGSt!
ze^x2qCqoqPFtM-SZ2KLAog>7!FU)@r4RbfsuV(vcQNZ7$$PTvQs#fy1K}3+q#exl%
zSvcSx2cS0s_B*cvN7X#YKy&<<ObOASn`Fx}CQ?8g?cOGY5P;H}Jh%h*k8=hiI1te>
z%u-@L-uJag{2*Jk?=RKATYT-ZRlF)m@EV%JuT6Q95SE)j!o<y&=*_5<0r<3UP0;&7
zh_LBxAhAXt!iXfMeo?Skz%Lpe3U&wI>WqrWz&pqY7@}kigMRaE_ybFrBq^4&GcOn9
zRH@0ju&jsXY1Z<TJuwTbXEhoRhK32N_KAgjEeI_A5}~EgaBjsj2DKjqoFOjCo`y(w
zw`}=2&6=*m)>3UhU!NlO+;_Fp&noNFP6uA?!I(Aiy-$H|$Z<l$FE*xX51wKqD@Bcx
z!Q90gf2#L9y+8F+o7U(#$Pw&vndyjqc8PUpR5#YF#fPpD%@xf-6bZ;g0l|B!sUFhR
zoThO5)!h>mMkCl8CgSgI_bG5yH@bO_5%MBNr$JL&k5DJ|%yBs5;QJr{$NN&=kHrY>
zWX4oFLLx%p4JH#~7{ag5gDG|*5p+_pbe@L7>|7%mQ_C5$NyUl@Y|f`J2r8g36hSx{
zh{=Pk6P=&uM$um-wP7-Z7*X8&UFRhU8m#4d6Cr!D9I@=0R-j$4>v{$^vS<p~XD49y
zT++GFQ3Arm9$9B@0!Cf@@h~Q<B=<i*dP&^zk@#_;y8exQV~+QlLB2O6A!$f`<!fHa
z>hb!&!rC3eyqTee)pO0A?<K)2uxp-$UTj>2wP-&N`FY|V?Mfkh=^a|R4z5~mk%w*E
z@zF7VfSCUoAm>MyI$VR-;aI#dIq=ER?jd{@k@9YT3^5l&%m+hEVBl&9s#htG;ILF)
zpv#aqX<$JoxN0Z%g#$d(7ZV(uQWK3E9|I$VfZ&8MMG+<gg%jOD=f(&OG#H&@8b}iw
zp#(<z-Jj|uE#$?<<sCjiEHe3|c2Y^|g`yv5Yo+xs5>eWjGHk1oEwtRDpMMzjY`#!a
z+<uR6V}H!HYTN2E?p0`R%e8i*svy!>RVf$|TJOymusfw;x1pB)y}1J#dr0qpYt7wd
z0wAN23oT79j|E0P2B9TOgi5p3J)w}BTQRi$iH=Z6d-8jGX2$8*u3kMOH)jhGb#{WJ
zVI<Yl>$qfyXwez>1mD6wZoB|{i<QgSUs?h0PT3dwSFfD=G++j{m>Og!MnN~Vq=MJ1
z5&WfyBl^6M$v2)5?L-TOLf%7d!rgR^CsaM3$)}kfU0rxlp22Vp>*Y;{NAT5&$>-L@
zinn!rs`m5|Ru!rENhi#UojVSK{M6Go9+){_onCQlk9UJ1B=b92DjAt^P8c4Rp2D-*
zo;S{l*nT!7xR3q37z!ym2wGDN#PVOK!r?q`PFsD=?T*M36aZR3qFNV>1kNTslV@@9
z5wOD5aU1P%c0PEWY25j_)A$aZO=z%J+=Apnd;oS5FiK=JJ5#}CK^S+b$xJGx=~}*D
z+sdIfVcy=cHci<&Dt5thKyeapgBWA>v}6ni6`Yazw5KJXJQ$kg8(k<<VA3-rglV`p
zPO5j-@!=3(lv-A0NG1Wic{DW)&O8(*ZB3#feo}4L<h{zXE6pK;vPGD{5lcp;vGRqz
zA6ndRcd})3ZCmpm2}BI_3s=B>Xu}QlxW+=-uZ~jK@}AWFw}_s$3Ts%7?+ImElvzYu
zX?IsN5)C)Bl{VfRd%K5Py(sZY$}jTTal9rS5|j@Ur$qS@wV16#4v(F3E(N!bDm<0%
zhKeN6EJzmh@U3I@Ep3B;WW=udCeMK&u~}5Y$bXBG*=ZjEgyCitPPK6(@NW{QLHy!g
ztKjBL_S1DtBuB*?f^G>{h;%B$iu<~nPiqP!J90hgZPF3)K=&oHCmA@?KrCR$^?C@a
zg5DXq&n0+SgFJQLeBnd_oC4Owcn?|x9x3$dB!yp<mW?F>v(c~rjV{G;In^81gFZJ>
zU7ui`=u(a9<*bRQ3j)*fYL>sVmiQtwSW@#_nHswGc4-kOm-o#vI9ZLUq0JJe%^WNs
z>~NN%!6?K<3UXY`s{NX17lE@$G~0#9$*ghG?JJUzbvM5X*u0(`2#YUZq-X>rO$Uz@
zY`r+&2lf%N`z9LGQ`v-hROsLnMM>oCo7bxsuh?^w_L-T=>Lk9K_V9*1az_$ydQr&3
zHyZc&-2;kSZ1385TjC5!rAOqTbSH*iJ)$p!8~&tf<#mQm^|mL~8==>g3_ZW8gh`pa
z*}0F-&K6Ysv$r|G^H3U5i})6SlJc^U*R|<d9(ZBmXup)C=kJmz>stZ)>lR-_UYtFi
z1FW8XVvn4n3KtA+1!ygf+mG(H0AwJQwY4qB?aNE{+`Ir0BJd2GFbVcx7Y9#fi8XeG
zsR%1^$^o~75Oti2ui^?ZR{y}#cOSwR6R{&Koyp_<F~_@2MYdJ5F4b^%5#c_!iJ4kw
zf=v0jRaksjd%Pr~B)2+NTfWIVh4B#;#=b!sM@t49%x<x7P+>|HV)i+rrqfDKkk+qJ
z{2;gq_Imwbi(8PXN#5B)SC=!p1jPA;6DZ*_Ay^#O;rUC{$H9~<FhlsB5MCvtF&TYQ
zS--T!=@7R}j!j?`1>0e-2=RvKL%LfQ;mU>TUD!=tUG&)CN_h=r)q1+bt}}e&nm1e4
zhWJsE`uI7>hZ%OP?)omUP{uxbx=f2ADnu3!o^maaMUA_Nmb1VSS+j0i9u}1MTC0z1
zBNyQ;?VnDex0GJWCDzwd2-AsgGEZz)o_$j(kHBvkrsZ1oky(J|hJd%sU$q{rTI*}o
z`{%2Z_b>pz_!!?wUwvFze`a4=L0e8C?S*l?d7<waw_Jev*2?-4Ut@L}<~O{AfUW!E
zt%~JudcsiO0RFYO(%nVKVKtL~J&2__Q4YuzW?&bwX1;_2tB+iTpqPC8r%(QL_-NwL
z2Ol0J)#(S;t<%=bw^-+zxJ0X!s-8ZMO`4Srh;5fB2KJ^MHun#USSBc;p&bzrpJ@5x
z8C~wuz=m4tL=<<l$1@YHjUM#r6-MKo@_J~83M?%~z>8DgOf7$w`&>akdfg6pcJoX_
zk+v$37?QWxg)h@*@jZ_?P6zi*-|7&vki28sU^nfc0f3C|cCDa=1{h?DyGXn_V1X=+
znecRTByucFjsv|g7BPz13?u^Gz@?7nJ{q>Qo@2Bvi0K)=(Fx!9>yb2z$a(Vzrh-z0
zsrDK8smOB}5YCcc+!U^k_DDNg-x#s_;RCQnWYSxp!I0~1X<ID6^`^yYR%mCW8<CuP
zOoJ%cXb}>JutuCQ%tYZBPcC<mn>WWba6pWa5Et=Aj-dCPctgWB+vZS=`%UFXjDy8n
z_09E0H3M?@pD1Dn={-anUhS`(XX+dm6C0{IG>vLFfBG!Q`l3SE*l<MjO}n%DAO>oj
z1<#Ohxq~6-U2O-AF3mkIZkdj1uES?u%xvH@yHDMQ%?!ww%k;K&C(~no(|2&EE!~T?
z+8Mod^V`=>M92m7BT!U)v3#`+qK2!EmfPLo9ie<dxC0)r65Xp|f+*NvCne^0o~dFG
zCU!*|CKediJciFx5fobEH>;6ed`$B{I!=?mCPmx2TIpA)3{!L-wi*q!;3i9vYGXYw
z{3TE_gr_j<{NJLS|FH*Y7w>DrU-a1FlTp{^<h?sdt2~WR^!-zBCglt6yh|BL(RYC`
ze^jIc^63|1?=q2pKl;%Tei_a>Hm#WzYw3J-?KXRD6da?DnF)3H$x^OpC_@g~Hp5iS
z(*Lq|wXscHYkFopW6yXz{)pow&Nm?>I1o~j5JE~@+JumjLZN&c+FZzxKxm;I$52wc
z4OdsU&K7MPwYwMr!rfN8?v}eZU0tdC*dOfb{y^G4GuV}VkMg6^qLsS3T4h_hx?=zA
zIcGfMaXf*vy&rp5wZu7-bIzR4ci#6o@0@uEHTA)(@F=oCP3sGpXr^EmiVUgh^1$@`
zIu&x@h?@HJWy9&w_1cTD<Dt>hr`7E*GH`Y%8zUpgSdEO}XJTXoA+?c_{8ORAwm%^~
z!H>T`@K5{xe`MtKi=)R!PeDQWy>mf)cm9V>olpG8|L{lonXSJDgGW_1fcCD9of}~8
z_ip3S{10%gdo~rXf9J-o%|G&g{!xDaDhOD;|Ki)1C)R8wh<`I?()susJ4p~%2$2X9
zA}La@l0g#3RFJ~4C1}Ah9i(y01Q{G#gH{+@UbVf(23Zo_OTWs!#sxW2t;ZhZaqI{R
zICchIA|qNwo5)_b1l<&I(sR;1v8J@wpK39w#zF5mA##612Ft|?(f%_s=z}p2;|kI7
zGa^_i3NWk^oiOx^E*MseZWz{x9vIe&WiYH0D@8B7qh2hBe>I$J6n$sb5`@IUUNmU!
zk;FynA{D3#zYQDyks#wwNs#dsT94fU63~d_0>-TNHZkmiIZMbM!VZ*XV0pj^MtX;U
z1<>eYOc?e&u@Q%z0UMyvS7SH~qHtF~u&M+MY7A2uOuhiDD)quBx#an=py%f{etstP
zya}pQrbSYd<|r8<C+N}=5)~o$5irkInm4ImO|T}<3#rMkp${goPNYQ3E&9`Z@1m5U
zw8NxLr=IIxzRjp5Vwv!k=0~VPui%(`CM7I}5`?tH(3=XA(uBy|vg&=OBlP!dl{S{x
zM+9tP)S5Y)17I*$!_-RI6CQd>UTs6swdENV9x~Y5_z4td%>fy309&5c0B|}*1Yk8X
z2*7Gw4}eA_4}jBY41iEXy#Oq5`vR=$h6Ol`ga^Q0U^f7~fjuF>09M;c0IRxJ0v6a)
zK!r{LYV1^qfm4H@+lRYVu2;=Yc5KkEE}WCjzzyl>$&QQX&W?v**W0H~o*kb!4KK#S
z$H&iwJ4$J&fW6=tlh&eE;J?tvFqkIp3*}j%UJ>dO7t%uO^kB{=y!*~i-^sarF(F6W
z-+l4Yiy!Q`^1|E;8G6-%-IKM~DE1n;cDwxIiHv=8x_^OlXSsTXtG^z)Idyd^&9zNG
zm!qxdzIcCzu7Y)BX`e#-;wLh+|2H<*0&1{P5gL<K^E)!a))<vz?X#{KSHhBEYjdnK
z%lZ}8pJ>XkbvfQMz3Uq!KcmK{T-OED1kwc&V#25?^1p%(ko*x+s+_dw5^tgpB2YeW
zQJTc>BOX^;s~gI+C0qtCN@*w_GT4<0C6G}+h_WD2RH-O{W&Cv!aQc~@00n3y7*ac0
z843%6)CwcO8@{uWI$_F2RN@J|c3;}~ccYi0@k1F~-G^%yTxGM-nP_5v%9n9<$ZW>~
z@0@*V=BY$^ny;7XdMSYB(Ku7AHjn?V&130Mu?7`MKN>Z}6)2uUvrA3!SgYyt6oL8<
z=g=ren;(%uX;5#GB9pJfs@FlDIi_kCi4koPeD_&eWYw!EGs$RmYnUNwjabbzw#dQO
z?bm3Lk65k|q9ejwBW_p=x<?dn3I$io86Jd5zaI(?2`$uDIQSI`tSCT@euZY`^x=14
z@J}!Zm*WL3Wv7_3Z0H=DvgJuB&;_o4iosEU1`Z>k2e3d`fuCPMqxW#oq*}(sF<{MS
z+Y4M5fYk^t0IM-<02&clNUM?j*oeHU(ZtNXls!Uv4kq`avNu#&6Wm#$P7&&$iYTi}
zc;*6$vE<&Xr<JnS)JDa#etIXAzxd9iS1E5!k%~8vcD1GUDy|LFLr|B*kI!68^dwhZ
z-J-Zqwkt3_^qYp(*rs^zJ-#M&VE(Yu{?zoqf>8P1G1vgqJHwZT<+7IK)ck6>rDuNp
z_Nm*!I~VRw-W`+O2Qu`*oKP0;j<>{~mTA8#*rp2N?;j!9(x7>|8K97!C<tTW;g3@0
z074`q6cRV1$(`Uz%M}bLgd#|lP~k54$M@h42{+&ceI%&|G(Z-(us|qc62a>%J&9z%
zsvdX_Sfs5e4G9Xc<;5XnQ6)x|7Ej4paA^pp4<K>bp^A&m{_cgJUWj*H**v#7>s@=#
zyLO?hDcO}SYmM2UlAq;gc)4QJ{P_G4x%`<7yDbN`*VN2ZqBZSmlG&yl@0{LQy0E&&
z8lo`@)H4Q%mw$^Fgo>Ch;CWIaFPFjPVoLK$my1L80P8wae5j*B?G7`XIKh{u>3R&|
zIFox2W~?SA1d}#(e5eZ_LS>#uw82Cekv;?0L>hskq48);P@F-xQW2}jptX-!^5vKB
zV|i%~yxY1pALF|#mv3#t|A?O0BB(F6G0P}L!pGwBf70jysx#?}?G~%Afyt^vsrV07
zDhiTvo`90z_<y#-7yj&4*!)T4(32?Ik+l~28;n7L7E8o$Y1B*VG%;q`N313$tF;=M
z35h0&BSf6K%v~Vfv>YKWkf!<MWv7-OaD-=J7WUHr-yt6gqkZ5~wx8W!A=8F&2jC57
zDuKOlq7qoulaXjMl&5vUyDoETk1(=JI<d#MsCsnBE8qgWf#e0)^2b2)mO^D|f!m--
z2{4$V1wxY0>9c^-zkY|A;nx6-4o(8AN?`&VMphDF%af7-PJaOfgd$p+0BbNsj5<;p
zF{A(*QKSF{*95HU$pv6lD-N-Qj&AY>tNziH0g?y{?SlV8?G&<J4;(~U-PeSR@x_=N
z?@m~5@>h9<U!8N5>$a;U8BUJN<r^|=$0EgZzMQur9-iyVdYct*^T+$r-ayvdu6WzC
z-Y&)4mG*9m4dlF)S9<4q6VD~~#<%7R9!YyUU`~0}m7%$z_jks&=UnCSH~!*c&R?6@
zf7O4p{%U>NABbCXYun@2cvNw($+fOa4SYEDcW;0A_WX(YH{{kW+171J>o%AjQQXZr
ze?!ut_}jDoKE>ac_HWJl2NeH6+CK;@sjN$clVft_`g`v6i#*}>FM8l07n_MPU)Iy4
zc$!kZvZpES>CSq#DxR%{DfnPm;Tz!a7b(Y|$xz)DSNJ(zuJ4{dFu&vWmOBUT?6})4
z*Y3%<_ZClsfTqEPrlE;&xZ~w<38Tv2Jy-`bz9vY{2jI9ZF#ewh%}bLWhLB73ANnmr
zH1+F-{+gk3>MqR<d6>Iy5{2d5&?@F`BMC!ugenu*9ZZ^IQE04IM+kNMHDM{u*Oje_
zGS%`V3SrAh-OQn_uI}1FhSkipnK-n}jFzZ_Ce%x88q_Uu4ZM5=_YSb-&wNYWuwcL@
z?Iq9?j=6P!<@ni)7*7ILJ3<a<R30$+<A(rM&<aI+WlsN^MMD4o8lvA01ql_7y5s#<
zhUbPATSKxx#eKLt+uW-(_okcsvdvE`%}-~V5dw|b7wk2O{$xekzIM7l=c~GMZ0=ar
z*P{4ZvcC0-Z+&WS+Sf7NALB0VHXu~v<fDnh3;U9#*iTvTz#+SBwwNM$&EQ>JrVJIK
zp>4E82#8qe$CfUt(YS?)kOtqRb(y;vgV5-^Y1q2;hQXUS$rZ{(7_H7@7G&k4Y?O=G
zE;$YL*mTc;4Zjl*VcRA+N@OQ>FDY!1D{9fD<(T<OMf*4j5rlQ3tTgWy|7m`lQG3L0
z>anCm$1On{1#C-4I1^PSZ8;4qz)jSYme3wTh(vqftiXDseAF>ntJm6K(g$92>1!h_
zvo`Js59h@b5wJ8kU|4p;R_M*;UYeL(rPpXHtPI3nHqna}%If9&8+9(zZn@c7!q$Be
zr&w`~7AuQTn{6d?yi7|~C7WpEzww9H!#}<YfzrdRPyK*j*s3@8km_rcQ0p3;GT|}`
zW6z5EP6orYFggP89VvvuXwpo{gQU|i2X4@jy2RHK+H(UMeP;Cu39Ra8SJGzGFzz@$
zDTFW8?=Bvajna$gzA;q9j?*=Ens3owqvRvbhcx<HQw4^v{mVaoTcOvr^vkS38Fo%#
zh=^U4!1C%z6<YxYp9+!=qAWBcz+vPoQb%FSe@Q0Q;W44pL4d)mK%lJjBL`p^hn4~9
zdFc?E6c@1N&nRLiD&P%=(WNHb2Yf;in+?@~oIt@VC>T}yh**yhWQ8EwG75nujb1wd
zdx4iH<RS{tvIt>1jHe}~7)<&({1^HwQgyOdCfKyFN_TC=d*W*|bglMROTw41C8T6)
z(l5K$W$5+A#qQgE@$uw=jC~z6N}6wM?JaL=-=|e6`+TeH?#s|m-mmSu-F3&BuHB)~
z{^`Nk_IOXuCQQG5-J3XgttLg@^9K-kb<8;A%E4b>ymLn0J0|<aGwiARZg1AzsJI(b
zZK>Kjc6raf?4Co)o<kY;VF><MVOE$C;%8DVF(JclTwwWGc7~06F7u04qM{a=$oF6Q
znxL4?GhB?0jV`!WWnF6&*P5i58cMr1#ppa{j*lhwCkB(Nk^{+}R9UJ&wP~J~%Qt<&
z_AD~+?guR4tc;&eoKFiaaH2ew2@V3%tfxisw4^<2r*~?T(sq4cvf|p&#4b6|cf0S7
zDEAFRWar+Kq4$0Xi|oF%JN{(y{PgY&y`eB$shX4^yL&S9<{V!k(-o@aHv0C_lZ!34
zV|mMscI1z@+~#Gd$Cg;eT2EL@%Nm@J%wtLx-i2D1Y1d+~;fz72(@-;;0Kp~5!@N$$
zhzpA9w635|({|18JOh=TR0H!&;Hh~%h50bdoU^oCgh74|>|@Hy39_AL9vN7M(`&U&
z>OL4~cWHy&HJ>TFVWx~0m4|lPf$g-34b|+lC_tO*Tw-%gyfV<{x^B4*HrKo_9#flJ
zh7FSGNe1aG%n7$H$$kX)Rv3Xj0~b~pEn5+XRZ)q~zDD3Ml1J_TgaVdR=Pm$eG35oc
zel+^%JnuFsCV7&O`wyEuDU5>iD0m$OZ=e7jL<%?z#tT^}3B6;J-X!uiRzk0hfGvN(
z1HG2hKk_RWA?Ynt6AP^BVe)+3DP2N&$OfsN3>wj6;jk0iAQOj$!+17b`fHf<x2WLG
zm6#LRy(L5UEiobacYzgYzHxeBq3AEI<htvl$>wXP6Gv0Ry!iR4&rZz;W%srW{p@24
zKtsSf-hK1ZoAGY6>%V8K$+^6<Z_T`w=)T!|wKwf*j?oXSgwr#7bmnNHJTaN7ze{I?
zU2p(`H!CzKLPK)DEHtErw&}qIp+dI;Uuych>Eo{C`H!AVw#yyc?o{06lxJVcK6_Ak
z_F($i!}5`r<!!IXubhzGqZwMv{Qw()v;ejN^VUCz4S0NsVjH0M%}leKzEaY)eJB81
zFxy#C5wiFhpVT#{nHCg{J1<@LneDEOHffpB(UZC^E3!09(f+QoP{yD?nQ^GGMJwD%
z?NZ+nn0c5d^*z|6{=Q(M)}-8{XVz?2MndhJ?~B&gr7_%hy+}hl6>5XyUsYO%XhU-^
zvc}Ou8)dQ&gZi4gWjE++sE3!CQ{FJAMoa4`Y6Z<(7o0bCnOcO$%!tX!!3bMu%_vWo
zH_>QRt)+M9yGw2|px43~)|f;6jFT@np`QMb9!2fk<D|7vmXGigZ6+l+v(!A^;3wKw
zdPkuh!%7wKf+NC-PHg8*PK)a^CKyIza+nv+A}_i`_brc(RXEt2uF{gQOqt#`XL#0%
z*TG@d*C{*p3<xXfF)r&`M%ug)JCxA!h!cAbT&m~5tVNu^g4S>iRX6y86voNRYiP70
zD|x!ids(|okLp&`qk7n*h9YD};dfQtTFGXJ38L<0G=(3}P<6yzjEs*EVolWjwi}M_
z1NG_xslUf+3?m=aaz{i}sD`xSBN-1>=Yfi%ipdGF;~+p!F_<B6%Bu1MPH{RrAWA}U
z{_bKjzA$=TS46;9KxX*mlunlyu2gF<<Hnn^BJ9-Oh(R8(7r4>(OYfoD{<1KrVyZCg
z11YeCh4CRg#Oa0n9#4#6c&ajk<=MZ&bUr2Q)v7Bd=2C0aIk|S{PydvXXQ`-E2c&_T
z0s;lpxsZfTdI1KrSxU&4jr}50geV252!H~sO1cQFD(fP|m#>~X0cgZ918jL}3=B{K
zt5JtZhtb{0)n353Dv*t%U=9VxP}e-bUihF5xR!iuxCGx~=L+zY(EkAo>>;3}G1c+%
zA?}J+3RqQ+C%_xn;D9Yp&jcJszBWLf40(tWM~sG@E0_7uQTheU{b%?u<g=^@Mz>|?
zbxUEiFuQqXGX~EM8NT39>7E~vySK}}!3?{@;32th7iQaM+7tG)eYL@NgV1c~quyj&
zYCv{x%+Q?+c5l2N!f|`!bUy^s{u-I~<y>37CP?Nm8MEYAf0k`j*v5}5ZqdKwKH+}V
z^m*%Nt)H*|Z2j$b<URY+o%_>k52RNe{DM8SXhHY<mz8kTq24)XoU%WV>W~9{a^Hwt
z`I8KLJO?|7{4v8RYmCpZ9XVIE%vR^Ic06#sUvAu(uI-d-y5z2%vU@l~zmRh_pw4VU
zr#17=oa-;$JaqL?wys^NYfqg>*KL&HE#+Nbu-&K&n-($NKRk2zy&%H$^@;OIi0L;#
zOz&*?hVU}J7`5Q_smFs>E8f+~=hEKR7`sS$INw5L^_A#cH2G4xa&2rV=c=5I%tR9B
z(E-zYu0XD4_04ClK9lPCW#1=#>6)H6d*4@;sCfTqY%o{dnDi^v>$24yN_9uNdSh&O
z!B=(V=-km{MXG8(BJV$(@g0HXd^K5LtKw@-?U#M6X<uh-aKWc`(YFC+SJYk)C3~*D
zrBt-VcI0YWQWf)@66nhYo>l@+rvuxRnr*QczI2q$4$lnB)orQXRNMT{+ZS%1z7xKC
zOb)&xzc8{$kUt@t7Mj*1$39w@tV^AoAH2Oq>E530eqQN*KHWWhcR=piEAKr3XUl(x
z{BMF34wH)waAFqQ2!=P`!nAxGO7(p7R`NBuYv9hoyU#1z4`#OqmF>aw_G9wP$K`<&
z@`*9oJ)WUYsau@iJrM9KON#U5V0Ztwo$HTj<GSN__Srt0^KTo2!Tc~k4GDw<LLej|
zNg;#~TADT_tP3T%godQh9MCk`k<v<&Oq&VZM#4<h;<VM`ZjtJB`@_?wiE5kLNnLd2
zNO7)IirW1%X@Af(Tczfoz4y+(GjN2gO{M+eZ=&P*-Me@1-MxEvpYQkcKAjVtoM%%c
zD`JRN%y!N0<#O8Qsr{c*j{8d*6^v%Wx7tSqdv|qIs6MICA9QLy+2T5A(fr+_0qJFd
zE(k3URPyUqKqXZQ4ER2QE=c86ulQ7cbZ#@Pl7nhYH<BUerUXs2;Fp$0lo?hHX{+e<
zui$c$BC_Vy`w*!<c_xz8rL=%t_KBrmqo+Dn2jPR#aWF`Ru9$(<Bl3(s<(f#Jjx?oo
zEKNbeZL+T{NgGilmsBDX(Iz=K1xi}PlH!$7a#Yb)DMC!rHdy+>AvNjxDFoAKyZm|2
z8bE!HL^~CzVJuEN>C79hWWetMVWt&un8lUgADvCRrO^NY1B;{Z<&Kd8sW;x$M=61n
zb&gPgO){vbbL2MkgAdfVaR494f%4VI0Jv9;!iPqddJ!MSFO#g1_MnuJ1!;7@6u=6d
z_aX+$<g=6_z%i!4ycz;i98LL%l4b{Iq1Jm!2JIu#C$ti*k;9R+nKjV`@P4CIp_KAU
zsUn!tW2{LD_5gTPOqbj!U3WZ`!NBUW@jwnVD3r5E{k5T9ZdWK%L09%mH7U&qQvWP#
zv}Qv}J1he`)~*5E@23+1X6u-3WM^`}(lb1=JDFl_v<0?a>Wh!w=*H+OX_VR5Xje6B
z2eZtkk-aH3RixGg_SS&BN|;A)22g0rvhZ{_vDWpkeTa8M-MS^VqVRbAL%cn<ty`jD
z-4c!KmRQ~^d=(JP2T)lvP_J-gIcr`%9svaOjhQc*0R-$6Qz;N|JMzUpq_(Yy2(2p)
zK!UVlXjr1`AXxR3h@z0gK2@C+t$$f9%__nzF<X*fMfeISuZm%4*keI`>NL%Cz!WAB
zP_7QKf&|s`0}O1OKtPRBgbD?_f`q6~#W1igIA4wCViMR=ln-SpVJZ0HSxeLt!OS)G
z`P9t&@KnEnn-9RorxC;O2jB*FL5woMLZB?96aoYX<tYVSLYz`$1=v=cLzIPFl`V=u
zA{Z=5KBDE8sREC!0-01@^r_(K=#wu2G>8NlN4?VZc#z~Gatl0AG_jNE!xG?7yx*%t
z44HN;!8#B%1}_o_k_{zuMP>mrgC+3p#iqjVVWV(u%6C-;95rDgdX$v_MdF5dT}HbT
zyMF-5ztN6F-H%G*Z?q~=;loS%`Mhh#x=<_&xdj4f@FIO<=4F_125t!Szzz<b7s#Z>
zu;^6KC1i7%pTlC%cfnjm`3VWH1jP>#dj+j1z$&QVDmaM-B9KX)QRH=CLg@H}^roQh
z>+cuHq!1=^36?knH;wRNb1-kiwB(5p$fU-OU|O#FlR^LqV@Vdsuh6Ly63Bu<mW33I
z%amXR{FI=VbWxcb@GvS>sZ+pv1vO|OgXkX)Ou<VbIVv(z1cR6u$hi-X`Q<}7Cu~Kb
zwoqfZhqG*(CmL4rMN+*Ei$=?o{q-XNgr_(e?!e5cqk+88wwSdtX56q~tce+Gc;l7@
z<F1%-*X&cAaTjlF2d6hVt_%mazWI}w$veI;=_J?~(gt7R${PMy6g~Ba%BTmR-F<j&
zJhwR57<xKfANEC{7r@$Z55qN$_ejj?4i<{Ejd0COFYnkm-s-P=<siV}(>;Nn5EVYo
zJ8MAR)>k^?_Uxc9W-nf_SH$cU;Tql!3)jBVv1rPiJUemr^}z*GNz7Ejl{WCE#_@e|
zXRdhQk~Rf5b9L<yIv<d(S5^2Y;paI^(>$>Q6tsC3ti>^Faj2WCYU8c@$J^sN>(r0t
zbOnn}*K~KF8`4_d^<FQ#R(8GmT6Lt7_wMAHkMhnQzu^mqZ~WjLqkZz&#Id`WT2t=F
z;~CV{!MTjWP#;%4FjqQ|Sn@Xc^2E!Ww|UOo9JhEDEX6TPaj5H-rEEzDj=E>SYHA~`
zzs-pDePB)yo);_?F-yg~Wy4>d{OFkvpZW8#@%{1A^6~wg$r~ezQj?zi@%{cYF~S`$
z@PgD}jPNXJ2~AU+v@VeD80ii+hM$f$a-@5nY?Vms+oJ9$#W`9)-jV@xxTkjpc7{rM
zyDwZEvsa7!?MQd@Bv<<wZ*N(!cg5^oy#3I4hjbZ>=ZVs|EpNeA7PFP{HXkJZ9Kt8G
ze%!xBV*b2)B3u+c^<E`kx;4T?ajvv6+Qyq(WF`<DtPPfiRL~6;uD{0eMYWu_E^;`!
zcUJTJ_Sq3`d&kF*agMHA<dcvZbO;BH1-DkgMb|Fh;;E{T8w{j<oTYZ2sFNCedY(8V
z{WvyHJS+WZo+tJ!e{k+3d?sZW<%I68Cdq@gHF7R;f^#)4klW|T?f0RFtLfWe_!;px
zZFy%lb~9W5IIh0wb#ykVZ#HQ_dYQ%xe1RipNjV9fC8x15jK))R29#E;Va!h1nepJV
zB0)T%IFR`Kti#4|nX9Qp@=jp}VKn+DM??lM%8n43iK`SJmRI^B*jJWvP6c-EgX<F4
z2+OS+tY5D@jkzK8DTJ%MXf3Nk;FV-mimnSSE%TEhP*THKK8vwx_%gU$gjSg;xiEAz
zDMz7G?u3GJf>KO+jav1zVP!mdO#8@c&!DMQwd+z;AKa^tS?AnJP8P|XMH_FJQdrZH
z14U*ujgbJFYiKian~XlPep=*vBS$q^{cwiNewJNOui|UeVw1=IOZcz~2hYNoVU4nO
z+95k(Mx^K>t7i?YmesKtzfL^0Fa|m32YnXRZw4P(pEDKuOg*yNT?(}uM<hoUOR+}S
z&t#k#_0kD?nT2}!{3B~ewn95hkF0jLLha@e$@#{bSo88;4!#5Wn1lK_`2+XSx6V<y
zv}e^(mNm-e$s>Er%37KH@6dMX`5DB)b738+l=hg=h_`?e(wHAKr#yN2XeDQX?1p(2
zACa6KDQ}8x3N80S{w1rUi{2vW;u|F^dyHcjup+Dm%=hP2&?$n$k*K=-g|y~#L$an~
zj5$1JV@aiGJzdD!R=N8W%8k@WZQZ<j?@A4=Wt^(o)#Yi;GItI(X1ip&{4?=yiZiFQ
zQ^-hB)E01tfv+_Nq6*|W`4-<H5KzsMiC{1SEK=Sty4Fy27Sv~l`j_ubk!gMk3O^;~
zl?7dr%03QiM*L4jfnwxzVP1u0kStWnVg@mW{07{-iHS=9`i4)6Fzp@G+ywIZ2>kdb
z3O#=e>w)k2<qD^x@Ff3}z$tQ3tkx?PbO}*V9|6^Spi(IFDLnB4+{D4A6<h==Rp8U-
zL4X!q4^iV>AfT^BARx~{AfVp3h?~VV3s4YT&>$6ZL4(x2Md*w|`MQLbEP$p7nQP`x
za3@mqFvtq8X!QW%e4>c8ph2D^1Q11B49Nh3UQ$>U46?*3q6Y{Xq>u{0HdIIT*#)hr
zfy;aXcd9`n1^kSFznIc@Id~86w|rOS`-g5EjKj+U!RXF`3j&i62@fY35sGM+fC{M0
ze?W2{82rO^YEf3i3^j)wZ>g>tW(+Z%chP9Q;+?JtR0Mw<vsH!<MeFB{djY!CduPur
zI<u#b295^Lz4?s4^&2(FLp%PVr!F5tk!DZOI;VFpTAkDGfcvezSKDUVLR~Wlcx(BD
z_MR4Vh-*@J2lq@h0pO<3^Y4wDZBr*+yBK%5r!NF9gzTX{-c=qTrtqnnc&-PavU6|f
z0b0}N0R%$XCL1Rjr%nZr^BILuti(LoJ<&bcJJB0Fe6?q$Cromi5AdE29uiv@&DP2D
z6XziU6a&}c60WA5&p*JMJI32UV>CJlJg%swbphS9IbaTHc}rQO_!dzwL;oeBK;wvK
zn5P=%GIAF)v!?q4{nJB%q0p)81J?#3+0i|G=`KFA*}w0$-MyGuz!h%eGaEQdLmWW6
z_CR~^e3$@Od(Kf8arj&0OcU>@_qQxkSwUB*B%B$^iq?Hf?Y?8QP92>y<|PquZ=g4H
z_<GN^o(LKB@+CWX=S~<Vf6z`)V`eB5*7zya2<z*9t#>gy=c;SQbv1t`Kb-l#`#pE0
zW7f<2TKMc%fBRyN=W6{-y&TWyw(k4b##eUmIS2h6i&RmlJwioK%|3UFdK}hwWX|Y$
zaD}ei#N9AN_VAv1-nGr&7Pn_lHwGHHys8L+mM7odu<6QX$UT_e71$N3owt`qX#Zo2
zlq+ZmwT9aw&&?9Is1|US?4IDIFdaV3ZF!dS^zpXi9CiGj1=AQV7f)79RK#sLaf>5v
z%Zi&_U+FcOS};tU?t2-UjGTXC8cpsuLpT_qU%^p&|K4#N!AiBpA<;&B;li+;b8eg`
zs}gxvjxIQTF{h7rR?d?fBq|9+&G(_(z#*<U{#|BgiTbZQ+EksJv705D&T8Gw&H5*-
z>YLklQ%_K8j?#c10-4kx_TimARX=*;jCVD6`U(YIZ!bM`qPLg%GC6yRj=NvL4g4@*
z{)Of&IANVVeVl=)Qgr_5`%ezN&^Ji;4T~$ew22=dqM2t9OgR?>)!+-rK?i0?@IDcU
zGbb+$i;KUWNX~d@g_x*5g9u1)nV9H(Tt5AY^B$ZwV0T1)<@Yc|-A41uw`3D|1CE|O
zxsw6k0T_I*4#TMmKy)1ciqPZ4Hx3N1Ucjor|Ie|k&oR%JnCmml{uyRl(qj3WqIq0b
zZ_N2DhuLmZ7S39<U@ebX%Xw=BPgRcV?`kn~CqAC>waJ9rzV;G${#}zEci+v@;^w<H
zA4na+VbblvU3Z)DEU?{LkGsAe!CP>UKB8*I@!X|W6-H40i*xGi+XVU3)+IGYxEbj?
zFKuZ{TiVi=wzQ=!ZD~te+R~P`w52U=X-iw$(w4Tgr7dk~OIzB~_Wxx61+bqyn*hiF
E03Mn#&Hw-a

literal 0
HcmV?d00001

diff --git a/workflow.yaml b/workflow.yaml
new file mode 100644
index 0000000..91105b2
--- /dev/null
+++ b/workflow.yaml
@@ -0,0 +1,105 @@
+name: 永善县增减挂钩后备资源调查数据库生成
+description: 从三调及年度变更数据中提取建设用地，逐层擦除限制要素，生成增减挂钩后备资源调查数据库
+version: 1.0.0
+author: SuiteForge
+tags: [增减挂钩, 建设用地, 叠加分析, 后备资源, 永善县]
+category: 城乡建设用地增减挂钩
+
+params:
+  db_2019_path:
+    type: string
+    required: true
+    desc: 2019年三调数据库路径，支持 .gpkg / .shp / .gdb 格式
+  db_2025_path:
+    type: string
+    required: true
+    desc: 2025年度国土变更调查数据库路径，支持 .gpkg / .shp / .gdb 格式
+  restriction_dir:
+    type: string
+    required: true
+    desc: 限制要素目录路径，目录下存放所有限制要素矢量文件（.shp / .gpkg），包括：永久基本农田、生态保护红线、村庄建设边界、国土空间规划重大项目、城镇开发边界、矿权、生态修复项目、历史遗留矿山图斑、长江经济带项目、增减挂钩已验收项目
+  output_path:
+    type: string
+    required: true
+    desc: 输出 GeoPackage 文件路径，如 /home/user/output/ys_houbei.gpkg
+  layer_2019:
+    type: string
+    required: false
+    desc: 2019年数据库图层名（GeoPackage 格式时指定），不填自动检测
+  layer_2025:
+    type: string
+    required: false
+    desc: 2025年数据库图层名（GeoPackage 格式时指定），不填自动检测
+  min_area:
+    type: number
+    required: false
+    default: 70
+    desc: 最小图斑面积（平方米），小于此面积的将被自动删除
+  orthophoto_dir:
+    type: string
+    required: false
+    desc: 正射影像目录（可选），用于生成影像核实清单时记录影像文件路径
+
+base_image: gis-base:latest
+
+steps:
+  - id: extract_2019
+    name: 提取2019年三调建设用地
+    type: python
+    script_id: process
+    params:
+      mode: extract
+      db_path: $params.db_2019_path
+      layer: $params.layer_2019
+      label: 2019年三调建设用地
+      output: /tmp/output/step1_2019.gpkg
+      log: /tmp/output/process.log
+
+  - id: extract_2025
+    name: 提取2025年度变更建设用地
+    type: python
+    script_id: process
+    params:
+      mode: extract
+      db_path: $params.db_2025_path
+      layer: $params.layer_2025
+      label: 2025年度变更建设用地
+      output: /tmp/output/step2_2025.gpkg
+      log: /tmp/output/process.log
+
+  - id: overlay_intersect
+    name: 叠加分析取交集
+    type: python
+    script_id: process
+    params:
+      mode: overlay
+      input1: /tmp/output/step1_2019.gpkg
+      input2: /tmp/output/step2_2025.gpkg
+      output: /tmp/output/step3_common.gpkg
+      log: /tmp/output/process.log
+    depends_on: [extract_2019, extract_2025]
+
+  - id: erase_restrictions
+    name: 逐层擦除限制要素
+    type: python
+    script_id: process
+    params:
+      mode: erase
+      input_path: /tmp/output/step3_common.gpkg
+      restriction_dir: $params.restriction_dir
+      output: /tmp/output/step4_erased.gpkg
+      log: /tmp/output/process.log
+    depends_on: [overlay_intersect]
+
+  - id: filter_area
+    name: 面积过滤与输出成果
+    type: python
+    script_id: process
+    params:
+      mode: filter
+      input_path: /tmp/output/step4_erased.gpkg
+      output: $params.output_path
+      log: /tmp/output/process.log
+      min_area: $params.min_area
+      orthophoto_dir: $params.orthophoto_dir
+    depends_on: [erase_restrictions]