From 0ebf57e79874686ec7b4de2b8381c29c8595962e Mon Sep 17 00:00:00 2001
From: Karim Bahgat <Karim.bahgat.norway@gmail.com>
Date: Tue, 4 Aug 2015 16:46:47 +0200
Subject: [PATCH] New subdivisioning of package, and added red warnings during
 testing

---
 pycrs/__init__.py                    |   2 +-
 pycrs/__init__.pyc                   | Bin 269 -> 265 bytes
 pycrs/crstable.txt                   |   1 -
 pycrs/elements/__init__.py           |   6 +++
 pycrs/elements/__init__.pyc          | Bin 0 -> 408 bytes
 pycrs/{ => elements}/datums.py       |   0
 pycrs/{ => elements}/datums.pyc      | Bin
 pycrs/{ => elements}/directions.py   |   0
 pycrs/{ => elements}/directions.pyc  | Bin
 pycrs/{ => elements}/ellipsoids.py   |   0
 pycrs/{ => elements}/ellipsoids.pyc  | Bin
 pycrs/{ => elements}/parameters.py   |   0
 pycrs/{ => elements}/parameters.pyc  | Bin
 pycrs/{ => elements}/projections.py  |   0
 pycrs/{ => elements}/projections.pyc | Bin
 pycrs/{ => elements}/units.py        |   0
 pycrs/{ => elements}/units.pyc       | Bin
 pycrs/parser.py                      |  78 ++++++---------------------
 pycrs/parser.pyc                     | Bin 15345 -> 15294 bytes
 pycrs/{webscrape.py => utils.py}     |  11 ++--
 pycrs/{webscrape.pyc => utils.pyc}   | Bin 2357 -> 2590 bytes
 tester.py                            |  62 +++++++++++----------
 22 files changed, 65 insertions(+), 95 deletions(-)
 delete mode 100644 pycrs/crstable.txt
 create mode 100644 pycrs/elements/__init__.py
 create mode 100644 pycrs/elements/__init__.pyc
 rename pycrs/{ => elements}/datums.py (100%)
 rename pycrs/{ => elements}/datums.pyc (100%)
 rename pycrs/{ => elements}/directions.py (100%)
 rename pycrs/{ => elements}/directions.pyc (100%)
 rename pycrs/{ => elements}/ellipsoids.py (100%)
 rename pycrs/{ => elements}/ellipsoids.pyc (100%)
 rename pycrs/{ => elements}/parameters.py (100%)
 rename pycrs/{ => elements}/parameters.pyc (100%)
 rename pycrs/{ => elements}/projections.py (100%)
 rename pycrs/{ => elements}/projections.pyc (100%)
 rename pycrs/{ => elements}/units.py (100%)
 rename pycrs/{ => elements}/units.pyc (100%)
 rename pycrs/{webscrape.py => utils.py} (92%)
 rename pycrs/{webscrape.pyc => utils.pyc} (66%)

diff --git a/pycrs/__init__.py b/pycrs/__init__.py
index de33765..4e0ec82 100644
--- a/pycrs/__init__.py
+++ b/pycrs/__init__.py
@@ -1,4 +1,4 @@
 from . import loader
 from . import parser
-from . import webscrape
+from . import utils
 
diff --git a/pycrs/__init__.pyc b/pycrs/__init__.pyc
index 65e754e04709eefd01119fd50a713f8f499347e8..ee49d64b898f5ca769a8df3711e477aaa95d582f 100644
GIT binary patch
delta 28
jcmeBW>SSVP{>;nOe)_;f_G(U61_p-GlFXdqiT!l|d-4e}

delta 32
ncmeBV>SbbQ{>;m@w0put_G(^E1_p-m)TH9%qQrvKiQRPooCpeT

diff --git a/pycrs/crstable.txt b/pycrs/crstable.txt
deleted file mode 100644
index d18a776..0000000
--- a/pycrs/crstable.txt
+++ /dev/null
@@ -1 +0,0 @@
-codetype	code	proj4	ogcwkt	esriwkt
diff --git a/pycrs/elements/__init__.py b/pycrs/elements/__init__.py
new file mode 100644
index 0000000..a17b4e7
--- /dev/null
+++ b/pycrs/elements/__init__.py
@@ -0,0 +1,6 @@
+from . import datums
+from . import ellipsoids
+from . import projections
+from . import units
+from . import directions
+from . import parameters
diff --git a/pycrs/elements/__init__.pyc b/pycrs/elements/__init__.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..28e8a4ac3e61d85bd130ce5929f04c6961357a38
GIT binary patch
literal 408
zcmY+A&q~BF5XL9{x66v)GkEH0A3#Kf74#(1WiKI6YMVpWqzO%W=-YYsaeM)1k_v0W
z_l4hNGLt0#x!&D=f1P?*M~3SY$`5F$fq@1v3^)TuoJG!nXNmL3Iq*Dj5xD?fBrYSD
zz{|u{<O+C&yn$oA#YT#u)nuWfzFnA&PrO5xS}W(R8rIkq#@uAjvJc9RYp}%NJgRVU
zZcx*|6sqOAsWURj$-mVb<bZ-E8uf4J1Dc$yAfu`^>cadZ;^9u5JQfk3%2@Gi``O4b
vc=4jb>--_!mqT+DZs{>YYMJ)hR^k2HcJG!zuf8{Sn6-SMzMk%O!~W(k?<h?D

literal 0
HcmV?d00001

diff --git a/pycrs/datums.py b/pycrs/elements/datums.py
similarity index 100%
rename from pycrs/datums.py
rename to pycrs/elements/datums.py
diff --git a/pycrs/datums.pyc b/pycrs/elements/datums.pyc
similarity index 100%
rename from pycrs/datums.pyc
rename to pycrs/elements/datums.pyc
diff --git a/pycrs/directions.py b/pycrs/elements/directions.py
similarity index 100%
rename from pycrs/directions.py
rename to pycrs/elements/directions.py
diff --git a/pycrs/directions.pyc b/pycrs/elements/directions.pyc
similarity index 100%
rename from pycrs/directions.pyc
rename to pycrs/elements/directions.pyc
diff --git a/pycrs/ellipsoids.py b/pycrs/elements/ellipsoids.py
similarity index 100%
rename from pycrs/ellipsoids.py
rename to pycrs/elements/ellipsoids.py
diff --git a/pycrs/ellipsoids.pyc b/pycrs/elements/ellipsoids.pyc
similarity index 100%
rename from pycrs/ellipsoids.pyc
rename to pycrs/elements/ellipsoids.pyc
diff --git a/pycrs/parameters.py b/pycrs/elements/parameters.py
similarity index 100%
rename from pycrs/parameters.py
rename to pycrs/elements/parameters.py
diff --git a/pycrs/parameters.pyc b/pycrs/elements/parameters.pyc
similarity index 100%
rename from pycrs/parameters.pyc
rename to pycrs/elements/parameters.pyc
diff --git a/pycrs/projections.py b/pycrs/elements/projections.py
similarity index 100%
rename from pycrs/projections.py
rename to pycrs/elements/projections.py
diff --git a/pycrs/projections.pyc b/pycrs/elements/projections.pyc
similarity index 100%
rename from pycrs/projections.pyc
rename to pycrs/elements/projections.pyc
diff --git a/pycrs/units.py b/pycrs/elements/units.py
similarity index 100%
rename from pycrs/units.py
rename to pycrs/elements/units.py
diff --git a/pycrs/units.pyc b/pycrs/elements/units.pyc
similarity index 100%
rename from pycrs/units.pyc
rename to pycrs/elements/units.pyc
diff --git a/pycrs/parser.py b/pycrs/parser.py
index b94176e..3d418ec 100644
--- a/pycrs/parser.py
+++ b/pycrs/parser.py
@@ -6,12 +6,12 @@
 # especially: http://fossies.org/windows/misc/saga_2.1.4_x64.zip/saga_2.1.4_x64/saga_prj.dic
 # also: http://saga.sourcearchive.com/documentation/2.0.7plus-pdfsg-2/crs__base_8cpp_source.html
 
-from . import datums
-from . import ellipsoids
-from . import parameters
-from . import units
-from . import projections
-from . import webscrape
+from .elements import datums
+from .elements import ellipsoids
+from .elements import parameters
+from .elements import units
+from .elements import projections
+from . import utils
 
 def from_epsg_code(code):
     """
@@ -26,7 +26,7 @@ def from_epsg_code(code):
     """
     # must go online (or look up local table) to get crs details
     code = str(code)
-    proj4 = webscrape.crscode_to_string("epsg", code, "proj4")
+    proj4 = utils.crscode_to_string("epsg", code, "proj4")
     crs = from_proj4(proj4)
     return crs
 
@@ -43,7 +43,7 @@ def from_esri_code(code):
     """
     # must go online (or look up local table) to get crs details
     code = str(code)
-    proj4 = webscrape.crscode_to_string("esri", code, "proj4")
+    proj4 = utils.crscode_to_string("esri", code, "proj4")
     crs = from_proj4(proj4)
     return crs
 
@@ -60,7 +60,7 @@ def from_sr_code(code):
     """
     # must go online (or look up local table) to get crs details
     code = str(code)
-    proj4 = webscrape.crscode_to_string("sr-org", code, "proj4")
+    proj4 = utils.crscode_to_string("sr-org", code, "proj4")
     crs = from_proj4(proj4)
     return crs
 
@@ -290,40 +290,6 @@ def _from_wkt(string, wkttype, strict=False):
                             item = itemclass(value)
                             params.append(item)
                             
-##                        if subcontent[0] == '"Azimuth"':
-##                            item = parameters.Azimuth(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"scale_factor"':
-##                            item = parameters.ScalingFactor(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"latitude_of_origin"':
-##                            item = parameters.LatitudeOrigin(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"standard_parallel_1"':
-##                            item = parameters.LatitudeFirstStndParallel(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"standard_parallel_2"':
-##                            item = parameters.LatitudeSecondStndParallel(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"Standard_Parallel_1"':
-##                            item = parameters.LatitudeTrueScale(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"Central_Meridian"':
-##                            item = parameters.CentralMeridian(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"Longitude_Of_Center"':
-##                            item = parameters.LongitudeCenter(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"false_easting"':
-##                            item = parameters.FalseEasting(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"false_northing"':
-##                            item = parameters.FalseNorthing(subcontent[1])
-##                            params.append(item)
-##                        elif subcontent[0] == '"satellite_height"':
-##                            item = parameters.SatelliteHeight(subcontent[1])
-##                            params.append(item)
-                            
             # find unit
             for part in content:
                 if isinstance(part, tuple):
@@ -338,10 +304,6 @@ def _from_wkt(string, wkttype, strict=False):
             else:
                 raise Exception("The specified unit name could not be found")
 
-##            if subcontent[0].strip('"') == "Meters":
-##                unittype = parameters.UnitType(units.Meter())
-##            elif subcontent[0].strip('"') == "degree":
-##                unittype = parameters.UnitType(units.Degree())
             metmult = parameters.MeterMultiplier(value)
             linunit = parameters.Unit(unittype, metmult)
             
@@ -408,10 +370,6 @@ def _from_wkt(string, wkttype, strict=False):
                 unittype = parameters.UnitType(unit)
             else:
                 raise Exception("The specified unit name could not be found")
-##            if subcontent[0].strip('"') == "Meters":
-##                unittype = parameters.UnitType(units.Meter())
-##            elif subcontent[0].strip('"') == "degree":
-##                unittype = parameters.UnitType(units.Degree())
             metmult = parameters.MeterMultiplier(value)
             angunit = parameters.AngularUnit(unittype, metmult)
             
@@ -445,7 +403,7 @@ def from_proj4(string, strict=False):
     # parse arguments into components
     # use args to create crs
 
-    # SLIGTHLY MESSY STILL, CLEANUP..
+    # TODO: SLIGTHLY MESSY STILL, CLEANUP..
 
     params = []
     partdict = dict([part.split("=") for part in string.split()
@@ -458,9 +416,9 @@ def from_proj4(string, strict=False):
         # first, get the default proj4 string of the +init code
         codetype, code = partdict["+init"].split(":")
         if codetype == "EPSG":
-            initproj4 = webscrape.crscode_to_string("epsg", code, "proj4")
+            initproj4 = utils.crscode_to_string("epsg", code, "proj4")
         elif codetype == "ESRI":
-            initproj4 = webscrape.crscode_to_string("esri", code, "proj4")
+            initproj4 = utils.crscode_to_string("esri", code, "proj4")
 
         # make the default into param dict
         initpartdict = dict([part.split("=") for part in initproj4.split()
@@ -512,7 +470,7 @@ def from_proj4(string, strict=False):
         coeffs = partdict["+towgs84"].split(",")
         datumshift = parameters.DatumShift(coeffs)
 
-        # if no datum, use ellips + towgs84 params to create the correct datum
+        # TODO: if no datum, use ellips + towgs84 params to create the correct datum
         # ...??
 
     # COMBINE DATUM AND ELLIPS
@@ -587,12 +545,10 @@ def from_proj4(string, strict=False):
         # create proj param obj
         proj = parameters.Projection(projdef)
 
-##        # collect param objects
-##        for key,val in partdict.items():
-##            objclass = parameters.find(key, "proj4", strict)
-##            if objclass:
-##                obj = objclass(val)
-##                params.append(obj)
+        # Because proj4 has no element hierarchy, using automatic element find() would
+        # ...would not be very effective, as that would need a try-fail approach for each
+        # ...element type (parameter, projection, datum, ellipsoid, unit).
+        # ...Instead load each element individually.
 
         # CENTRAL MERIDIAN
 
diff --git a/pycrs/parser.pyc b/pycrs/parser.pyc
index eb1c1cfc77835604fcb116dc52ab8104fe20d171..6441d4a4e492097d0a617ea251a30fb945af6e1f 100644
GIT binary patch
delta 498
zcmZvYO-lk%6o%h3Iy26fV@4mk5FsR&6>1+7ZCbPos+|iV7>mlRbi|Dnq~Iz_?muwR
zIv83+&?YFNUBAF(+vo?>og41ta_)KVIS=<<&Uf><sf2zkz46ggPL32Dt0<}<S49UX
zM&L#SXmB+FI?4uI9UVf^6l@5#1e=1Rf-S+e;Hco3U|aB5Z*Zo`02<ywrS6yLb!IAL
zkxeFF2Bj!L3TYIVOo=&}4I;mt`66ohmDv}K?)|xWL#8fUEu53a28C(%S-7BgwzZhP
zyn_cE5?ljqLTwFRD5$guZ3$HE4%=bT-iH?!6&bD^R9poKo~T+HimDKA;VGlV{7{Sp
zVy^ZGPaSQV@c2lzK)d2E>DbAdd<OMBw!YM-H0!Bv>|yyjs`29Xc(C!0gYB-|EIdiv
xo8%{0HpL_!vQ0^z|6?;7=PiCaP_C4Z%T=$=QcjYREay1XXQxg}o-lur{{tf)VI}|o

delta 522
zcmZvYJxc>Y5Qg7La+lm!5(xoKV?+*9*jSj02DK5av^cwf2KAz%m}qj6#xFDoHj3;&
zP_PwAZzHC(7O}81t@s0UZ$xY?X7-t#X9s5Qn;%Uz{AuY<=Tfb00$VCdD#(kX1(c%j
zVhk~O8iR(44o^pmP%@YeW|P@uwwNvEICGpi!JJ@DUcT*E8h~bHxv^4T+OD)q)EOjO
z4A^N>g=5c<(mu4iL^tjAp$Cm#+QTD;LU*E&-zP)7=6l6c{(!n-*2$(G;R724FN709
zeG-0HEIDc2<Z=OiL>xKkLxsyz_zL?fJT<tH+xE$Ws*As4Vl3XmkNz)?{)!rl&+uda
ziza+6_%AvBm`r#E=5v-@yHOF@!j-twI$9$l>)m>FT%3)a_v^B7Q$dm&!8c`Et=G27
zwY8P9l!QHgG5sEr12VD$D?!Q7Ut}cz#9XdNw%tH-1JR`*$#-zWL5c_}7U^0{6!(-V
H^9#;5jx}h+

diff --git a/pycrs/webscrape.py b/pycrs/utils.py
similarity index 92%
rename from pycrs/webscrape.py
rename to pycrs/utils.py
index 6b8e7e6..f27a6eb 100644
--- a/pycrs/webscrape.py
+++ b/pycrs/utils.py
@@ -4,7 +4,11 @@ import re
 def build_crs_table(savepath):
     """
     Build crs table of all equivalent format variations by scraping spatialreference.org.
-    Takes a while.
+    Saves table as tab-delimited text file.
+    NOTE: Might take a while.
+
+    Arguments:
+    - *savepath*: The absolute or relative filepath to which to save the crs table, including the ".txt" extension. 
     """
     # create table
     outfile = open(savepath, "wb")
@@ -65,6 +69,9 @@ def crscode_to_string(codetype, code, format):
     - *codetype*: "epsg", "esri", or "sr-org".
     - *code*: The code.
     - *format*: The crs format of the returned string. One of "ogcwkt", "esriwkt", or "proj4", but also several others...
+
+    Returns:
+    - Crs string in the specified format. 
     """
     link = 'http://spatialreference.org/ref/%s/%s/%s/' %(codetype,code,format)
     result = urllib2.urlopen(link).read()
@@ -87,7 +94,5 @@ def crscode_to_string(codetype, code, format):
 ##    return result
     
 
-##if __name__ == "__main__":
-##    build_crs_table("crstable.txt")
 
 
diff --git a/pycrs/webscrape.pyc b/pycrs/utils.pyc
similarity index 66%
rename from pycrs/webscrape.pyc
rename to pycrs/utils.pyc
index 79bc85609324610154fd19198690a10203f6ecb7..497b14421f05f6374adfc902cecf2c14c74e951c 100644
GIT binary patch
delta 368
zcmY+9u};G<6h)ID(h%qXLP%ZEqaz|&x}b;w2^ayXbn50LPGjLFRs526>jwancmAQO
z*qQkShJFM);>BXA2m4vq=icj&qxWNX=cntxe68G}<Ak;zJM8v}HCzF*SaM)y2vRP^
zR3J|Q`2ryml7~LfdUzb&MQ|@BITG<q0L*hzb+^-7Dl(ce%+Dhj<s`+blr?fFE2y|+
zC>n08+Bbk@OD1M$K|n&1>c7YZh$4}7DvF5_`XMf`4^)O1N|Z$irVcumWzA#s=t$Pz
z&aYz+^u`N!YZOxHL-)+~^}v0#$LmiQ>`Ma;Io4I7B6A7ghB6gaZIUQVn40lKWP*O5
jmDQAC2zt|Lww&L*-u^X(Dwj6vK(D-DOF!8=wXgmFI8}3+

delta 157
zcmbOyvQ>zU`7<vU=dz7#K}<|V43ndn68S<Bvr~%|5*5lbGILV(Hn%bhFf*D@KEtWU
zS)Q6yoLrPxkUIG!n<69I<kxIHjEb8p*cBKjU*sr-DAI*1(krMe5oBOsNG>W)&QD2=
iFUgNDE-A{)OSb_E7PA2f2F}UzICU`$|IaDMXaN8b4=~~Y

diff --git a/tester.py b/tester.py
index 66ad539..bad2fa1 100644
--- a/tester.py
+++ b/tester.py
@@ -1,5 +1,6 @@
 import pycrs
 import traceback
+import warnings
 
 
 
@@ -67,15 +68,28 @@ def render_world(crs):
 # Source string generator
 def sourcestrings(format):
     # TODO: now bunch of randoms, instead add only most commonly used ones
-    yield pycrs.webscrape.crscode_to_string("esri", 54030, format)
-    yield pycrs.webscrape.crscode_to_string("sr-org", 7898, format)
-    yield pycrs.webscrape.crscode_to_string("sr-org", 6978, format)
-    yield pycrs.webscrape.crscode_to_string("epsg", 4324, format)
-    yield pycrs.webscrape.crscode_to_string("sr-org", 6618, format)
-    yield pycrs.webscrape.crscode_to_string("sr-org", 22, format)
-    yield pycrs.webscrape.crscode_to_string("esri", 54031, format)
+    yield pycrs.utils.crscode_to_string("esri", 54030, format)
+    yield pycrs.utils.crscode_to_string("sr-org", 7898, format)
+    yield pycrs.utils.crscode_to_string("sr-org", 6978, format)
+    yield pycrs.utils.crscode_to_string("epsg", 4324, format)
+    yield pycrs.utils.crscode_to_string("sr-org", 6618, format)
+    yield pycrs.utils.crscode_to_string("sr-org", 22, format)
+    yield pycrs.utils.crscode_to_string("esri", 54031, format)
     # add more...
 
+    # Misc other crs for testing
+    #crs = pycrs.utils.crscode_to_string("esri", 54030, "proj4")
+    #crs = pycrs.utils.crscode_to_string("sr-org", 6978, "proj4")
+    #crs = pycrs.parser.from_sr_code(7898)
+    #crs = pycrs.parser.from_epsg_code(4324)
+    #crs = pycrs.parser.from_sr_code(6618)
+    #crs = pycrs.parser.from_sr_code(22)
+    #crs = pycrs.parser.from_esri_code(54031)
+    #proj4 = "+proj=longlat +ellps=WGS84 +datum=WGS84"
+    #proj4 = "+proj=aea +lat_1=24 +lat_2=31.5 +lat_0=24 +lon_0=-84 +x_0=400000 +y_0=0 +ellps=GRS80 +units=m +no_defs "
+    #proj4 = "+proj=larr +datum=WGS84 +lon_0=0 +lat_ts=45 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs"
+    #proj4 = "+proj=nsper +datum=WGS84 +ellps=WGS84 +lon_0=-60 +lat_0=40 +h=2000000000000000000000000"
+
 
 
 
@@ -83,33 +97,23 @@ def sourcestrings(format):
 # Testing format outputs
 def testoutputs(crs):
     print("To:\n")
-    try: result = crs.to_ogc_wkt()
-    except: result = traceback.format_exc()
-    print("ogc_wkt: %s \n" % result)
+    print("ogc_wkt:\n")
+    try: print(crs.to_ogc_wkt()+"\n")
+    except: warnings.warn(traceback.format_exc())
 
-    try: result = crs.to_esri_wkt()
-    except: result = traceback.format_exc()
-    print("esri_wkt: %s \n" % result)
+    print("esri_wkt:\n")
+    try: print(crs.to_esri_wkt()+"\n")
+    except: warnings.warn(traceback.format_exc())
           
-    try: result = crs.to_proj4()
-    except: result = traceback.format_exc()
-    print("proj4: %s \n" % result)
+    print("proj4:\n")
+    try: print(crs.to_proj4()+"\n")
+    except: warnings.warn(traceback.format_exc())
 
 
 
 
-# Misc crs for testing
-#crs = pycrs.webscrape.crscode_to_string("esri", 54030, "proj4")
-#crs = pycrs.webscrape.crscode_to_string("sr-org", 6978, "proj4")
-#crs = pycrs.parser.from_sr_code(7898)
-#crs = pycrs.parser.from_epsg_code(4324)
-#crs = pycrs.parser.from_sr_code(6618)
-#crs = pycrs.parser.from_sr_code(22)
-#crs = pycrs.parser.from_esri_code(54031)
-#proj4 = "+proj=longlat +ellps=WGS84 +datum=WGS84"
-#proj4 = "+proj=aea +lat_1=24 +lat_2=31.5 +lat_0=24 +lon_0=-84 +x_0=400000 +y_0=0 +ellps=GRS80 +units=m +no_defs "
-#proj4 = "+proj=larr +datum=WGS84 +lon_0=0 +lat_ts=45 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs"
-#proj4 = "+proj=nsper +datum=WGS84 +ellps=WGS84 +lon_0=-60 +lat_0=40 +h=2000000000000000000000000"
+#############################################################################
+
 
 
 
@@ -130,7 +134,7 @@ for wkt in sourcestrings("ogcwkt"):
         testoutputs(crs)
         
     except:
-        print(traceback.format_exc()+"\n")   
+        warnings.warn(traceback.format_exc()+"\n")
 
 #render_world(crs)