ó f÷ Yc@suddlmZddlmZddlZddlZdd d„ƒYZdd d„ƒYZdd d „ƒYZdS( iÿÿÿÿ(tPolygon(tdistanceNtGridcBsMeZdZd„Zd„Zd„Zd„Zd„Zd„Zd„Z RS(sW x y are gps cor_x, cor_y are the grid number cordinate in two dimensions cCsr||_||_||_||_|||_|||_||_||_|j|jf|_g|_ dS(N( tstart_xtstart_ytstep_xtstep_ytend_xtend_ytcor_xtcor_ytgridCoordinatetintersectTransRegion(tselfRRRRR R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt__init__ s        cCs||_||_dS(N(R R (R R R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytset_corrdinatess cCs||_||_dS(N(RR(R RR((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt set_start_x_ys cCsD||jkr@||jkr@||jkr@||jkr@tStS(N(RRRRtTruetFalse(R ttest_xttest_y((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytisInGrids<cCs|jS(N(R (R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytgetGridCoordinate scCs&||jkr"|jj|ƒndS(N(R tappend(R t transRegion((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytaddIntersectTransRegion"scCs|jS(N(R (R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytgetIntersectTransRegion%s( t__name__t __module__t__doc__RRRRRRR(((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyRs     t TransRegioncBsPeZd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Z RS(cCs||_||_g|_dS(N(tverticest region_idtgrids(R RR ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyR*s  cCs|jS(N(R (R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt getRegionID/scCs&||jkr"|jj|ƒndS(N(R!R(R tgrid((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytaddGrid1scCs|jS(N(R!(R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytgetGrids4scCs|jS(N(R(R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt getVertices6scCsÆt|ƒ}t}|d}x£td|ƒD]’}||d|k||d|kkr¸|||d||d|||d||d||d||dkr¸| }n|}q,W|S(Nii(tlenRtrange(R Rttestxttestytnverttctjti((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytpnpoly8s  +?  cCs|j|j|d|dƒS(Nii(R/R(R tpoint((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt pointInRegionHscCs|jS(N(R (R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytgetGeoIDKs( RRRR"R$R%R&R/R1R2(((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyR)s       t RegionHandlercBskeZd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Z d„Z d „Z d „Z RS( cCs8|d|_|d|_|d|_|d|_dS(Ntmin_xtmin_ytmax_xtmax_y(R4R5R6R7(R tminmax((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyRPs   c CsGg|_||_||_|j|jt|jdƒ|_|j|jt|jdƒ|_ |j|jtdƒ|_ |j|j tdƒ|_ x£t d|jƒD]}g}|j |j|}x\t d|jƒD]H}|j |j |}t |||j|j ||ƒ}|j|ƒqãW|jj|ƒq°WdS(Niii(R!t grid_widtht grid_heightR6R4tfloatt x_grid_stepR7R5t y_grid_stept grid_start_xt grid_start_yR(RR( R R9R:tiitone_rowtcurrent_start_xtjjtcurrent_start_ytone_grid((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt buildGridsVs   $$!cCsTg}x>|D]6}|d}|d}t||ƒ}|j|ƒq W||_dS(Nt geo_arraytgeo_id(RRt transRegions(R tout_edget all_regionst one_regionRGRH((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytbuildTransRegionsgs   cCs^xW|jD]L}|jƒ}x7|D]/}|j|ƒ}|j|ƒ|j|ƒq#Wq WdS(N(RIR&t findPointGridRR$(R RLRtverticeRE((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytmapRegionToGridps    cCs|jS(N(RI(R ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytgetTransRegions}scCs|j|ƒ}|jƒS(N(RNR(R R0t target_grid((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytregionCandiates€scCs]|d}|d}t||j|jƒ}t||j|jƒ}|j||}|S(Nii(tintR>R<R?R=R!(R R0txtytgrid_xtgrid_yRR((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyRN„s   cCs-|GHx!|D]}|j|ƒr |Sq WdS(N(R1tNone(R R0tregionsRL((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytfindPointInRegionss  cCs^|d|jkoI|d|jkoI|d|jkoI|d|jk}|sZdGHn|S(NiisINFO: point not in range(R4R6R5R7(R R0tresult((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt pointInRange”sLcCs‡|j|ƒsdS|j|ƒ}dGHdGt|ƒGH|j||ƒ}|dkrƒdGHdGt|jƒGH|j||jƒ}n|S(Ns%INFO: find point in candidate regionssINFO: length of candidate sINFO: find point in all regionss&INFO: length of transportation regions(R]RYRSR'R[RI(R R0t candidatesR\((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pytfindPointTransRegion™s cCs+|j||ƒ|j|ƒ|jƒdS(N(RFRMRP(R R9R:RJ((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyt initialize¦s ( RRRRFRMRPRQRSRNR[R]R_R`(((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyR3Os       (((( tshapely.geometryRt scipy.spatialRtjsontpdbRRR3(((sH/home/zf72/Dropbox/projects/off-peak-trans/src/exploring/TransRegions.pyts   #&