From 7105cb6dab2c151f58f46c2898bfb837f934c8c4 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Laurent=20Gr=C3=A9goire?= Date: Wed, 24 Dec 2014 17:59:19 +0100 Subject: [PATCH] Initial add. Modify licence to EUPL. --- .gitignore | 14 +- LICENSE | 340 ----------------- LICENSE.pdf | Bin 0 -> 38338 bytes README.md | 22 ++ build.gradle | 43 +++ .../com/mecatran/insee/extract/CRSUtils.java | 35 ++ .../com/mecatran/insee/extract/Carreau.java | 32 ++ .../mecatran/insee/extract/Insee200mConv.java | 355 ++++++++++++++++++ .../java/com/mecatran/insee/extract/Main.java | 107 ++++++ .../mecatran/insee/extract/Parameters.java | 44 +++ 10 files changed, 640 insertions(+), 352 deletions(-) delete mode 100644 LICENSE create mode 100644 LICENSE.pdf create mode 100644 build.gradle create mode 100644 src/main/java/com/mecatran/insee/extract/CRSUtils.java create mode 100644 src/main/java/com/mecatran/insee/extract/Carreau.java create mode 100644 src/main/java/com/mecatran/insee/extract/Insee200mConv.java create mode 100644 src/main/java/com/mecatran/insee/extract/Main.java create mode 100644 src/main/java/com/mecatran/insee/extract/Parameters.java diff --git a/.gitignore b/.gitignore index 32858aa..81631c6 100644 --- a/.gitignore +++ b/.gitignore @@ -1,12 +1,2 @@ -*.class - -# Mobile Tools for Java (J2ME) -.mtj.tmp/ - -# Package Files # -*.jar -*.war -*.ear - -# virtual machine crash logs, see http://www.java.com/en/download/help/error_hotspot.xml -hs_err_pid* +/bin +/build diff --git a/LICENSE b/LICENSE deleted file mode 100644 index d6a9326..0000000 --- a/LICENSE +++ /dev/null @@ -1,340 +0,0 @@ -GNU GENERAL PUBLIC LICENSE - Version 2, June 1991 - - Copyright (C) 1989, 1991 Free Software Foundation, Inc., - 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA - Everyone is permitted to copy and distribute verbatim copies - of this license document, but changing it is not allowed. - - Preamble - - The licenses for most software are designed to take away your -freedom to share and change it. By contrast, the GNU General Public -License is intended to guarantee your freedom to share and change free -software--to make sure the software is free for all its users. This -General Public License applies to most of the Free Software -Foundation's software and to any other program whose authors commit to -using it. (Some other Free Software Foundation software is covered by -the GNU Lesser General Public License instead.) You can apply it to -your programs, too. - - When we speak of free software, we are referring to freedom, not -price. Our General Public Licenses are designed to make sure that you -have the freedom to distribute copies of free software (and charge for -this service if you wish), that you receive source code or can get it -if you want it, that you can change the software or use pieces of it -in new free programs; and that you know you can do these things. - - To protect your rights, we need to make restrictions that forbid -anyone to deny you these rights or to ask you to surrender the rights. -These restrictions translate to certain responsibilities for you if you -distribute copies of the software, or if you modify it. - - For example, if you distribute copies of such a program, whether -gratis or for a fee, you must give the recipients all the rights that -you have. You must make sure that they, too, receive or can get the -source code. And you must show them these terms so they know their -rights. - - We protect your rights with two steps: (1) copyright the software, and -(2) offer you this license which gives you legal permission to copy, -distribute and/or modify the software. - - Also, for each author's protection and ours, we want to make certain -that everyone understands that there is no warranty for this free -software. If the software is modified by someone else and passed on, we -want its recipients to know that what they have is not the original, so -that any problems introduced by others will not reflect on the original -authors' reputations. - - Finally, any free program is threatened constantly by software -patents. We wish to avoid the danger that redistributors of a free -program will individually obtain patent licenses, in effect making the -program proprietary. To prevent this, we have made it clear that any -patent must be licensed for everyone's free use or not licensed at all. - - The precise terms and conditions for copying, distribution and -modification follow. - - GNU GENERAL PUBLIC LICENSE - TERMS AND CONDITIONS FOR COPYING, DISTRIBUTION AND MODIFICATION - - 0. This License applies to any program or other work which contains -a notice placed by the copyright holder saying it may be distributed -under the terms of this General Public License. The "Program", below, -refers to any such program or work, and a "work based on the Program" -means either the Program or any derivative work under copyright law: -that is to say, a work containing the Program or a portion of it, -either verbatim or with modifications and/or translated into another -language. (Hereinafter, translation is included without limitation in -the term "modification".) Each licensee is addressed as "you". - -Activities other than copying, distribution and modification are not -covered by this License; they are outside its scope. The act of -running the Program is not restricted, and the output from the Program -is covered only if its contents constitute a work based on the -Program (independent of having been made by running the Program). -Whether that is true depends on what the Program does. - - 1. You may copy and distribute verbatim copies of the Program's -source code as you receive it, in any medium, provided that you -conspicuously and appropriately publish on each copy an appropriate -copyright notice and disclaimer of warranty; keep intact all the -notices that refer to this License and to the absence of any warranty; -and give any other recipients of the Program a copy of this License -along with the Program. - -You may charge a fee for the physical act of transferring a copy, and -you may at your option offer warranty protection in exchange for a fee. - - 2. You may modify your copy or copies of the Program or any portion -of it, thus forming a work based on the Program, and copy and -distribute such modifications or work under the terms of Section 1 -above, provided that you also meet all of these conditions: - - a) You must cause the modified files to carry prominent notices - stating that you changed the files and the date of any change. - - b) You must cause any work that you distribute or publish, that in - whole or in part contains or is derived from the Program or any - part thereof, to be licensed as a whole at no charge to all third - parties under the terms of this License. - - c) If the modified program normally reads commands interactively - when run, you must cause it, when started running for such - interactive use in the most ordinary way, to print or display an - announcement including an appropriate copyright notice and a - notice that there is no warranty (or else, saying that you provide - a warranty) and that users may redistribute the program under - these conditions, and telling the user how to view a copy of this - License. (Exception: if the Program itself is interactive but - does not normally print such an announcement, your work based on - the Program is not required to print an announcement.) - -These requirements apply to the modified work as a whole. If -identifiable sections of that work are not derived from the Program, -and can be reasonably considered independent and separate works in -themselves, then this License, and its terms, do not apply to those -sections when you distribute them as separate works. But when you -distribute the same sections as part of a whole which is a work based -on the Program, the distribution of the whole must be on the terms of -this License, whose permissions for other licensees extend to the -entire whole, and thus to each and every part regardless of who wrote it. - -Thus, it is not the intent of this section to claim rights or contest -your rights to work written entirely by you; rather, the intent is to -exercise the right to control the distribution of derivative or -collective works based on the Program. - -In addition, mere aggregation of another work not based on the Program -with the Program (or with a work based on the Program) on a volume of -a storage or distribution medium does not bring the other work under -the scope of this License. - - 3. You may copy and distribute the Program (or a work based on it, -under Section 2) in object code or executable form under the terms of -Sections 1 and 2 above provided that you also do one of the following: - - a) Accompany it with the complete corresponding machine-readable - source code, which must be distributed under the terms of Sections - 1 and 2 above on a medium customarily used for software interchange; or, - - b) Accompany it with a written offer, valid for at least three - years, to give any third party, for a charge no more than your - cost of physically performing source distribution, a complete - machine-readable copy of the corresponding source code, to be - distributed under the terms of Sections 1 and 2 above on a medium - customarily used for software interchange; or, - - c) Accompany it with the information you received as to the offer - to distribute corresponding source code. (This alternative is - allowed only for noncommercial distribution and only if you - received the program in object code or executable form with such - an offer, in accord with Subsection b above.) - -The source code for a work means the preferred form of the work for -making modifications to it. For an executable work, complete source -code means all the source code for all modules it contains, plus any -associated interface definition files, plus the scripts used to -control compilation and installation of the executable. However, as a -special exception, the source code distributed need not include -anything that is normally distributed (in either source or binary -form) with the major components (compiler, kernel, and so on) of the -operating system on which the executable runs, unless that component -itself accompanies the executable. - -If distribution of executable or object code is made by offering -access to copy from a designated place, then offering equivalent -access to copy the source code from the same place counts as -distribution of the source code, even though third parties are not -compelled to copy the source along with the object code. - - 4. You may not copy, modify, sublicense, or distribute the Program -except as expressly provided under this License. Any attempt -otherwise to copy, modify, sublicense or distribute the Program is -void, and will automatically terminate your rights under this License. -However, parties who have received copies, or rights, from you under -this License will not have their licenses terminated so long as such -parties remain in full compliance. - - 5. You are not required to accept this License, since you have not -signed it. However, nothing else grants you permission to modify or -distribute the Program or its derivative works. These actions are -prohibited by law if you do not accept this License. Therefore, by -modifying or distributing the Program (or any work based on the -Program), you indicate your acceptance of this License to do so, and -all its terms and conditions for copying, distributing or modifying -the Program or works based on it. - - 6. Each time you redistribute the Program (or any work based on the -Program), the recipient automatically receives a license from the -original licensor to copy, distribute or modify the Program subject to -these terms and conditions. You may not impose any further -restrictions on the recipients' exercise of the rights granted herein. -You are not responsible for enforcing compliance by third parties to -this License. - - 7. If, as a consequence of a court judgment or allegation of patent -infringement or for any other reason (not limited to patent issues), -conditions are imposed on you (whether by court order, agreement or -otherwise) that contradict the conditions of this License, they do not -excuse you from the conditions of this License. If you cannot -distribute so as to satisfy simultaneously your obligations under this -License and any other pertinent obligations, then as a consequence you -may not distribute the Program at all. For example, if a patent -license would not permit royalty-free redistribution of the Program by -all those who receive copies directly or indirectly through you, then -the only way you could satisfy both it and this License would be to -refrain entirely from distribution of the Program. - -If any portion of this section is held invalid or unenforceable under -any particular circumstance, the balance of the section is intended to -apply and the section as a whole is intended to apply in other -circumstances. - -It is not the purpose of this section to induce you to infringe any -patents or other property right claims or to contest validity of any -such claims; this section has the sole purpose of protecting the -integrity of the free software distribution system, which is -implemented by public license practices. Many people have made -generous contributions to the wide range of software distributed -through that system in reliance on consistent application of that -system; it is up to the author/donor to decide if he or she is willing -to distribute software through any other system and a licensee cannot -impose that choice. - -This section is intended to make thoroughly clear what is believed to -be a consequence of the rest of this License. - - 8. If the distribution and/or use of the Program is restricted in -certain countries either by patents or by copyrighted interfaces, the -original copyright holder who places the Program under this License -may add an explicit geographical distribution limitation excluding -those countries, so that distribution is permitted only in or among -countries not thus excluded. In such case, this License incorporates -the limitation as if written in the body of this License. - - 9. The Free Software Foundation may publish revised and/or new versions -of the General Public License from time to time. Such new versions will -be similar in spirit to the present version, but may differ in detail to -address new problems or concerns. - -Each version is given a distinguishing version number. If the Program -specifies a version number of this License which applies to it and "any -later version", you have the option of following the terms and conditions -either of that version or of any later version published by the Free -Software Foundation. If the Program does not specify a version number of -this License, you may choose any version ever published by the Free Software -Foundation. - - 10. If you wish to incorporate parts of the Program into other free -programs whose distribution conditions are different, write to the author -to ask for permission. For software which is copyrighted by the Free -Software Foundation, write to the Free Software Foundation; we sometimes -make exceptions for this. Our decision will be guided by the two goals -of preserving the free status of all derivatives of our free software and -of promoting the sharing and reuse of software generally. - - NO WARRANTY - - 11. BECAUSE THE PROGRAM IS LICENSED FREE OF CHARGE, THERE IS NO WARRANTY -FOR THE PROGRAM, TO THE EXTENT PERMITTED BY APPLICABLE LAW. EXCEPT WHEN -OTHERWISE STATED IN WRITING THE COPYRIGHT HOLDERS AND/OR OTHER PARTIES -PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED -OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF -MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE ENTIRE RISK AS -TO THE QUALITY AND PERFORMANCE OF THE PROGRAM IS WITH YOU. SHOULD THE -PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF ALL NECESSARY SERVICING, -REPAIR OR CORRECTION. - - 12. IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING -WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MAY MODIFY AND/OR -REDISTRIBUTE THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, -INCLUDING ANY GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING -OUT OF THE USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED -TO LOSS OF DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY -YOU OR THIRD PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER -PROGRAMS), EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE -POSSIBILITY OF SUCH DAMAGES. - - END OF TERMS AND CONDITIONS - - How to Apply These Terms to Your New Programs - - If you develop a new program, and you want it to be of the greatest -possible use to the public, the best way to achieve this is to make it -free software which everyone can redistribute and change under these terms. - - To do so, attach the following notices to the program. It is safest -to attach them to the start of each source file to most effectively -convey the exclusion of warranty; and each file should have at least -the "copyright" line and a pointer to where the full notice is found. - - {description} - Copyright (C) {year} {fullname} - - This program is free software; you can redistribute it and/or modify - it under the terms of the GNU General Public License as published by - the Free Software Foundation; either version 2 of the License, or - (at your option) any later version. - - This program is distributed in the hope that it will be useful, - but WITHOUT ANY WARRANTY; without even the implied warranty of - MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the - GNU General Public License for more details. - - You should have received a copy of the GNU General Public License along - with this program; if not, write to the Free Software Foundation, Inc., - 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. - -Also add information on how to contact you by electronic and paper mail. - -If the program is interactive, make it output a short notice like this -when it starts in an interactive mode: - - Gnomovision version 69, Copyright (C) year name of author - Gnomovision comes with ABSOLUTELY NO WARRANTY; for details type `show w'. - This is free software, and you are welcome to redistribute it - under certain conditions; type `show c' for details. - -The hypothetical commands `show w' and `show c' should show the appropriate -parts of the General Public License. Of course, the commands you use may -be called something other than `show w' and `show c'; they could even be -mouse-clicks or menu items--whatever suits your program. - -You should also get your employer (if you work as a programmer) or your -school, if any, to sign a "copyright disclaimer" for the program, if -necessary. Here is a sample; alter the names: - - Yoyodyne, Inc., hereby disclaims all copyright interest in the program - `Gnomovision' (which makes passes at compilers) written by James Hacker. - - {signature of Ty Coon}, 1 April 1989 - Ty Coon, President of Vice - -This General Public License does not permit incorporating your program into -proprietary programs. If your program is a subroutine library, you may -consider it more useful to permit linking proprietary applications with the -library. If this is what you want to do, use the GNU Lesser General -Public License instead of this License. - diff --git a/LICENSE.pdf b/LICENSE.pdf new file mode 100644 index 0000000000000000000000000000000000000000..fd83fc241301c303a7ebd9dedc1a128daf05efc1 GIT binary patch literal 38338 zcma&tLv$`ouqf=^j=EP3YEpx)dTKS9 zqNq3>6Fn;o+4xOD5eyq4BcX$_H4HB=gRGgog{vhY`+pH-25~D}S2Jfq260;>S2IyF z69-c>etsAiS7$RLI~dR0TfOOcJdWhubq(Do^S^i?w@838>V+*lsJXTJe{4=*F|)L-p81N{cNw?EG}C({C3yS$XwH9L8|+*lD_?#8S# zu#Es?T=LFlZ~U6gJiWXdkBH{}+zBQ4wJ>oD42u1Cn{xM;3~^+Cn9>tg*`M=(-0#mx zD|;X7FofAG{H~0?uj1_}E6mb<-AEN?!1u)t<{GErY|10qdTj+7V)i*J~Ekx1E!g zE%r{03SsFX0lRZ93x+#G;QQ$GNFlivF&%;>d#MW{ds&-c`kI-BP$Ft68!^L|oU!8b z*6Kf3tdbUa#Tl>}S3=mA`?uwrK#*ku1z`femyrzUVe|c<}i5}Y|J4}Q-9yLMgGXbT703r8_aob z)qjbxcrLmemKzY$@ILMu6hyRRUymQJb8_^-V)R;OQq@@W73ddI04&kmk+dnFR(KR@ z&T&r5;rD#HkxXW39W@vb=LGD?+xz?C`)gA@4V`fhg><}QXpqX<`X`_xV%vr*g3Qvz z?su5G8bmFAdw>R5p{`4RcnKF7<5J?^SULAS2>EJynZ0wr4pJR>sPV!ZS=a1;_a37P zCc^TjM{T|~b_W@Isg*_xZ|^HMqu^f8ld@Y3zw~b&dw`IzD>r`q5zN7z%QQ<*;A3+`XRW z^=pBdIE#B>Ee@ZQ&!qxy6VAK;o4Dp8cP_v&^BjX~&gU%;2;lTgyC%{z<$ulF9B7J~ zmB?VrzT>KJB#{G39ylWgImouXSOkw{Ru-KRdVKd!k_Ti$q+j8R`{!&Z;BaSTm~CRS zN5lW^R?*t~X7PDO$UHCGKHM!b6CXiRH&J?Hf}5tcg3XaE?LGTM;#wyUXL_vgLkyzp z0YUV3!zS67_a8gM3i|fQ!WNhQiuv}&JdJ}%yvd>Mu>Qs72!9s_NiYPORm(rIWEDHj z?VuM_j$ziPHZ-EYw_{xtdX!p)@*_{Tla-tm%|n-_=bRywjC!l-uc}bO|7YoP-U`;3C!x z85IZy-`K>)%vi5B<|#BLt9OX91b2p&8=rq^%H+)x9W^Aov=A!syirX`DLZm$J}#y7 z$}YIU*gmst@n(Js_CwNbeGtrnZPqyVdd#yFWzYD4p?#4%55xPeBi*VH9*?vrFbqX3 z8wgsOILX1 ziIf8YHjCun6DZGGeq1{_NDxrNLi(Xs1oR22BU>SG;6w}v=mhaQ@e`X&IU{=s8o!*r zl;)YS`88`wvNLY7bIo-!Wo@B{peoHKtTI=ba!LBk&a%6+<9|>o4_= zhIDmo$szxfC=hhGXvW-3NT$v0cl$%cN-CNXNji=}+gSX_affae$E;>mT;UAXJfo>7 zKs8220qL+K$FAQbXe`J$gU>bQi0ngo70`5E%hKYRh?SjrUI z)>`+Ua0oAoJyyqO4spfq zMt5_1@hnUlhf^3}+3zgvuR> z^gX$zr!u4lS;Zgdd=*5}_sz`{UZS+!h9k{F&v671q9uteBP2;8^9?5tD)(tfCFENN zrk3Olq+n~~O=QNpcUxeuHSpN!7^C-VR3ulx^Fm#BK+B>$1g6*|Y9k_Xr+G-f>lk5f zUP;!ofB+kA*v-JoypG@~H_)sjrEw@XE*fUz**yicL?ehRFYaiE6#!l&?UZXAfD0-K z{2&8V2{Yg<3{(AI+O=%Ua|qfawct-mRACm07GJ!wl!-Q5yjZ(TU_Y{YV`hXns6&;Q zhSpa1h>~tyTk~?@mIIDpYfuPQ?c||X*=~>0 zT{W%mRr?wflooNEwA#f7IE$U5DYTk_Eyzx)h4W-fR-4)& z*yDstk04K$bdp}RABV2us#%2Cgzu0$0IywiCvohHpzCiCNpf_1(%axZA7GDK2nyL# zkg>@tHDX5TbEgK)v%LDS36hJ|>=Q3>qnDiqfrw6(%xZO|D}4W!Gys)yaT;AT)hwxH z?st48EG!gha}CT)MJM5mtmL(y-h^P$b2MP+TG=ZEKAo1_AzM^h-5(2OpCC17ShA0c z=RIV-z_>iNsaxuDx5wx8;bb+Oc(xUW)9BmzwkH!&UD!?)+?xo6wj`!-1b-+9U*Ik? zt+w$BmrY#<_iA0&PWqzAyM z;eFQ}>~P)T@tXuX<|GQGC=6ccWb9(Aq>9Xw0!q{Fnt{jshp^!f>@Gu!`N@JG zyB%bDt;Jh(16o$9)6$u1Ymj`|g=LK{+$70!+P9mL$630D9?Mqqxa&}3<>mhLv}VLv zGG+p=bnJA_#|8PJvP3C#+4HjP-%<3&-@pdKG|8{A^%ABb%s-aN0aqW=2v8W90)Tsg zp?U^oGjF-vwe`_7R%ZY)@f)1f!uQE}!!hzqxtODb;_Wv^N=6J(h4q>;+GTaWBgJ+U z#@n;jrF-^WTEfhElZ4B8FH@)9b&3+vzN?IBoH&9@7d570iCLx;H+}eL?{*s>BLI@ z+i5g5_r>pb!ZRk)+Jhh$upj_zV;@4IJU9TJ4v_R2CPHUTGfd%Y99nJXyB#=)gfg#L zY(FE~WJPvAG`i{@cepou=hTwNoLeA8dVN3?oBy(5uCVVZu5H^YHsuCP*xI8PzWtw z#n)q(C#Y>J7cs7O6G}!ItPiN@QKGcZXjo&lc3)dMDw^58`oo|LVroqNddcff^m=Vz zI}h?2x8Q0!Qn09Yy-JG2CWiA|=l0{8g62u^=facA7yf;o^EAUC^VtmS zTrA@(uX1)w^$Bt>Cwj#ed7pJ((`YS?NXZbZb#$7M3Ra5>p>CwJ;EtRrNVC9GVJP2Vv5Nm9^$l$zc*1+X%lC@2rG2Z1qD>vOB6|FJd?SIO^&tyse4$JYGv7vzpJ( zNfokCnEsF55IZ;XLrT?CA4?_dy96};# zH*kkJxmdY=MpD}v{kC!)ED;F+VT#&T1oihm7ej3-%U$_pipS3D0pDW_W@}dSNca!5 znj|>psC{lyEwtpVZd2jkpw`d@D2NK&x98y1RHcU`2v%W46U}?|$Krr3>qLv2!$+|j zESJU|3A=Bn)&?LvyPwlb``psM%FJL{vz zX7UOJfzx8W%_VbX>Ze-r5P7`!r63Az-zwzAteBJ+gX0mx0g61iV{`m%n1+i$8!BJL z&2q)i6Ogx9_RztBGCFYc$jP1gcP@j#@l>=YG!q)vrv;wQE?ULo92c1Ho4iv&Q_ z9}=G1$)~(gI(tkW4Swr5K&H(snlmbmyl%xwq+Btt-J~W)Y0eVLrMoEjHEH-R)Zp|H zbp*la^C1rI62$Nj^B#tt{0<7lE@A1sRiu%<#iF3ZP_km2l@_(NP(v~y-fa0(QCsP~ z@**pW07d8%EbmZnL(@p_0nC9->6Bv_j$b^zu!3CaAK{H5GF%9a-H{gXbd{PuQ)~ zQjEc^!)DNKS9JiR!!mBC$m6fAYTNEDC~LbEb?DK=)q1x&_Q%H!5((CL+M{&g2km;I zo)=B{$9tl{2Ckq1A@L-{YLp^p%As0*lXUB4%SydE(oyhj<=$A)jMfHIgp8a32r~aF z4d0&Es7=*5nFG215RyMpGFs+y5W}?LT%#5O#JM=%ZFI_$L0ZXMhvs8bkzvO(=%&*b zZ{6`!cFS7)ZAbu91ga7i+#)KM_+w?5BC4b}`$v=7+$LKIdNXFMb3cUg=kteufqPYt!GkspdxfYHeM8s8%)d(1CGAssEL%GFG~)!@ zR+&8>!LLL1Uq6Da6g^weGzYg5SMkRe?0+M^FM0*h&dG7UwE#lVCNz(_)k!8UAo>@~-lYM&kR9ZaO<~E=LF4`*AQt@Li4MCCFkhfh$_= zz3s{tWRg5kxAVh1tk1)qRtBr6v;|O>@Yfd5e4o7v%0DuC|I`?=6wk3Xy(l;TGA;k~ zU*4bMN;Yu!>+zN4Vn@P&hyErs(F=gwhHFI~UGa(C4s*qXtzso8EwuS}6{#3t@Lbxt zY#1}uKpT-Y+Mt61mi&GF1rj1O<$CBnu@PP1hji-UgfOguyKk4hro9n*a^EjecJ13v zAG>iEz0oc5g)V}11RMKs(=_4wbsk$!6L>g;6|H2bT zL_F188XD8h^rWN&s4C_ealGF{vCSl5|_>@)85TVC~J=?d>K~y3Lcy$xkNUNiLRUrv7UpmnHD!j{gEJ zGw4=j`rxab>d#0z1diC0@C@SBZjMO#p;^J!WG zvmgB&o|N_;;%dHNy~G@VpyCN2^FwqAzcaE{>WioZnK-Ri`+ERL>xU*L6?<#BKL4J*G~R>L6Qp#tXzcU2_RTv&kT5<&SpUV zt^j2nVij9f4#q#zE7wzwTEp{IVeRUNt^Ynex#ODtq%gBL{lB{X-|c@y{C^z$e@-?| zR#uq*8)IVmUtG?_`v2nPy?V3pxb4V&cPcZpEJ5wNUzhgo@lJ4{K|sR~#26$It3MMc z>PyzWG?aIVDRW*4J@N5c{Xq^!99pJmTC04`*c3v+iwHS#=-21f0&|~o+myXl<-bD_ zF&q@PXFk9E85DO4x+Z16GX@9q7!K?K{|#bi?+sWuEFD4)b9D+^=D2y*pI)I`l&kD{ zdvNZy;r=c!XQthodA*I^F#}>PtgZH-kwg%5?=3bol=k?ze zc0U2m?p}p}p6x7g_}WJQgUgMJLcE(7@caYdUlM`H@cv#3)rpfyG*+>@7=EwZ?*sA;tSOetq`^ZM@g1+XkXVl?K#$ZB(6W;mv*?L3xc@R zQ}fwcf=;qF(VDv!*=?!4Zcp%hTt4EnDIeu=?$Nv(qZuWO!=C)BhL4l5Pfa{_NZEk74^xyu$7qbWS(5&s z5AX(%?%P;DOd0$W_<>HI&@O&0`letIB}Op!m)a|1s=G__BQ5qKk%{jz?=nC8a>!4> zB%(rGS>ZH>bF2F3Q_9^@Sl*CzccxX_p`*a}M0~g}(vCi3l=(ExLAc9v3sLg@?uFX&vOT=!o z`&Fm5MY;sMiR}4kxgR%t`*f0tkFSsqU_p#9Hl%33%q*kAr0cl#LRk|~+JM#wD(sQsc*wKE z?HiCcGt4|wCe`F0QUis~v3?sy4a3+B8jm@o=zwfK@F}lCF_|R}ahnSEqSU(N4AA7f z55WrBu}^`>pl+ggiw4Nom1Jz-7dcSxK@mrtJaDksPgZj=&wn} zD;d-wh|fck$4^hn!Z|9rficJ%{Z+rV z9}7FaskHqCNVIvHa%gqHpH5~@wLI~ML*0FL312c%-jaoe(6VnG}D6dU)1cevIWo&xEf_V&02-vH2$QN&%Ze zUwlTH2?j!T9j}Blo*NQ}I6*@;fj3~_%rRO8W%F#PPF|$;$-vJ@78$XFfDuD!ziL?E ziHC%2+2(-pB#R(Ryqm%o^oPAWWtPDgotD{gAbf?7SY0$F{SFl#=CT3Isl=soi2=2E ziHcYDLcUwG|CEc}BY6X&hON2dG6gji9L2L=O-J(aGLnvzK5@l62&CB~(n)aIIkyMF^>U|8=FZ_Oj}6t)f$LUOpz9fUqY8Uh2& z9qde4j$#H2O?@BY6a}RslT48|#65{uOI~XniG|(}sZQ?zG&Cyb+6y>fs=G3TAxiEa z2(Q*qNt`Ss(;a;)^BOJ4%7AS&N{1T(hsD3$XmK|MzsbXlkIJvq-5c&;9*UMeg!1V^ zi$G1o=j=nObX3CP+!SK$4G_n}JO8VqswRqGF)<;grtp}N6cRWHOw)lIa+gxI3(3Sh z$HW{_f{QZw0|B`1sU*#{ zm@U<~InTM(Y6hKr8qLl--5?y6HC8mo^K96(NlvIzYaU7GbWg|I|Av<^myPW5WK-Cz zicYN<6L4}xNUKtX+E>NNeuIdBUuS-9K>!|XN~-`y&7zMKSeY0}f>Hx9$%aWa{lwfi zaLg{4UP7HvuERP(Sye)g_fKs(?lBd!F~A}h5_x@yA{w95IKV@N1Y*JK4d7!J2`8RT z_M3<{`lnA6*SgH52q<|+CuUa!^G*VVr1VMK4g7-k!`olGoP$-Up|HW>S@;8dUCUa! zoBieY`(#?7%Q6yN=HN0rnX3WCsi zy~VT=K2}^n$`Lv$!rR9tRxX|-P?v*Wk%mzOMA`tee2z(2T}RNsrvETM3H&|R&EV^p zFOj4GD@pkJrOrx)_0ZS>nPk;@UI>K@23noN54o(VaH-yqM z%rWY&T%?gxnUTk&1?m)TQ^j76NVTn#^VprTe(KZG86~TerOqJVPc%GqD;gH|H6AbH zXE#hkrOw~{=h`B}`Ft;Lru75~IZJHKQAu=1l=*$5IuPu_5Nf^lB4Wb=W8JnNcO~(b zm#&q@N`YJ~T|^eu=S{U!05Nd(p-Iim`s^RpYX{`z>_?*3y8>uR@5&o2C<^Lw*NwmW z@vZPMnT=7#e{?K^T+_mm?OI|&hmRuEo%3d9*w_=*bGp z^7C@Pt*nO^^VTQPJ_3{98El6SXHL!Q#dE3taI|qj7mY;HEJqG&N}@h-cOXJmT_&tY zgmBvw8v6ss3=SdBrDwmfBrQynaiG4Cp)RaW zO;)wOqodhglx2+BapH{7QZ8IMx2lc6vmM5yGmCnO7i(`F-=yKk?>mqh!_ajG_q9d3 zx(&uYY@<4a)Xq%%9_EHebpedcuY0)J9T~4AqBTMy;1`lHSW^61myV63`#B;aOlM!K z7r-FSK>Zlra%ejCo`awVMRW36@{V)}`=)ri{Sd!C@`gQvV>FNMH5H>w%S-F}w={oa z&}l=AFiWzcA9dq5U!OJ%$2vYVUPfA>-yXHNc~qV8q2KVXD6;PA&+uskO1d%#{LE_| zm#eQ2t3POs2##s{5~^z_wDb6ezF7^P83>V4%HvYOXRi8uP|^9BNGy1g4^0n*N=g3m z2mf`hTmm0iSn&6ccc8CUsEcwLp&4VCb|zo6rx4D*%P9H7F5&*ezF<`$vUr#Z38HsbrRNY^vre!#$8`p%DI zLM9k6q38Uzs}Fd|NU9k_J0V;7f~n5aLMn$|D5aP$Tv{@!c4^N`k;aqZL@mOWWv|j$ zPSrhT2KU9FyhLf9wRX2-b;v!eyO63#5ov77uB07BGmfgfz_U8Mf-}kHoZ6NFz=8po z#yT{dNGi_F;IBkD*Xa&JOeq*cTy~IUXukGHa=J?9VTZmH+8?6*ze{wCVx(dHLqO_Z zDk$mCg$TM9n7n=J;}?u}az5K@A;@ItG^kaje$nn^8ZHo6ZWjd?d$L}b}ep8#cD z*?kXEhyaT%rtGmI%m&F6?j=T4x2{QXWVbL34V(vc*jf?Z=9%cBtotC@CrW;H&V%(D zgM!gJo3n+IC$@(1Mr9XH4?wy8)QlNj(idIifVWFZPS(b>H58a238r7;pZ$(i3~YJH z5Hr0N@)lj|^gpx2ivR)=5ohV-gnqDSqj`sUSkSO4xJd23QWTZ;a@E}s@N_hv@eL#@ zz~$kQ4H#j{6bfoWDAK;Z5KLqkv{gXmv`ANrLYIa(5_;r2U7P&?irI1-w#V$0j}$^M zg=HBKoi!!xFb-c75j`afrgLbty1BaO6QXY4doUT;);z-#lu(D_?FLaW<|KJnh_jfC z>?y%@zCxx05d@3~0?WcDSgi5OTQ5-I+oG^9qu+0slF+2QQ`WmNFw-}!UOUJUAvi=> zCd#XE%}X#MpTyfFY-qNlZ6lx%TJ*bMCpL^EhAeXS-ihIKvCtFO@CEDcq zGz<@oj47*@pFhO`Mc)pt4Us$di@8DJN|$>978XP;{rXs80AW8p!E-TTX<^;+dxVPK z`MJPWM-und4uUc(K03047U^q!Ffl}L&867viRFBa1K8%q);)Bvca$t8hH~7~7_$yA z##$Jt4XM5%v{{AEvw$LUn1oaas5SGcGiIXc3^-z>$fi9}dPzx9j9cIai1Heb=zjl6 zXQ)^sCCUgipLT{p%|m=rQZeg|F_cjz4vklC3hbPI=_ed_%FHS3j2l}Irg&~$(p0MA zOv>#t`8!uzOy~`N`SUxNAODR`7mO}Udee+$zzF(MRab71+`9;V55$!b!$Sg9&Y?X~ zMKef6kPv+~a@Flypy1+Ck3~2+RzEJBwy3rn#g#YC?ZHT{)>-72RA1U7Ay}EC!%(Aa zX1O<@)G0T)^&)<|DIup!zvc04=;*UgB+x&4TF2Me7UoRLJ?1g0FYD>SkV>|1{0j5> zGmREDlaM_~yv$qDBUVJ*e9mTUg1Jl2=h@nkmiC{&cdebS{@lY(^Px*_5tYNH`c~}E zZ}`1ryR8C8K$ZNgwCK!OMJ3+q+9u|DiRp~XU;g{8nv;i@cSr)1Ha+X~A*`3TR8Ls7 zuOC>IKVy+pik_f-JQx(C@EE3eqe3qbgPEZN~ z+vCZ26oWe(qw2d&GE?eywX*~nDLc|_-V`O()03oL0ST=rqW$SnXN{~DHx1y8?xlEY zAm|V-HIG;uyHoq>KS~+aJt4%j3)Y=m4b~V5l~q#eAWzlb6LF+g!6P7|*ZHDa>1%@V zy+tt2XSd~7!=@IKzVe|N9{RkI77Q|=2A2yUbL8isRRIp$)&Q<0U|WUBaF&&l+@)NO z^9?^I1BjBWx!+QlE_9=tq0wh2>O>s}TJYRWyhV89-W>i#HuK@A(?tHWd;2~W^MKx? z$p9cxeiQ8A0Y8oD=m5(+F|W@6Ik#axXg~c}D@+D2>lC!U7oL!ULg!av09L}NW7;`x zL(DXr_en8x7X==jj74AlH<_g0Ekmg2{mgOVoJTJ~dDxEn`nEc*i=zoCku|b(6z32f zTDpw))+$SZ==j3;byI7nvw98xQNcdnkI~mdPy2>{l zPul79`V-w%39~ULMgGVoQEv)72BLi%E1!(`ja$38yt}N!r}-Nq@3Kbl2I~87erDF1 zz;}JtG{2CX~K|k4>|U$rIKdZ*7_HbU)PtnlB3`4H-Piqm05l5 zmqXxB0RzAN&*k@I4ng>=>ZTEH_Uz|%N>_H-;>tJgkZo*bb=@Z1tt$NY+4p$Lt@-QE zzI^tU`tSv0Z$IyQ@#kOH%pHx^dxE%DISseJ$Nhc387V55KHQmWh#2OI{oaqm7lIYp zR@g|K{}nYhZs6N#`&RzC>aX4Y+HTuco_pG@?48kNl{xO1(`WZ?<#wIjZ3o{z+v+Iifv<8{5j7{rpAEXI^dWXBU~t-9_Lwr(N^6mB*I4-(foVl`uqP`)^~N`L8~fs zx7mGTeVu(z-YJxg*>+WCv5&=V$>u~ogR^|3J;V*6|2=IZ9Wf}wg@4mcFwTC-rAbd!16m)A+w8U)kd63@w^`m- zpmSJ|dCLO!NR6sP-oX4%smo|R)UlI@o5gL7ClNsh+!6CWsAnF&!Cn?+EqfDocB<$jKYU-_-2Cugt|oEFKLr~hyJ)LV{7!m7ZEIz{AGc3`EopU z+yvp5z-{<}Z!&cRKou-qTL@jc*}!b=ZF_B)%)iKoB7X#ycaw$#lPF`t+-QL> z=z~znnQYD9tzlaHau!>zIc@*KNiVIq>Yj5o8#MXvNLZ)>t@+;?)^1h4>gEBN zw#>4^Uf|E)A}}HBOq=rhN-3a;Knn0lx{r05F4AU7IWFm~L&h?>3~Ir~TVl7+1^Mqg z&mMaiX-ju!f!u>$4c;5;xY-MQFZOT-?Bkd#Euvvd?8z5YMq^mPrAe8uh+-yHL9<0m zE3aj63Iv)_ZnQ%u_O>KetkpUGW>Lp6-oZ{otNM3lno|NpI8e6J#hY5VVvIy#?J!Sb zC5f)ImM^p=fN;^fHR4nmS@Pr{m_<^%tEX46-&JV$lm>{-xWW&|mWeBQZF#%u;8>h; zwE_P%Gf^_MBF6cwWLDYz$Azf_s!YHY0c#nMdiX|ODbTdXMZiPi)f6`arja^BN#jv$ zqBzHGI*r{%5wfSNu zHP!(Xz6UYnfQz%*f`!Z`fFsAj@~ zbI*9I!rfx?pl946m=t0nHWeNl5`x62VvjP6-dA=$@)YM;UdE7eAUXtTv;$~Yccp;( zt#ve2imK*CD+oXf5;81~z1Vybphzd123*bPN-U;@W7)H?lW_yM>J_WxeWK3JLqrwZ zAusc51Z6M4V70GtQ`dCZQUrJ`y5h{U*F;TOUGOp2Iw1!d6iyFM)`zHT_nQH}f|Z20VtG+pI-?sz$J0o1n0^Iz6`F!(A&6>bmQsa)zVD#dF)h2m{2 zR#?cpkGYjFa`ju4J`Xhqy%s{AQ7x{lwA9?)3_=Shr~>VUbFO54e?Y3x!$~*bPBp5& zgp@(uz~o?^(`P;`q8?pZYZw18?i0hOqXf$mi)Gnr!1od*Ic&JotI9dF2i%wQXr5y` zWs=PYPJ3jk?Z0B6{EG_N{kN=EXIM0wdl^ek-JoTM9F`9~xTu;Ua``AC z7+je8o)?nSFe`h#L?;1EsFRKfJ1d*c`MoaNF~6}whw`Q+-WdF>V1E3!=H`9&QQLZh zCB6^NcQZk4nY2%oYSwOV>A|>2Cyo!F)eCq^i~to*$h$<0vHpi#bn?GVBKUCPJWB zjH6Bh*B^jwF38gS>4mucJcP-EZc;<+w2^`-$;6;>LxhweW5zU;r2`VmEPeVBd0U_o0V&$E=4K&ADL0Pnd9B(&h)4a=o6Jty zjARKGdTHmesYQ|_$Xm(G?8FB0$`Jf%G&NA}hS%d^QHM-<= z9i_M$N&lLhR$5*`B{>u96wt<`{>GQ|O_8chQ8znA(sE!VOVk29LC=-Vi{FH5pP^(L zX>-!_6hEmD^p26F7c*(W0TwE@d_o^bU)hZmvjJ3Gm?vC9`^|oPAeVkleYwp7NI6sj z2MtSFg!}dt;;2UtN&zXD5xCx!p&?=5|GvFB)DXb!ezB7;6>lEsous36bR9ijF^ zu|5`OtWr|@g0V#Fd5WSwkV+~JJnM7JtMhBW-ffRjFX`9R@My8S(kW%(@CkH=OVrJw zyAWU78KUYf!~`vmeh`1~`)Ay-79DlL5L?j8_*fK;tyUe9oqt7+hVPh$6Jq}R<;BUj zr(ldi{Tg5pd-jE0=%qRJ*h2Xv{mhUtXDBL5~$w7+(_;+C(ET+LFX7zMx+zDIyP!x)z^O; zQRUuKj(Z4bJu+NleoB|Hs~ce(v4%HRqDPaT24>FGY9>46iV=jGFThiv*h0{)j@I^L zd8suZ97)0>*i@BG9n21oZn&8Ao&1(dqEhxa;;?B8*jf_ale`o4D10nYV@nJsSJj^| z4^=pFKZR<_Z|j*=MK~n?^eg#zBGAZuhbF^iyGl^U23u(deV(IDDw$`QKr z&?{P#z)X@P3_Mn7Q8XbhV;BXtXnR7UMjtF7bf(4F0o2e#x+UWTK^tWqDJ-hVd8T)x z=_YrS!LTyMAIK>Nd;0E(*iYH^@3qJ}IjIjs%M41a;;UQMR?HQf4qw)QUJVgu6b88< z`bc?3xIia6y9zOBr>Knrn~%KyuM8I9=zwereYnt6+Nw%m;S3BYKPM7lqPcbYk{;Hb zB@vR5zqcGH`jUzJ%tPoQq%19yxFgvmBzd*uD(bVjNQX&{>~bpf!f+W7uukO(8A;}1 z_9tSvHpxw0+(l?eym6~wk+8@d@gKv1q%T?1sz{uxnkXyj6VSzrIUh(q2KjZ z*RaajqDKc0;jK~Dva2m45r&EIqg9WkrjyQepDb$!%ngGRYb7WkU@*~%NG&U=keDp% z+BZ;i?UgY`Ng8s|ym|_QsgsYx85u?2pgrxstrOW$;eF~#kckMho;}V=5++(qg=ZyA z-owgMlcg~dr~&l5uYuHJ_dp$u6FiICe7)TldnM3k4Z(xV1+-}fwxD2M=+q0SDZ#lJZHag|tJ zdX9obt$VB2>@HzuGN@b;lZ zA)>q_{EbKzy<|d$GY)0wd+_YU;3ASneI_0hWZ~yc-DR`#G6NdH0W`Fu>pg1x-H&+e zOt7e_=vFgBH?ERdL8Y;7O%NMw{1p$LI&F(~|AHR6@0F}?P{)i;5%J)eY_kAGJPGGg zmh^|B+2I*+5Jq6;vwS}3`DfN7K}X{t=4x5FiddtvMT0UB^`XnrSwg3piasn~)$zQy zfSX*<70kLN`Pq}3Zz8oLFXMG%G22MNFC$2aweFMDF4f^=H~vYw82U}TxZR)ixoqp3 zA&)pLnq~JUHi=Z&t!-_ZLRs9oM{khuM?<)OF7(SJm<)()j1rM)E+Y?co0nTMqlMmV zy16!kKy=7|FvbB=uP!fqMW`?{aQfiOysKsSUxg~@oB+t3yU+uA6CELqx~ad^klfp* zO~ivWL>8Qw+Q>X~X-Qf7-O6RkE0h?)Hl3#$R8rv|T~4(Q5%* z`;xqOk0HthdPe?bQDh=cmDDT(6_ve=pJYudsCLEJS4PW?6=4-w3*mmzXl8T31j?RS zWUoyRuxqBFf-h-Y;QwmRSU0Syd}TwZaR2G6h5;TQFs9uoofuX;xP-T7WFMIAM2@pk z2)Uwl=A$2|Vnfo!OEF9tpiiiopt7`26njvaDWFj5MfDkpKnP9pf~vF(FO<(~o_*`N z6FM6yRF6F7^b4k>eJjVB*wYe*`{-YnD_{*zH$Yc61@?*v^LV7o#)%z>NT5>xp;;Wf zM)SR>w_FE;VK=`>8swyHaUxp-~5*#pZeqC+^97MUgyuJm_FlRf#tjo`&fxcERwR3 zvS6DICaLXD0ROLy%!r8f9)V)E#~PTP+{IX6i&txlJ-%NV9{IK+;|TAr$=;3k(sTF; zujmsS~8M~K4AGOiq`fq*S4$!>$}Yl8EqA3brPAf7U3#Cl|Kp`vxNdX zym8cb(E1+IkOXtiQq!Mz&}Z(=^f^KTtmDv-94tDX*$cDvyv`f%a$YOF@IVHfC`y4gU)vCyx#EfN4`mV))?{>9L;X0MoA}-- zuyWxTVy6@%p~<21F$4VN^Iy|Qce<1`2!?js!~b@TB-J_XdqC#GBdN0mk00ElWr-E} z<5I&va0Fkkc%5>-*1JV2wOL+TgVv(!Ro%*nEzv@`*(G?uigxDT@)94HZYQkcS@6Ac z!T1h0lT~_Sk49-Z*(d5w(rAovlVC0E~auT80B5!9q51PfP&r)viiK?V$Y3iku z7%bg?Y?;z%TOKItg37wgz{_>&&*Cawe1+|lYoBk)ZB+9o{^_Np)1!FKjoE8iYV~`l zpPwg2e5m2>dn+H&P83rqX%4TE?J#FKIx|Kn)0;P(>T7`NgPY1Na3qHUkE4>3vI$a@8JY#Nzssm6-6 z2Tf|)RULO@z93&8Nf+~xA3Ed1a815iWUy)MCRC18@*>}LZB>SJ#p&*qT%V>o6!+zB zYpPe>m-{dpJ2rWeN~Z!LF2`x<%i75$UnhwO@7;j(LKPm=p`L2{Af}~INALu@0Y&r- z^Dl`Dm}$979#s&CN|CbN#8hkVPmA5vnaK2a&)&8~=UH&E(%ak_Q5?D@6Wxz*0<)& zh>SUMM9lnR#ftRn)6xj}()G$N+8u93a#c85G@=&oI{9LYr#R`E-* z$-pqwpqLks<6ak}p9i7wa+M_J>ZRmtBZ;DwnTiE%hhf z$j`O~2%(R&^EcMUhQV?Pwm61^mCcd-fW73y>qJX~?<9e273ZS$uSamh98gZ`{^5=8&5=QN#4 zBqB0it&5l)(ZLv$r3*h1w9-(Z6oKa#>_U}O)LqIVo%>GUm9$JV)VcH?4?D1*L>=Rg zTBfN;;dr#_M8V=1x6r$9(G6)?O_^OW#;A?_=1j6}AFUQ5Ia8z*uI zpb=X?+{^)i*CNehC6LpU?(*6-ASW&M8$qVRw0;kJ$V^4V5cf&1BJ>j$9xvtu z3VQhLfibWWDH|3tOx31y7}RguMOqxPPYOuqI%$)fw@5s&O9a;}N3Em7a#R(3J9Hzd z9b7RyrugPJL-shS`P>DyIdmoa4zHBr!29p1`9D(bzf<$<%yfUH<{23NEj7>ZM=A6F z*ym$Gd%}jeDPi-B{L-cpXzR4p;H|pR-VcBuIK>rf1*1-msG;0fgStsm;M{b0RCZZ% z%$BZAAJuKxarjCpfx{YaFVMHlL;aO~|LFtn>_t$U`64ItV}z?&Hzex1g;qc8~iB*JKLD)f(=0HxCunIioapH4Ubq!jiI1?(hwzh`ybPwz|EY z*NKnxhK3RmPF8nE=ZA`yOD37ivepz{mce?@g|?`Vx5Ft@QVS2PG#eW@ce$=q(j_Kk zPHDr`baOL@jFNKhX0IdXwh>Q{*XLIwEl19Vp{q$*XJvSkK@}{SiVp`69wVg^Eg}L?j zPTC6rFM-9@t`~c&5#M-Z|F+jo;XdAE6&Oj&ELE2_ZnSLHWZrus$Ll(5@cWsY zQrk8-{AqUO3vSZhic3YBxhN^SbeQz3{!srM!jyS1wou!}n|_(p52L7>En3Ac0}hEP z^P7PuE9s?hA(oG!oAu>8Cj+pxe&rXmU$5{YHP!3(S`RhOe){Z8*18@s2*F+sywd(f z?OG$X4TEALZQ*S%=~%i2gHzYP&WG8sS_62U7FgN&BPo7Ns0hhjaop5y{ zpFCet&an`KfJ6<$zUY zQ|;_RhRV|_T;BKenMRa`Dqxmo4i98|=s;09hiB<%Bf7?^Y1R+(#J65n{~_QbEhxt4 ziE$0irRZwur0_$;Gi|cj@PM;_D=EYo6CGOQ*HL*2$C4jUOf7s`ANxxxGTCKc4|QQf z9;QRUm*tfHR;9xJq={=(g-olgE3_8HwfxnfK(=z`8;m)rY#K>kUHwY*MBiM^k}yjr zHG-`ocHRBGHT`V@A8B1X)R)W#5siEGpEbB25;ONt=Sdq;t3_sBBIP`;%juv^S)4hv(U9 z_OtWju3Kb>k6Y_QrD=gy2=u&IL)A~P<82S1qe90R}1wOjpa zQ4xLlP+t8I!(jofm+ZTSKTEI}oSEJwB28?OpygQw42Yd+{ddcF_E+(6Sl!3)u+sfy_5W zjawc5LDWkp);M$!=hXehp08;ohXit5E+YwW%YuvQr}#4;)(q$YEGeQ5bL02t^RY_f zLSFC58EkoHMH$j>N%J_TJyP@Fe!spCi?YM16e@BM)8CY`JOx!kAQy}Zq$e@W{YlG)D z+pvYh%O8ghB$uigpUnRLX#&m;5Ua)Yaq z{!P!HvD>-&m`oLVY}CE4oO8#?TOKk)6mW6 z)~L1AW}-TseD0Z%`KE)=-}V?&`lET=lwq>{F^C`rSsr&v^V?-RO5{u24_OzhVrz$# zbR{Q8>>-?VNlx&#+Mg112Rsm;+1hJFHuq0|TAEOL*-NYI%9OTtv~$Vgmhxw3y5H%$ zUs?MaW8H?tWs?DXk=&GJwNN+wFv7C1+?04?j1kI1-?v&pc)pC(Os~iZ;q}e5wlU)1 z;R3k`-te2v{JiwSd|Yofg9b(_k4d_$U^U6GyWTol0D>c5%dd&`Zwf|B*XPiQ!G;(dD)ngXL_B%5_;B+!Q6DHa--hjjH(_WH(=8_gZU;xmyX#zm_J zODpnPHzTLFB&}hrLKm`wUfwp=t*uEmiV^fv#Dex2>G40Ph$J!2Rp@cD!v}@De^O!_ zLKsGJ5I5-H09Iu-g`e|Xpo68>giHV(-YMDWCiPF9pB)pA!Hnz@M0OjD0dJ|HLjrr` z2oMs;Pxl_4l5>=d1B91?0=CNMI^XV z{C>`ESEh*MEV%xRN}!%dJ(oK+|H>bu7exsZ^5sl7dZBNyWeoMDbwN4`@)FO$tXBXU zfss(FwpvTfNdFY1u6hf|#wv`D3Oc$=m_qv6_u$4qBrhld?#OYz{VV$#@eY9VBh4Z_ zFnZ}I?|nwCT%}9)12wC8qu4-d(mfmMy*c50xIj-xIPaifEbe)&;xV8RHPXjGRu)UV zPNJd2ZH`rpx*Vb#jeskS6fMaWL!_KAVNwD6p^hUV{p-{7`NUNx_q!>R-Mt~PL)7K) z)Y-IP$UH!*jns1jJN_}u;5!XJ?b|hqXbTV-WR9a)~zWIcbGj){2g_O#pT_k6`w$P#(4gJK&{rKPn z*ISZE3?DMfw_~y`{I{@Qd_mFE&E&06f`tML3RP4QBfIq$?Ia+UT>LGNWDUTWj!B`% zou)-uMW-2i%S*rk@j*xt$QVpv82f{})J4RQj*DG#F1&e>4SMA`k>iOo((X9yU*Dt1 zK8D=oYLkj=($zT_sa-&f^xGl9&QQ%jlQ@`0@~f7UR>wDmgXIi!#~!+)u%#G`3WYI< z=Db#9zQx(q(L{?F2)_23_i4#hQ~Mck?5v2;_kQ@ZLM{$3V;?;*vR=qD=CWJrH}(!r(ofXG6hCzmmGLca*iD`ENA(VL!TvC|&bPBf`6BwUtCDm|xU$ja}p^h&isV$mxH z)#^Wt5!1GibLQNv4k6Ev0?nk^#ArYO)hd-MMQNPT8KK*m0h*CE+q~Zf46Kha0{kuoq zHwqS#bu|e;wN^19*gglCKKu1Rs$&*I&7!3Q6dG>VS5$)v+K zSAKOgPyup}mFMyUA`#O@=T63+0gW8~N77`)qT?VD;v{W)`~;~#4v8&B!j38~CBkFE zg><2)MOVx=7%4(JU<>#;g4cQBZE(vvBcQZpRR@)gx6^~Icj+l$5OPr~hmZvb>Emv+ zfDoKQfb?W?#eT)H4kJ;lV~-z2bH;p`zAPysISnZH!=Jy&Hz5}_EOIMX-C_&WCD+J8 zlh-?xwos0VMURthto`aimgZSgaC6s7_T}wdWx4rwR@(B2IsLd`u{#O6luB%fDj>9h zg8%dsnRTkzNe+DOfSajqGNiA!f4p*6Avk2KtzvMc_#vQ;qY0;cd0}4x0zpm!1~LOD zAPDZZC-;sr5HlEM>5J-WA#b`Vs7!RAtCb5)Erm;$uc`)+9ENQ+!4SIpu2Y#=&ptHu zqZjwuTx*1-Q7~jspuPYzi8yE|x5n1^2kOY|ESql6gaR6i{iE=tMFpZT141WlI4d_U zc%&Ytz?30!vz&Z2^*M z(~4iOjys|I2G0OxmYek-mxQwP|!v#7YYUOILooCh%2cmX;Ute9Lv;RnabXM zVc0+=?0uE`hXm)Jr8F-jrJx;c!y(LlQTK7}cNgLxKR^0haM6!H6IAGGUbt?jl0(z2m8s(z!>5irCh$fctl5Yq- zz-csO%%GG|N7~}qVM|504wn>{e@`9>&9}If65nAU`&DJOa5&GXQdZ6*CEzG&q_nz- z%h)N4M^@eEcFApsmP`|FAcO66#iJ>FT4nXc#U8}tMb?;A zJSMjq1IGA#H)cw#M(p7eD0w(g#T3^xhI0y5-YJzZ5NL>{o|{GNtIdE)UUXaP z`H0Ey=Iws0kc4V$4uz%3(7}Faf~2|_H64+Nh_hE{e*pPOhp(CGnUIsl(}!7Zi4=r- zYRdMYFNTKl>eUcT{W5rQgnn3tgYE!5k3moWk;Qxz;w75$lj~cz;_Oof!^{atB=i(1 zVYpm!e<`uBG{|vYeRvgO%oWNUE?ipOcRC;X_jj-XX~q4jQ2Jsoixg2wU`aAW=Ti0j zlrz}lYz+9^UJ95BNI;4rb#$Qa93v}Ix`5na#ynT~eeyhU#+l}OsHyo>vti17VoMB5 z%D9ce5Q9-FEYrf0(Y|In@GeQoFlem?KAQCc*hkqS_E z=8zGb)O_9oqG`f5nF0X$ zF{7wj1c|w-PqE>;aTH2x?Bv_M0GTiaaybu>0D|k00JoAYwqjo}40HT`3U)zVnS5(K zD9ew&P>QHHz}F1X`|0rR%?Lr(Y$2HY%B}Gew4WZa)`k-T#jti_`eLZTApXp3> zGp--L(jQ=mkO?1ynw714Jo2^{U{v-+O?6!i)!xO4FhRB%-1x%R85tKG=&W-x0_N?3 zG|&UNxM5y^NNY~IlOT&pl=bk|9<b95cuoj~P4AVFPUJNDt}r3J z#8Yh(N&Ga8z7H|}k5QM(p7a6K(l7m^G{krzRd{M)JQvrrF~HLqN33F_2g;M%AiFgx zvugSQP&C_F_<^=%BO+xM@J6b>Zviv0VKu0P`>I#vd4&O^Fnh1CTUO9Eki9y!?vKjIca*3uLhv}#!P#XacBj#1b}yHv91#pDf({bCUzSq;4!OAh4g*--p+ zhsjA}MTN*(k`y0_6z0ApdyU7qQMH^Z^v4{CO5ewY7TdZp#8>$yRT9^kFq6fS=<5O` z3}19tZ?76#G<+?z(l25xGt?DHuR+>`x{*l~YQDyZeT9(@3;Y7W8_l+0mJz^UL2~q1>vnf) zQLmxPN_=#5MqnJq^5F-etxh2#(7Q|J2SbQyK`33Fl2_iBr}qF);ZjFvIl^nWVQ|?D z%O=Euq^W}0R4SiWx9EfEs^`WuZ{~WGsNRRR`9N3#3gqw>JkOMGMtaZ4vl_$Mp`j+4N#K#-qFV_M+fD4U)jywmz416`z*lrFz(mnn`t zJR5PCk~htUP*f4uHPFD+QiOwTK4JvUM8=awwHsod_k9cX5tG!ck>iLufMtBbP`odU z;5WWNMR2I>^sk!QtAL^P?4SmtweFe}%U%nf9KdXyH{UyVvEWxs!$4sTcch9*3(I3W|0ijK}1PwPRwijJ2!Q$SPO zpiG!*k0fbbGQF!H>3Cm=)D8nTV(Am2di;te>sS#FBfu?8u=1dW1aa|zMls5$Q_zh8 z9k(vl;6&a|eyAM4BUnI-&KIbomXGR?RcSz`k+g6&<7fdx2~0E!%5({d=u4r@1A+Yv z4u`?;^FS+y_=hC-|D&k>e@xRiXpcF39g5l5qr9}a%%Drv-dKpzaOaAe{((lA zvZF4h;{$N)>z6ME#b@>ktFy!U3i*zCNoBT^z*6U^tZHm*qHOXF_%Pg%pnGlMet7g~ zWPCA5a4q?AI`E?hH_3Ktwa@$Uk{^RE6d|+v=hxv_3c4Oz@ ztgg}ChUKC;oDSS@)#Mm1&Uu0fN2imU6Gw)&i}P%?Gv$!(_+Vj$^}XqKk7sd~;TIRK?J!4{fs2#x zQza=RoOeT(uh(znJ%vjMMu5aey&S_xL zI_m1+dO00;KBODZ=mPfc_BV&et(wxg+OBsyU})?4*k-S=PD|M&HuQ0vN|{*Z^iI-d zB*=PXjH;l|(57k0*dzX!#I35ap*r>u{kZNsS!db9DB)(>;JVN25@PuLu~+g(*TTE2 zH&iWsM7)!CeiCq9%n|SInN>toifT0JX69sv)l4UQ%&!F8w*l5!&S7+ks>3i_VtL2 zN1U5+tFxuB%u}nQgGN#cSfLjzRhU?e zO-jZHyJmHZP!_=)LTw!F_ z7ES38)!i)+?u8)1cLZ}*k<%q9rO}Y0t_&Med2Z(WiQ_`Z3&@UDDoich|MC_eU-DWiwVDKetT!G!+BF zsk@OLS?jO?~P6;zIbujF#^}?H?@TTMsrWDcGPp%n5}2 z>oKSMCf5K7QB-1{op~#ci@k$+zk&uss%z{Za_8?p3kFhe5L&<>`;%eCVXva!h|*Y{ zwm~92U4Vk@NnBW4rT47-mi=G1_ljFcXv3?`d3$fM)tn97ku!?t#xrPss$dFsi4Z(TJ4l+@LyBU1aybJ}}@&D0bMO8B?a?ZdTc()zVi6nH!K--(VWPG#0RF z#RE&R5u@KCR`{VGWHz!gMMR7f>8sbxJ9>3Mrp!Gwy1f8OF0MXQ?p8`UA$`m#l5Z~R zek1&0XOUnhZMjLJo@~-ggGzlApMwQ1t05ble1G@73PY)w8EWP>gy52lv|Wg?jKEAQ zeka6nh~KM?m7ia0A`GDViAb%|@iy14$mMMSJ*uqI2TcvJ7TT|BkS+u##$R+ixJqvz zQ(@Wkk`sAGgsN7b?O8VbQ?L!e(O76|S3OAE$BHuo->pQWPYMGg(rUF{qSFG8nrfZ9 z1b1tj@o>gRJHOznQ}g8`4p64#d*w&%2`$zV%Ol83g@}=iD#Wboi&Yto33cGEgx09@ zwYfX48%>97Bhw2=bG|U(nS3}U!n`M$(tbRp?RA-?m-Fc>R7*w|XE6#6|3ta9*m)z} zvtys#(%)!A7e>HTbJF`eYxeq1^%My=4sg&b|s+1#NO^Fr}&7Hj!D1m^A ziGC7?v)<3~nOUXjgPIB~D8=WhC*j=VXGm9cwAj}91}aXli#He4d}xJ=KWf}s5DW;l z5t|^IOQDQ|1qpyy?vKt-e=bM6hL#C!VyCGcL?aX+Ok_n9;UVZ&BBR+}Xom3&3G^LM zu@SlD6`Wk^O}@upu6Lk4!l&$vF0L7yzSWpecq)m$*h#z>^l^np$d{=yB?Hj6o$|03 z4_T*7={iyG2l`P;&$Fi5>v9@52G4YaZeSM zkJvEvw1_iaf(JNW<6};kYXnG{l_3&r$dkSBItEQ4Z8Oe7Y{q#@zrNhA?(rugenXU! z7GxaXo1?4BM+WIax~35{%+aL(xD_#cn{iKT!HKP9%)hoJ`VsPDI7tZ_g5}olbk|x7 zY8?UX2Q-_DE;$Aa*)Vfv5CgN&?H4&;>1)ZgjM8Xr;hy+l2^HjbZKbm8; zf=PqiO4#)AIJ}jHB#hX*4VE8-8L#q zNGPn14IWi2d1f%97$;khp77@q0&q)<@IkM8V@$3X*tvcPnUU)uo zVal2Q{*r&AD54>g$|x;48nu2 z;r(R;Gstd5Ik;aFI$o?UcqHgh^qCZd+(BOz#laARO0>tcCTIBr@UN6l4`6Ncm$DSYy?jRZ`+jVO1d-Dw={IVm|dzmDRs&U&B2rC<4%YhPOBY z1jHIuA{Eah<=@ap*nrrOmarbwi^_8&hy-cs^!M+-{?wuk?L-XABD~K^W=j%Obj!Re)R8accWbG_YvsJQpLXnv-+0T$eLhlTGB=#S8;!o#03@3cF)1l4vTso{rO8A| z?|1teDlZs|xMeXho=Ym#D(L~^I{f~@PR(s zkz8F@{{Ea?a7Nlwt=?4nt40M=NK@mg^zy12R#VKkkLb%fJ5-EEVLTqt7bLA4Atx+V z#>5YCaxroDYkw?URTJ`4@s9htqi4M&Laes&B^Vgl1bXTlK99xgNsJ*Mz$%8-*cf2} zRO?0k>gBk=&x1zRl~LMxG^>pza0=sIA`&zz84|D~!8%lLav< z6ZvhFtoj~l>HQxlpg+h*$a_Xs($qqYY^qnl;CEtLA zEGLYUbH)hUf;^}f?ubd`7)ep`<@+&I-gl5`QpE>uR44LxH@%mF>~W>S-Q^a>E}jdG zHwT>GSwT4+>pT!0(EOlWs6^U|U!=TbOq$Mx7DQliLutynJ_ZaMzCsn@^AZ@Lc4Pq{ ztW*}|uEOSIj41pdvB1-v-_G&GRKU+mI1wjFVuzKDD{Gw-d{GFQP{>BgApG;z?*mo2|<+WPl z^pYT7*`QG}f?9Gb6*(@6z%|cplH8-1@1p<$*M$}I_;SzgtbCepz9+iqU9Kn<%Um@$ zu_7pua&Yu?=^x864q8Zkb0aI)9!ho->P65{^_k3-2Q!8<;)45a_ zJ98OWvt`t1XnKbU%}t;c71KZplEC%j;Ddat*HDNKjCZai2PkvX48_7=p?dHam}wA< z7)CmlJFj|}KQi%|n!XY~l71=W4`aA*`kH1oC|QagqQ5#H?pSG4Mj3Hy97F+RAesMu zoOP@X3BcOJfg?Oxs{On$!6insrO3M`3@twU)$RuoAf|%#qZi~26(aq@j zh)w6#^d%^a|CKZ0V zmgXTE{Xc#@ZjunI>4s^BI2Iq2?>rr;OwnBtDjNZ144`~*>*6~Jh_tyj0`P^04Frvn z@cNoq8k#w}84Z0Ra83Q73?72^*>p}>^NCZr{(cmVCEQa4#$|Ff>u2A0Eyg|*IY{J* za9YnSU%(c6i=56^VjA)2ci(e_l)OrS(5o8nVN|}S*;AzK8)jll9NvLN3*20JO3_64 z<{D^GRX?l1ZOijHH6bbsb7X!Hb!anl8}@CZ4iLYiT+bStwq}skxHshJO--l@2+gf- zu5v{lS+UG5NxmXAtm61C3rtXHz-gA4Qw-U1M^%cp1gsN=HB#$SU~T z7kMiDbT?5H(`LP#tHCD)DbKd|HJ~SCk?2IrEt6et=2k8rwaJbiEc8c&{XB^kJS5AT zUJI)DM!2&EiBVUp0uCi^QWvpvJfX86Cu?4ntws))dX~>eF;4iQB3R16@y8Ku?=zM7 zplz>4SN7RXQoxKU$;x<6kutBE03MyGj%Ug@1dClGb?Uu)k_E4|mw?=0U8ptI9G3u1SO+AWTe=;VijFHOts)BvaPmF9ThzZSW=Q#iHF3j5gc*g6AWXsAWbpbfw>$Eyo~ zu~Qh+TZZ^IMF{!`pNDR|f49s5v2SBY1p-z&h6MZ0HnI`hE7!0P_wTvT{eHr}7<(vHZO9?j8G zSYC9Nc=j8nF&JWJfSz%;LfmO?179T4E~o`grroB38ZQO4{~GRh{}Y#d4}&tlsAW#o zB)XG^*;xw@9snZk-T=Nv&yt&D6J#E!q|791L^lVyRtx(8pK*GG0900ek4lInhI(h_ zsf0^kiQA8?3J?&XA(Q<^7@Rz=s!g>Q+b`2(7I`gZrq^!imQ7(){bmO5*`qkMA4IXo z^EyzbS}$5e)fb!=Q>V1UT!-`ZXMq|&FqdZW)Wp6@M}{vujRy=q-g6X3AeZU3AU@3O zr>1^!8N5yEkdeF|Q8o1jS8jsu{P&dRU;RIr85r6Asyt?-|F@JT;~&My|6@wCLsK$t zy%Di7XIheSlfl}6pQZL zn6>}(BxeJ+N85=IgpP{~27Z%*ORlP>sh^#ThfC9n^)FLiYo(%qTl49ag~K4%t8Iby zg^@c&iwgIdp{x&{7Z=B2LwK1~ru38?H%dg(0y%s(Jee`6^qJll>vYx&PNZLbmsZ}@ z0xlUV9b32@GwRmTOIMMPyY{pJjXcF0_wU(1nmUx)J}%)!)Xv&58xc$!7q4udaylSn z+yg{;KI`<-rGX33JWXFO_jI67IIO2BuO5Q%;?J~%cep!?gJ)vuC=eBOJXMcVHyQcC zJKB2CMZfzd3daGH>h;1;uuo620c~7b?EBiglDT_9atWP(2bm?{mT<~wdUVsTUF%pG z(-GU!cR$H@Q|?9c8ZP7KPH60P4Patd^nkPF4zDA(&X+{DPpf2fL8`77PJ~&?gi_2@ zXTY#_aO<%Tyw%8?^$V8EM?Z?Up7XJNkTLUDvPu6%(^2G%+whB~rf#k$D%`o27Hx*o z-@EtrS2FppJ(`Tjv(ip%SZ{+1d?@DSm@cu(p^zvILP30%YaiOoJ(6_B4wL@zWCV@l zh4*RbjWVG%{7Y3L6GB*JBrvn@MS?`8+C5}S1}ebgW$;*hPxQgNHH{d^Y~wj za84;YLbaPhgV3~sX4t3A{>b%Vjxe;T38N`>V72y42J3D2vr>RU_cMI%S=1;M(i#2z zU;<$05Rexot2|=Fw)dWW&Pcaujh+j3F?{lKJVQsGCZUv8jxrIJE7KV!p2t9#+x_fT zK22h^0aUqntn9Hd+2TNA4{A|bLeJAM?LRoheKQ(UyvKu&&-%7EF<^;q{a%*Q6h(PH zcwNt5a(&E1d6Uxe^x3aS55mrzy*wUBc4Byi*>>r_gizpxE4;D96f6wug{U$<3h431 zdx;-?`)sl?bCylzEQlMDD^>N@Oo<)ACMwE-q{;)N$_Ih&L7Xll;;>G-n=V;Ny~Q~s zxnV$>RNw#^htztB1iFEdlMbPM0FMjp+uwfM8a3eMwSMIsDw83l<>lGP{y=gU z(?^rrY4dXJe)|COc#4jYdurTUu_&q6QWgN+JjS0F2}i(RcITuHpmfF|=#a{sf)R7JI%sq{5Q`J-Xp)FF z6q%21%k7^71i-C$q>W8UyP8ovqElArqm_$a8e!!&vMHRx3Nf z=uRq9mP1r=qtQKu!&G9-_k+UVXsT79@J^ayjS(x*2rXmZ#%rGM$H46$QNLLGcfFbj zn@xV%3!a8F^xyCt6P-=i9X_Cf{A&7j@jelriAl2-pCaOjr71cA*X>r6=25LhPmW4# zH6K-F6bHH8rtt?G5=7RL;4F7i&qOG$JblP2JQY;s$j!O*X>U}Zo)e`Q zp3_oB&64#ltjXsAg|oR^x5p-KFvpf{OfHS&q62_*gOs_&)QVY$?o(Q`Ga(=5uIx0{ zeT7=^r68CE0m!KK$8imKV;tEz0{sa1f-^eb7%4kD9 z(8Xssyk!Nzj@kC>hZ=dKETJY1qpX)|EVmf4S6;7i>pWo&?#zP=Db=zCQq77I-5S&} zm|r)A2L_!?vOy_`cz}&;!!C20%$@REt8XjTB~kc5UN*-Z^Y5;6G-rSzS3L!QI%&CO zk?{HBDpBJ{v5({8+h%pL#bG@Fbd%q)nU_RB(XHV{xPEp>ktmGY9pPUwFDz&Dl5*)R zE1aYxjyB7v%&|t)_HdGz#mWo|@_&ObFz&%=<=e$a*0QV8zKWpFAJcbQSHBb0Egr%5 zmo7ieBb3+(Q%NH1jYH_z-Y%wsz>Bq2oHZ9aaSAt?1aBmO3&In-eRqR%UZHp_C+Sx` z@C`&JJY}PG9PN+Q0WW}+`_6K4iZGz$ZPUn+tnp5s0au7Lf@9AIlXD;4qRL)^!l2Bp zQ4oC(kcq9`C57pScHpNvluwsqwnEdpOl@^1A+mI6=2Q}t~MJq^jy}Y0kl*PJ8Fl;6MMGlJ3btv1AXNsK z=F3y^p#X7sDOOqdafH=wp?Sbo8${XoxUA{usI@cY4 zGsjKtLL!O>o2673((Cm*$LshQ)bJiHZWj@Lvg?me&aX!gy!bhfhgyoYb#W0e*S`f`o@ z`79W{ghf%RM>=CtL{>YrxJBz;mNa zSM#~u5AJu)Vzx(qp<*e0g;Hwb=-Gt6ND247MZb*%WN%9+wPdZ+mwK6Z!6 z^((?A>*gxy&@-<4aH<;b;-d4-Wsb`YM-vEJk?G>FbcQcpE(jk;QwVmmeg?6B+U6m` zqHUroVg}7%ls}Dig81cG?0hy^x_;61r&Z_r_pLicE8Llbp9wIc+OnBNXw(>eK2U7_ zglvLZaQ=cISl!yCwj9uzOu%kK$#^-02FN}A@!_lV96InSZ!^SrmLqjj>`5`1)~Tiu z_Tk;g29_O`@=f2`4o{KAxB})lV8y#p8tCbA)qbg^BjYyrZ|@h9v8qE)Vr)H_;fby8 zLb4rG<j0fyr?t=| z_@)gaI%6u`EbJLNSx#M*)tr%N+nU+!y4SJ#EpSXQ=rnPz8Ygl|Hz6FCBJ^w$xKK~N z!;39Fmi5#cHX2lPnmD!2&UkyHFIm6Wj z%r40H!IR;E7>zJD#Z`}eE-*`?GFzPCdVHmK;zKTkK+YdrGK!sn@0E{>_DTts+-;~X z9mQJNT-aF93sc7Bw4+qKcjX1@b-`${1RcnY0MGYkv#zStg_i8Gd+V;9s{LXL+nSPd zbpA}(02L-1sIN__U;T!(<#hE55UqJJ`hKnC<2>YC6ClgN?$!oNj?UKTz zd-l|_;UcM4eDw|`+nA085mPG5Qv!>%^~e#_7-H3FTX}p)dV>Np6*>~(owuB88c_ww z584?Ko+5S^_dwb)&tLqUN#kT7bc=nNU?R!!^q8CD3AMi3x%XroW<2=zTzmLB@f-@P zGAtVMs5r@`6X%ipL|~@{<~fgkmcgM%jlLTCI<-YxDUa4(k)CwKz^exA!S_jt+F0Lv z92_LTbYhDT8(FJ}dug>QCo`UqJg>s7Sh`XKI5oVh3p*;a2@0Src$Xjb?w4gpj4Lh5 zyLWNK-#E#%r7Gf%b~Q5wlBh+#&*O`_@3v7aO09VAlxo!Z2$#LCsnclJ@r=W4iFk3B zgHLK2v>IG;!XO}cnnT#~x((Y+zSY08JsCePDzOQ?tG8e_qRy{Mf8(@uWfyJOOF}Z^ zr~vFLcY?ar=Nu;~QZPfn7deX3#VIO`y?@OTdW=HG3td@UY@tz7Hs4gC;rV#hL1SjV zbSJ2bS+fDNCG2^eZmDt%B{|Xs@9HaC{7~$fgRdTa4 zrj^q-HKvs|HZs>2uyw`Lr2D+hLeGH5%Fdt-MXO-zr0-;mNB23(!Pv$L@ArS1KQ9|Q z+B!QJ8asX}k#n#$R5W(dq?HpA!J}0&c6Iuc6S4gqO#eqpLIH}02Z~nk^ZTc~Bii4BjTmP`)bEf|@bAAWLe>La#IaB{)&p%a${V&7A@CVfY0Ur9_RcL?M_&0nQ{wb~V zzirRo*XAF~_51(q|BMgAUySbH%w3T#lnUFDzV{wHSs_oDR7 zQ1riX`cJg{cFMo;=|7$KzZL#N>tDhQOi&EJG5ha~{>|~fB>w)q|Mf%u=HOrB|C?d| z6vh65#eY5iuT~bHPA#DCX#7X?rd2YxHg=>I|6D(F!$02?e_V#56}B<7H8Qs`#bf-l zVph@F!0C_o$RF&Y|HpFtzWX;ze+RC=WC%vw{DGLBARFAH~Ittz3+q z%nkMbKNS2ErhlLQPmmNo_s%x|tl_U`#bf-F!T*IxzwiFf7%{T{iH+kY>Q=`8Ib?s1 zFfy>RGB$Dg`!0GiN4@Jsqc6N5Yc6OOrOF`o{-lr-`%V;r8GPMFSpQgKg$Dw(nFX=+h-dQxg z(b`_Wd79qm*{;K0anSaMj&59@r>1NM}Ri??3-CNE#1~Ph!l*#i@Eszke>(&GV19(eosd zv2hd}430v`la#p%3o}VG*dTMF=F#xo>GkFLJ0w{3SJPq&Cn{M@v)R1A8X9^^A7DgN zgJ_a!GBobSSES1GY?4;Q?#Z5ITKqebl;Kn@^KwG0o@)ThGEVx}kjUX)L1dPl7uqgz2E zFsDq5D$QpE1&8+KAmDf zhQ)Z?ALVegKvi%t`j5-~hAR|DwVK7E%ydM@5mlK&jpndO7U&sA{Uk4^A;e=?6!0Pp zF-$20lMxJHuu7{;1$w_6XVn_WQjGt87S0}q#>$LU5-zgJ2tISqg+K2zM6tqVflYx; zh7GBkKWa$htk=3`J#dcDb=^%6O0k(?6Jqls;|q3eY`&~RAxncbplzTpb^#P&80O{L zv>WyN*2bMUf|Ca#lGf*1a6sPZc{@fghdCs!Lt}|YH9*Y| z8}F31+oR`j-){W?^)b2`LJRl%;i!*?0kbS=xx%r%Z<%c{W&jx+7<(yT=oUlq@c$TO zTl4a=q23in$p{-Ay{8D3*yzonW46L>UJvNq@N-R~iRakWp|0t4ylXm!13heH4RKKA zcu_sJG|>Xv;R~ccs3&+U#`8mfA)ok7Mb+VlAlQ^Dt5PKgj5ziI$k;r@6#`1IgsgGg zVH~Owe8_BzbA7ji;}VOlCLO%yAeCw0po95T6DK4V5e*!7Jhy`rkZO16fUtw}7?dq7 zeinnmtcl~yg#d5hn2=6eu4Qp7S{ig*Fb`6@CeE?kj(*%`ZT*0sXzK@g!X7(ptX!Ki zI7n{cV6KGBzCmZh|1;>=>$sTVHR!lw%MOk^vaKKR6K!AFzSr>|i0#VYE|yXaxiHu7 z_|)>TKx)tdKfk>moO^y-7jF4ZJFd7bA!^(F9J{Rxm?k carreaux; + /* Index of rectangle ID -> carreaux */ + private Map> carreauxPerRect; + + /* Normalize per individual */ + private static final int INDIVIDU = 0; + /* Normalize per household */ + private static final int MENAGE = 1; + + private static final int NDX_INDIVIDU = -1; + private static final int NDX_MENAGE = -2; + + private static final String[] INSEE_VAR_NAMES = new String[] { "men_surf", + "men_occ5", "men_coll", "men_5ind", "men_1ind", "men_prop", + "men_basr", "ind_age1", "ind_age2", "ind_age3", "ind_age4", + "ind_age5", "ind_age6", "ind_age7", "ind_age8", "ind_srf" }; + private static final int[] INSEE_DATA_TYPE = new int[] { MENAGE, MENAGE, + MENAGE, MENAGE, MENAGE, MENAGE, MENAGE, INDIVIDU, INDIVIDU, + INDIVIDU, INDIVIDU, INDIVIDU, INDIVIDU, INDIVIDU, INDIVIDU, + INDIVIDU }; + + public Insee200mConv() throws Exception { + envelope = new Envelope2D(); + sourceCRS = CRS.decode("EPSG:3035", true); + csvCRS = CRS.decode("EPSG:4326", true); + transform = CRS.findMathTransform(sourceCRS, csvCRS); + } + + public void run() throws Exception { + + /* Compute output raster grid size. */ + sx = (int) Math.round(envelope.width / GRID_SIZE_METERS); + sy = (int) Math.round(envelope.height / GRID_SIZE_METERS); + + /* + * Compute data CRS to grid coordinate transform. TODO Is there a + * simpler way to create a GridGeometry with a given raster size + * (width/height) *and* giving the pixel orientation as a parameter? + */ + float[][] rasterData = createNanFloatArray(sx, sy); + GridCoverage2D gridCoverage = new GridCoverageFactory().create("TEMP", + rasterData, envelope); + GridGeometry2D gridGeometry = gridCoverage.getGridGeometry(); + gridTransform = gridGeometry + .getCRSToGrid2D(PixelOrientation.LOWER_LEFT); + + /* Read the data. */ + readCarreaux(); + readRectangles(); + + if (outputGeotiff) { + outputGeotiff(NDX_INDIVIDU, false); + outputGeotiff(NDX_MENAGE, false); + for (int varIndex = 0; varIndex < INSEE_VAR_NAMES.length; varIndex++) { + outputGeotiff(varIndex, false); + outputGeotiff(varIndex, true); + } + } + + if (outputCsv) { + outputCsv(); + } + } + + private void readCarreaux() throws Exception { + System.out.println("Lecture de la table des carreaux..."); + + carreaux = new HashMap<>(sx * sy); + carreauxPerRect = new HashMap<>(sx * sy / 2); + + DBF tableCar = new DBF("car_m.dbf"); + NumField nbIndField = (NumField) tableCar.getField("ind_c"); + CharField idField = (CharField) tableCar.getField("id"); + CharField idInspireField = (CharField) tableCar.getField("idINSPIRE"); + CharField idRectField = (CharField) tableCar.getField("idk"); + + for (int i = 0; i < tableCar.getRecordCount(); i++) { + tableCar.read(); + + String idInspire = idInspireField.get(); + String strPos = idInspire.substring(15); + DirectPosition sPos = CRSUtils.parseCRS(sourceCRS, strPos); + + DirectPosition2D gridPosition = new DirectPosition2D(); + gridTransform.transform(sPos, gridPosition); + GridCoordinates2D gridCoord = new GridCoordinates2D( + (int) gridPosition.x, (int) gridPosition.y); + + /* + * To be sure we are on the raster grid, we check if the *grid* + * coordinate is in the grid envelope, not the data CRS (there could + * be rounding issues). + */ + if (gridCoord.x >= 0 && gridCoord.x < sx && gridCoord.y >= 0 + && gridCoord.y < sy) { + String idRectangle = idRectField.get(); + + /* In bounds: create new carreau */ + Carreau carreau = new Carreau(); + carreau.id = idField.get(); + carreau.idRect = idRectangle; + carreau.nbIndividus = Float.parseFloat(nbIndField.get()); + carreau.position = sPos; + carreau.gridPosition = gridCoord; + + carreaux.put(carreau.id, carreau); + + /* Add it to rectangle -> carreau index */ + List carreauxForRect = carreauxPerRect + .get(idRectangle); + if (carreauxForRect == null) { + carreauxForRect = new ArrayList<>(16); + carreauxPerRect.put(idRectangle, carreauxForRect); + } + carreauxForRect.add(carreau); + + } + + if (i % 100000 == 0) { + System.out.print("."); + } + } + System.out.println(tableCar.getRecordCount() + " carreaux, " + + carreaux.size() + " dans l'enveloppe."); + tableCar.close(); + } + + private void readRectangles() throws Exception { + System.out.println("Lecture de la table des rectangles..."); + + DBF tableRect = new DBF("rect_m.dbf"); + CharField idRectField = (CharField) tableRect.getField("idk"); + NumField nbIndRectField = (NumField) tableRect.getField("ind_r"); + NumField nbMenRectField = (NumField) tableRect.getField("men"); + NumField[] varField = new NumField[INSEE_VAR_NAMES.length]; + for (int i = 0; i < INSEE_VAR_NAMES.length; i++) { + varField[i] = (NumField) tableRect.getField(INSEE_VAR_NAMES[i]); + } + + int nProc = 0; + for (int i = 0; i < tableRect.getRecordCount(); i++) { + tableRect.read(); + + String idRectangle = idRectField.get(); + List carreauxForRect = carreauxPerRect.get(idRectangle); + if (carreauxForRect != null) { + float nbIndividusRect = Float.parseFloat(nbIndRectField.get()); + float nbMenagesRect = Float.parseFloat(nbMenRectField.get()); + for (Carreau carreau : carreauxForRect) { + float kSum = carreau.nbIndividus / nbIndividusRect; + carreau.nbMenages = nbMenagesRect * kSum; + carreau.varsSummed = new float[varField.length]; + carreau.varsNormalized = new float[varField.length]; + for (int j = 0; j < varField.length; j++) { + carreau.varsSummed[j] = Float.NaN; + carreau.varsNormalized[j] = Float.NaN; + } + for (int j = 0; j < varField.length; j++) { + String varStr = varField[j].get(); + if (varStr.length() > 0) { + float rectVar = Float.parseFloat(varStr); + carreau.varsSummed[j] = rectVar * kSum; + switch (INSEE_DATA_TYPE[j]) { + case INDIVIDU: + carreau.varsNormalized[j] = rectVar + / nbIndividusRect; + break; + case MENAGE: + carreau.varsNormalized[j] = rectVar + / nbMenagesRect; + break; + default: + throw new IllegalStateException(); + } + } + } + } + nProc++; + } + + if (i % 100000 == 0) { + System.out.print("."); + } + } + System.out.println(tableRect.getRecordCount() + " rectangles, " + nProc + + " touchant l'enveloppe."); + tableRect.close(); + } + + private void outputGeotiff(int varIndex, boolean normalized) + throws Exception { + float[][] rasterData = createNanFloatArray(sx, sy); + + for (Carreau carreau : carreaux.values()) { + float var; + switch (varIndex) { + case NDX_INDIVIDU: + var = carreau.nbIndividus; + break; + case NDX_MENAGE: + var = carreau.nbMenages; + break; + default: + var = normalized ? carreau.varsNormalized[varIndex] + : carreau.varsSummed[varIndex]; + break; + } + rasterData[carreau.gridPosition.y][carreau.gridPosition.x] = var; + } + String suffix; + switch (varIndex) { + case NDX_INDIVIDU: + suffix = "ind"; + break; + case NDX_MENAGE: + suffix = "men"; + break; + default: + suffix = INSEE_VAR_NAMES[varIndex] + + (normalized ? "_norm" : "_sum"); + break; + } + String fileName = baseFilename + "_" + suffix + ".tiff"; + GridCoverage2D gridCoverage = new GridCoverageFactory().create( + "INSEE 200m", rasterData, envelope); + writeGeoTiff(fileName, gridCoverage); + } + + private float[][] createNanFloatArray(int sx, int sy) { + float[][] ret = new float[sy][sx]; + for (int iy = 0; iy < ret.length; iy++) { + for (int ix = 0; ix < ret[iy].length; ix++) { + ret[iy][ix] = Float.NaN; + } + } + return ret; + } + + private void writeGeoTiff(String filename, GridCoverage2D gridCoverage) + throws Exception { + System.out.println("Sauvegarde fichier GeoTiff '" + filename + "'..."); + + GeoTiffWriteParams writeParams = new GeoTiffWriteParams(); + writeParams.setCompressionMode(GeoTiffWriteParams.MODE_EXPLICIT); + writeParams.setCompressionType("LZW"); + ParameterValueGroup params = new GeoTiffFormat().getWriteParameters(); + /* + * TODO Force GTRasterTypeGeoKey to RasterPixelIsArea (but this should + * be safe anyway as it is now, as this is the default.) + */ + params.parameter( + AbstractGridFormat.GEOTOOLS_WRITE_PARAMS.getName().toString()) + .setValue(writeParams); + GeoTiffWriter writer = new GeoTiffWriter(new File(filename)); + writer.write(gridCoverage, + params.values().toArray(new GeneralParameterValue[1])); + } + + private void outputCsv() throws Exception { + String filename = baseFilename + ".csv"; + System.out.println("Sauvegarde fichier CSV '" + filename + "'..."); + List headers = new ArrayList(); + headers.add("x"); + headers.add("y"); + headers.add("ind"); + headers.add("men"); + for (String var : INSEE_VAR_NAMES) { + headers.add(var + "_sum"); + headers.add(var + "_norm"); + } + Appendable out = new FileWriter(filename); + CSVPrinter csvPrinter = CSVFormat.DEFAULT.withHeader( + headers.toArray(new String[0])).print(out); + for (Carreau carreau : carreaux.values()) { + DirectPosition targetPosition = transform.transform( + carreau.position, null); + List values = new ArrayList( + INSEE_VAR_NAMES.length * 2 + 4); + values.add(targetPosition.getOrdinate(0)); + values.add(targetPosition.getOrdinate(1)); + values.add(carreau.nbIndividus); + values.add(carreau.nbMenages); + for (int i = 0; i < INSEE_VAR_NAMES.length; i++) { + values.add(carreau.varsSummed[i]); + values.add(carreau.varsNormalized[i]); + } + csvPrinter.printRecord(values); + } + csvPrinter.close(); + } +} diff --git a/src/main/java/com/mecatran/insee/extract/Main.java b/src/main/java/com/mecatran/insee/extract/Main.java new file mode 100644 index 0000000..7969cb3 --- /dev/null +++ b/src/main/java/com/mecatran/insee/extract/Main.java @@ -0,0 +1,107 @@ +/* + * This software is released under the European Union Public Licence (EUPL v.1.1). + * https://joinup.ec.europa.eu/software/page/eupl/licence-eupl + * + * Copyright (c) 2015 Mecatran / DREAL PACA + */ +package com.mecatran.insee.extract; + +import org.geotools.geometry.DirectPosition2D; +import org.geotools.geometry.Envelope2D; +import org.geotools.referencing.CRS; +import org.opengis.referencing.crs.CoordinateReferenceSystem; +import org.opengis.referencing.operation.MathTransform; + +import com.beust.jcommander.JCommander; + +public class Main { + + public static void main(String[] args) throws Exception { + System.out + .println("Convertisseur/extracteur données carroyées 200m INSEE."); + + Parameters params = new Parameters(); + JCommander jCommander = new JCommander(params, args); + + if (params.help) { + jCommander.setProgramName("java -jar insee-200m-extract.jar"); + usage(); + jCommander.usage(); + System.exit(0); + } + if (!params.outputGeotiff && !params.outputCsv) { + System.err + .println("Please use at least one of GeoTiff or CSV output."); + jCommander.usage(); + System.exit(1); + } + + Insee200mConv insee = new Insee200mConv(); + + CoordinateReferenceSystem cmdLineCRS = CRS.decode(params.bboxCRS, true); + CoordinateReferenceSystem dataCRS = CRS.decode("EPSG:3035", true); + MathTransform transform = CRS.findMathTransform(cmdLineCRS, dataCRS); + insee.csvCRS = CRS.decode(params.csvCRS, true); + + DirectPosition2D minCmdLine = new DirectPosition2D(params.left, + params.bottom); + DirectPosition2D maxCmdLine = new DirectPosition2D(params.right, + params.top); + + // Transform to EPSG:3035 CRS + DirectPosition2D minData = new DirectPosition2D(); + DirectPosition2D maxData = new DirectPosition2D(); + transform.transform(minCmdLine, minData); + transform.transform(maxCmdLine, maxData); + + // Align to 200m + minData.x = Math.round(minData.x / Insee200mConv.GRID_SIZE_METERS) + * Insee200mConv.GRID_SIZE_METERS; + minData.y = Math.round(minData.y / Insee200mConv.GRID_SIZE_METERS) + * Insee200mConv.GRID_SIZE_METERS; + maxData.x = Math.round(maxData.x / Insee200mConv.GRID_SIZE_METERS) + * Insee200mConv.GRID_SIZE_METERS - Insee200mConv.GRID_SIZE_METERS; + maxData.y = Math.round(maxData.y / Insee200mConv.GRID_SIZE_METERS) + * Insee200mConv.GRID_SIZE_METERS - Insee200mConv.GRID_SIZE_METERS; + + insee.envelope = new Envelope2D(dataCRS, minData.x, minData.y, + maxData.x - minData.x, maxData.y - minData.y); + System.out.println("Rectangle extraction (EPSG:3035) : " + + insee.envelope); + + insee.envelopeIsSourceCRS = true; + insee.baseFilename = params.outputPrefix; + insee.outputGeotiff = params.outputGeotiff; + insee.outputCsv = params.outputCsv; + insee.run(); + } + + private static void usage() { + String help = "--------------------------------------------------------------\n" + + "Ce programme extrait une partie des données carroyées 200m\n" + + "publiées par l'INSEE, au format GeoTiff et/ou CSV.\n" + + "\n" + + "Ce programme est un logiciel libre (licence EUPL v.1.1)\n" + + "développé par Mecatran pour le compte du DREAL PACA.\n" + + "Pour plus d'information sur la licence, veuillez consulter:\n" + + "https://joinup.ec.europa.eu/software/page/eupl/licence-eupl\n" + + "\n" + + "Les données des rectangles sont fusionnées aux carrés selon\n" + + "la méthode préconisée par l'INSEE (voir méthodologie en ligne).\n" + + "\n" + + "Deux formats sont disponibles: GeoTIFF ou CSV.\n" + + "En GeoTiff, le CRS est celui des données sources (EPSG:3035).\n" + + "En CSV, il est possible de choisir le CRS de sortie.\n" + + "\n" + + "Les variables disponibles sont celles des données publiées\n" + + "(population, ages, revenus, mégages...)\n" + + "Elles sont disponibles en sortie à la fois sous forme de somme\n" + + "sur un carreau (suffixe 'sum') et normalisé par individu ou\n" + + "ménage (suffixe 'norm').\n" + + "\n" + + "Pour plus d'information sur les variables disponibles:\n" + + "http://www.insee.fr/fr/themes/detail.asp?reg_id=0&ref_id=donnees-carroyees\n" + + "--------------------------------------------------------------\n"; + System.out.println(help); + } +} diff --git a/src/main/java/com/mecatran/insee/extract/Parameters.java b/src/main/java/com/mecatran/insee/extract/Parameters.java new file mode 100644 index 0000000..3777d1f --- /dev/null +++ b/src/main/java/com/mecatran/insee/extract/Parameters.java @@ -0,0 +1,44 @@ +/* + * This software is released under the European Union Public Licence (EUPL v.1.1). + * https://joinup.ec.europa.eu/software/page/eupl/licence-eupl + * + * Copyright (c) 2015 Mecatran / DREAL PACA + */ +package com.mecatran.insee.extract; + +import com.beust.jcommander.Parameter; + +public class Parameters { + + @Parameter(names = { "--help" }, description = "Print this help", help = true) + public boolean help; + + @Parameter(names = { "--left" }, description = "Left side of extract bounding box (min lon/X)") + public double left = new Float(-6.0); + + @Parameter(names = { "--right" }, description = "Right side of extract bounding box (max lon/X)") + public double right = new Float(10); + + @Parameter(names = { "--bottom" }, description = "Bottom side of extract bounding box (min lat/Y)") + public double bottom = new Float(40.0); + + @Parameter(names = { "--top" }, description = "Top side of extract bounding box (max lat/Y)") + public double top = new Float(52); + + @Parameter(names = { "--bboxCRS" }, description = "CRS of the command-line bounding box. Please note that the extract bounding box is always a rectangle defined in the source data CRS (EPSG:3035)") + public String bboxCRS = "EPSG:4326"; + + @Parameter(names = { "--csvCRS" }, description = "CRS of the CSV output (default EPSG:4326 aka WGS84)") + public String csvCRS = "EPSG:4326"; + + @Parameter(names = { "-o", "--outputPrefix" }, description = "Output prefix for filenames. Default to 'INSEE_200m'") + public String outputPrefix = "INSEE_200m"; + + @Parameter(names = { "--geotiff" }, description = "Output GeoTiff") + public boolean outputGeotiff = false; + + @Parameter(names = { "--csv" }, description = "Output CSV") + public boolean outputCsv = false; + + // TODO Option to select the variables to output (ages, revenus, ...) +} \ No newline at end of file