function ge_ldb_poly(fileDescr,plotFile)
% GE_LDB_POLY  write a kml with ldb polygons

[pathstr, name, ext, versn] = fileparts(fileDescr);

ldbFile = dir(fileDescr);
world_coasts = [];
for n = 1 : length(ldbFile)
    [p, name, ext, versn] = fileparts(ldbFile(n).name);
    ldb = readldb([setBackslash(pathstr,'set') name '.ldb']);
    [x,y,log]=convertCoordinates(ldb.x,ldb.y,'CS1.code',28992,'CS2.code',4326);

    world_coasts{n} = ge_polys(x,y,...
                           'name',name,...
                           'polyColor','9933ffff',...
                            'altitude',150000,...
                             'extrude',0,...
                          'tessellate',false);
end
                  
kmlFileName = ['Westerschelde' '.kml'];
kmlTargetDir = [''];%..',filesep,'kml',filesep];
ge_output([kmlTargetDir,kmlFileName],world_coasts,...
                              'name',kmlFileName);