2013-02-042013-05-02。61399。1987-。1673-6338201303-0284-04金安1,程承旗2(1.北京大学遥感与地理信息系统研究所,北京100871;2.北京大学工学院,北京100871)、GeoSOT。GeoSOT、R-GeoSOT。GeoSOTP208ADOI10.3969/j.issn.1673-6338.2013.03.015SpatialDataCodingMethodBasedonGlobalSubdivisionGridJINAn1,CHENGChengqi2(1.InstituteofRemoteSensingandGIS,PekingUniversity,Beijing100871,China;2.CollegeofEngineering,PekingUniversity,Beijing100871,China)Abstract:Thecodingrulesofspatialdatafromdifferentsourcesusuallydifferfromeachother,whichisdifficultforspatialretrieval.Inordertosolvethisproblem,acodingmethodforspatialdatabasedonGeoSOTglobalsub-divisiongridwasproposed,itscharacteristicsandadvantageswerediscussed,andtheefficiencyofcodegenerationandspatialretrievalweretested.TheresultshowedthatthecodingalgorithmbasedonGeoSOTwaswithhigheffi-ciency;thespatialretrievalinrationaldatabasebasedonGeoSOTwasmoreefficientthanthatbasedonlatitudeandlongitudecoordinatesandR-treeindex,andthesmallerofthequeryregionthebetteroftheperformance.Keywords:globalsubdivision;GSG;spatialdata;GeoSOT;spatialcode;spatialretrieval、、、1。。GRSWRSGoogleEarthWorldWind。、2。。GSGGlobalSubdivisionGrid3。GeoSOT。1、、4。、5-10。2nGeoSOTGeographicalco-ordinatesSubdividinggridwithOnedimensioninte-gralcodingon2n-Tree11。2013303JournalofGeomaticsScienceandTechnology2013Vol.30No.31.1GeoSOTGeoSOT。GeoSOT。、GeoSOT3360°×180°512°×512°60'64'60″64″。1。GeoSOT0512°×512°。32。。1256°×256°2128°×128°……321/2048″×1/2048″1.5cm。1GeoSOT1.2GeoSOTGeoSOT。GeoSOT0GlGd1d2…dl。1≤l≤32d0123。01232。。2GeoSOTGeoSOT。。GeoSOT。2GIS、、3。GeoSOTMBRMinimumBoundingRectan-gleMBRlevelGeoSOT。。2.1PLBGeo-SOTlevellevel=32Gd1d2…dl…d32。l=0Code=Gl=1Code=Gd1d1PLBd1=0L>0B>01L<0B>02L>0B<03L<0B<0.1l>1l<levelCode。1LBD°M'S″LBGeo-SOTLGeoSOT=LD×64×64+LM×64+LS×2048BGeoSOT=BD×64×64+BM×64+BS×2048.2GeoSOTlevelCellsize。3LGeoSOT/CellsizeBGeoSOT/CellsizeLGeoSOT/Cellsize2BGeoSOT/Cellsize2level-10。4LGeoSOT/Cellsize2BGeoSOT/Cellsize2PGeoSOT4。BGeoSOT/Cellsize2=b1b2…blevel-1LGeoSOT/Cellsize2=l1l2…llevel-1PGeoSOT4=p1p2…plevel-1。pi=bi×2+li1≤i≤level-1。5PGeoSOT4Gd1。2.2582303MBRGeoSOT“”MBRlevelGeoSOTlevellevel1、24。。1MBRleftrightdownup。2leftrightdownupGeoSOTleftGeoSOTrightGeoSOTdownGeoSOTupGeoSOTMBRsize=MAXrightGeoSOTleftGeoSOTupGeoSOTdownGeoSOTGeoSOTlevellevel=32-log2sizeMAX。3MBR2.1level。3。3GeoSOT2.3GeoSOT、。1GeoSOT。GeoSOTGeo-SOT。2GeoSOT。GeoSOT。3GeoSOT。GeoSOTCGCS2000GeoSOT。3R-。GeoSOT、R-。CPUIntelCorei5-25202.50GHz8Gb5400r/min。11001°×1°。。1leftdownleft-180°180°down-90°90°。2lengthwidthlengthwidth1°。3leftlengthright=left+lengthright>180°right=180°up=down+widthup>90°up=90°。4leftrightdownupcodes。10031.775s0.031775ms。21MySQL6822013B+R-。、R-1°×1°2°×2°3°×3°…10°×10°1000。31000GeoSOT。43。a3b0~300ms4GeoSOTR-4250msR-2150msGeo-SOT。1°×1°2ms8°×8°R-。3°×3°5°×5°9°×9°GeoSOTGeoSOT2n。R-B+。MBR、、、4。。4GeoSOT。GeoSOTGeoSOTR-。1.J.2012426805-813.2.J.2010275382-386.3.J.2009285657-662.4.GeoDNAJ.2009454663-667.5WHITEDENISKIMERLINGJONAOVERTONvSCOTTW.CartographicandGeometricComponentsofaGlobalSamplingDesignforEnvironmentalMonitoringJ.CartographyandGe-ographicInformationScience19921915-22.6DUTTONGEOFFREY.EncodingandHandlingGeospatialDa-tawithHierarchicalTriangularMeshesC∥Proceedingsof7thInternationalSymposiumonSpatialDataHandling.DelftNetherlands199615-28.7FEKETEGYORGYTREINISHLLOYDA.SphereQuadtreesANewDataStructuretoSupporttheVisualizationofSpherical-lyDistributedDataC∥ProceedingsofExtractingMeaningfromComplexDataProcessingDisplayInteraction.SantaClaraUnitedStates1990242-253.8.STQIED.200424-28.9.G2ISJ.20075625-29.10.J.200824411-15.11.M.201251-54.782303