Changeset 7127


Ignore:
Timestamp:
01/20/11 15:25:26 (14 years ago)
Author:
jschierm
Message:

kml: replaced mapxy by xy2ll in all kml functions.

Location:
issm/trunk/src/m/kml
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/m/kml/kml_mesh_elem.m

    r6785 r7127  
    139139    kring.coords    =zeros(size(md.elements,2)+1,3);
    140140
    141         if ~isempty(md.lat) & ~isempty(md.long),
    142                 for j=1:size(md.elements,2)
    143                         lat=md.lat(md.elements(i,j)); long=md.long(md.elements(i,j));
    144                         kring.coords(j,:)=[long lat alt];
    145                 end
    146         else
    147                 for j=1:size(md.elements,2)
    148                         [lat,long]=mapxy(md.x(md.elements(i,j)),md.y(md.elements(i,j)),'s');
    149                         kring.coords(j,:)=[long lat alt];
    150                 end
    151         end
     141    for j=1:size(md.elements,2)
     142        [lat,long]=xy2ll(md.x(md.elements(i,j)),md.y(md.elements(i,j)),-1);
     143        kring.coords(j,:)=[long lat alt];
     144    end
    152145    kring.coords(end,:)=kring.coords(1,:);
    153146

  • issm/trunk/src/m/kml/kml_part_edges.m

    r6473 r7127  
    218218%  if first edge, write out first node
    219219                    if ~elast
    220                         [lat,long]=mapxy(md.x(edgeper(j,1)),md.y(edgeper(j,1)),'s');
     220                        [lat,long]=xy2ll(md.x(edgeper(j,1)),md.y(edgeper(j,1)),-1);
    221221                        kline.coords(end+1,:)=[long lat alt];
    222222                    end
    223                     [lat,long]=mapxy(md.x(edgeper(j,2)),md.y(edgeper(j,2)),'s');
     223                    [lat,long]=xy2ll(md.x(edgeper(j,2)),md.y(edgeper(j,2)),-1);
    224224                    kline.coords(end+1,:)=[long lat alt];
    225225                    elast=elemper(j);
    … …  
    309309%  write out first node of first side for half-edge to midpoint
    310310%                                disp(['segment j=' int2str(j) ' unshared half edge from node ' int2str(elemp(ielem,nlast)) ' (node ' int2str(nlast) ') on side ' int2str(slast) ' from element ' int2str(ielem) ' written.'])
    311                                 [lat,long]=mapxy((md.x(elemp(ielem,nlast))),...
    312                                                  (md.y(elemp(ielem,nlast))),'s');
     311                                [lat,long]=xy2ll((md.x(elemp(ielem,nlast))),...
     312                                                 (md.y(elemp(ielem,nlast))),-1);
    313313                                kline.coords(end+1,:)=[long lat alt];
    314314                            end
    … …  
    316316                           
    317317%  write out midpoint of first side
    318                             [lat,long]=mapxy((md.x(elemp(ielem,slast))...
     318                            [lat,long]=xy2ll((md.x(elemp(ielem,slast))...
    319319                                             +md.x(elemp(ielem,mod(slast,3)+1)))/2.,...
    320320                                             (md.y(elemp(ielem,slast))...
    321                                              +md.y(elemp(ielem,mod(slast,3)+1)))/2.,'s');
     321                                             +md.y(elemp(ielem,mod(slast,3)+1)))/2.,-1);
    322322                            kline.coords(end+1,:)=[long lat alt];
    323323                        end
    … …  
    359359%  write out half-edge from current node to midpoint of unshared side
    360360%                            disp(['segment j=' int2str(j) ' unshared half edge from node ' int2str(elemp(ielem,nlast)) ' (node ' int2str(nlast) ') on side ' int2str(slast) ' from element ' int2str(ielem) ' written.'])
    361                             [lat,long]=mapxy((md.x(elemp(ielem,nlast))...
     361                            [lat,long]=xy2ll((md.x(elemp(ielem,nlast))...
    362362                                             +md.x(elemp(ielem,nnext)))/2.,...
    363363                                             (md.y(elemp(ielem,nlast))...
    364                                              +md.y(elemp(ielem,nnext)))/2.,'s');
     364                                             +md.y(elemp(ielem,nnext)))/2.,-1);
    365365                            kline.coords(end+1,:)=[long lat alt];
    366366                            nlast=0;
    … …  
    393393%  all different, so cut through centroid
    394394%                                disp(['element ielem=' int2str(ielem) ' centroid written.'])
    395                                 [lat,long]=mapxy(sum(md.x(elemp(ielem,:)))/3.,...
    396                                                  sum(md.y(elemp(ielem,:)))/3.,'s');
     395                                [lat,long]=xy2ll(sum(md.x(elemp(ielem,:)))/3.,...
     396                                                 sum(md.y(elemp(ielem,:)))/3.,-1);
    397397                                kline.coords(end+1,:)=[long lat alt];
    398398                            end
    … …  
    409409%  write out midpoint of opposite side
    410410%                        disp(['segment j=' int2str(j) ' internal edge from side ' int2str(slast) ' to side ' int2str(snext) ' from element ' int2str(ielem) ' written.'])
    411                         [lat,long]=mapxy((md.x(elemp(ielem,snext))...
     411                        [lat,long]=xy2ll((md.x(elemp(ielem,snext))...
    412412                                         +md.x(elemp(ielem,mod(snext,3)+1)))/2.,...
    413413                                         (md.y(elemp(ielem,snext))...
    414                                          +md.y(elemp(ielem,mod(snext,3)+1)))/2.,'s');
     414                                         +md.y(elemp(ielem,mod(snext,3)+1)))/2.,-1);
    415415                        kline.coords(end+1,:)=[long lat alt];
    416416                        elast=ielem;
    … …  
    442442                            end
    443443%                            disp(['segment j=' int2str(j) ' unshared half edge on side ' int2str(slast) ' to node ' int2str(elemp(elast,nnext)) ' (node ' int2str(nnext) ') from element ' int2str(elast) ' written.'])
    444                             [lat,long]=mapxy(md.x(elemp(elast,nnext)),...
    445                                              md.y(elemp(elast,nnext)),'s');
     444                            [lat,long]=xy2ll(md.x(elemp(elast,nnext)),...
     445                                             md.y(elemp(elast,nnext)),-1);
    446446                            kline.coords(end+1,:)=[long lat alt];
    447447                            break

  • issm/trunk/src/m/kml/kml_part_elems.m

    r6473 r7127  
    183183
    184184            for j=1:size(elem,2)
    185                 [lat,long]=mapxy(md.x(elem(i,j)),md.y(elem(i,j)),'s');
     185                [lat,long]=xy2ll(md.x(elem(i,j)),md.y(elem(i,j)),-1);
    186186                kring.coords(j,:)=[long lat alt];
    187187            end

  • issm/trunk/src/m/kml/kml_part_flagedges.m

    r6461 r7127  
    8383
    8484        for j=1:2
    85             [lat,long]=mapxy(xseg(i,j),yseg(i,j),'s');
     85            [lat,long]=xy2ll(xseg(i,j),yseg(i,j),-1);
    8686            kline.coords(j,:)=[long lat alt];
    8787        end

  • issm/trunk/src/m/kml/kml_partitions.m

    r6473 r7127  
    220220%  if first edge, write out first node
    221221                    if ~elast
    222                         [lat,long]=mapxy(md.x(edgeper(j,1)),md.y(edgeper(j,1)),'s');
     222                        [lat,long]=xy2ll(md.x(edgeper(j,1)),md.y(edgeper(j,1)),-1);
    223223                        kring.coords(end+1,:)=[long lat alt];
    224224                    end
    225                     [lat,long]=mapxy(md.x(edgeper(j,2)),md.y(edgeper(j,2)),'s');
     225                    [lat,long]=xy2ll(md.x(edgeper(j,2)),md.y(edgeper(j,2)),-1);
    226226                    kring.coords(end+1,:)=[long lat alt];
    227227                    elast=elemper(j);
    … …  
    311311%  write out first node of first side for half-edge to midpoint
    312312%                                disp(['segment j=' int2str(j) ' unshared half edge from node ' int2str(elemp(ielem,nlast)) ' (node ' int2str(nlast) ') on side ' int2str(slast) ' from element ' int2str(ielem) ' written.'])
    313                                 [lat,long]=mapxy((md.x(elemp(ielem,nlast))),...
    314                                                  (md.y(elemp(ielem,nlast))),'s');
     313                                [lat,long]=xy2ll((md.x(elemp(ielem,nlast))),...
     314                                                 (md.y(elemp(ielem,nlast))),-1);
    315315                                kring.coords(end+1,:)=[long lat alt];
    316316                            end
    … …  
    318318                           
    319319%  write out midpoint of first side
    320                             [lat,long]=mapxy((md.x(elemp(ielem,slast))...
     320                            [lat,long]=xy2ll((md.x(elemp(ielem,slast))...
    321321                                             +md.x(elemp(ielem,mod(slast,3)+1)))/2.,...
    322322                                             (md.y(elemp(ielem,slast))...
    323                                              +md.y(elemp(ielem,mod(slast,3)+1)))/2.,'s');
     323                                             +md.y(elemp(ielem,mod(slast,3)+1)))/2.,-1);
    324324                            kring.coords(end+1,:)=[long lat alt];
    325325                        end
    … …  
    361361%  write out half-edge from current node to midpoint of unshared side
    362362%                            disp(['segment j=' int2str(j) ' unshared half edge from node ' int2str(elemp(ielem,nlast)) ' (node ' int2str(nlast) ') on side ' int2str(slast) ' from element ' int2str(ielem) ' written.'])
    363                             [lat,long]=mapxy((md.x(elemp(ielem,nlast))...
     363                            [lat,long]=xy2ll((md.x(elemp(ielem,nlast))...
    364364                                             +md.x(elemp(ielem,nnext)))/2.,...
    365365                                             (md.y(elemp(ielem,nlast))...
    366                                              +md.y(elemp(ielem,nnext)))/2.,'s');
     366                                             +md.y(elemp(ielem,nnext)))/2.,-1);
    367367                            kring.coords(end+1,:)=[long lat alt];
    368368                            nlast=0;
    … …  
    395395%  all different, so cut through centroid
    396396%                                disp(['element ielem=' int2str(ielem) ' centroid written.'])
    397                                 [lat,long]=mapxy(sum(md.x(elemp(ielem,:)))/3.,...
    398                                                  sum(md.y(elemp(ielem,:)))/3.,'s');
     397                                [lat,long]=xy2ll(sum(md.x(elemp(ielem,:)))/3.,...
     398                                                 sum(md.y(elemp(ielem,:)))/3.,-1);
    399399                                kring.coords(end+1,:)=[long lat alt];
    400400                            end
    … …  
    411411%  write out midpoint of opposite side
    412412%                        disp(['segment j=' int2str(j) ' internal edge from side ' int2str(slast) ' to side ' int2str(snext) ' from element ' int2str(ielem) ' written.'])
    413                         [lat,long]=mapxy((md.x(elemp(ielem,snext))...
     413                        [lat,long]=xy2ll((md.x(elemp(ielem,snext))...
    414414                                         +md.x(elemp(ielem,mod(snext,3)+1)))/2.,...
    415415                                         (md.y(elemp(ielem,snext))...
    416                                          +md.y(elemp(ielem,mod(snext,3)+1)))/2.,'s');
     416                                         +md.y(elemp(ielem,mod(snext,3)+1)))/2.,-1);
    417417                        kring.coords(end+1,:)=[long lat alt];
    418418                        elast=ielem;
    … …  
    444444                            end
    445445%                            disp(['segment j=' int2str(j) ' unshared half edge on side ' int2str(slast) ' to node ' int2str(elemp(elast,nnext)) ' (node ' int2str(nnext) ') from element ' int2str(elast) ' written.'])
    446                             [lat,long]=mapxy(md.x(elemp(elast,nnext)),...
    447                                              md.y(elemp(elast,nnext)),'s');
     446                            [lat,long]=xy2ll(md.x(elemp(elast,nnext)),...
     447                                             md.y(elemp(elast,nnext)),-1);
    448448                            kring.coords(end+1,:)=[long lat alt];
    449449                            break

  • issm/trunk/src/m/kml/kml_unsh_edges.m

    r6461 r7127  
    8989
    9090        for j=1:2
    91             [lat,long]=mapxy(md.x(edgeuns(i,j)),md.y(edgeuns(i,j)),'s');
     91            [lat,long]=xy2ll(md.x(edgeuns(i,j)),md.y(edgeuns(i,j)),-1);
    9292            kline.coords(j,:)=[long lat alt];
    9393        end
Note: See TracChangeset for help on using the changeset viewer.