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

kml: replaced mapxy by xy2ll in all kml functions.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • 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
Note: See TracChangeset for help on using the changeset viewer.