Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  726 1700 (N48E30 L=342)   364     81        1 ( 706/1668- 743/1729) -45.0 E33-E27 N45-N50
      2  683 1617 (N41E30 L=342)   489     81        1 ( 649/1588- 710/1649)  59.0 E32-E27 N38-N43
      3 1061 1482 (N32W01 L= 14)   534    126        0 (1004/1442-1122/1515) -35.5 E02-W05 N29-N34
      4 1334 1417 (N28W21 L= 33)   263     65        1 (1311/1384-1358/1443)  36.9 W18-W22 N26-N30
      5 1587 1322 (N22W39 L= 51)  1898    416        6 (1423/1238-1740/1414) -24.6 W27-W51 N16-N28
      6  637 1325 (N20E27 L=345)   534    165        1 ( 612/1243- 677/1388) -60.9 E29-E23 N15-N24
      7 1165 1300 (N20W08 L= 21)   944    207        2 (1113/1185-1219/1387) -72.9 W05-W11 N13-N26
      8 1326 1268 (N19W19 L= 31)   708    161        0 (1265/1192-1404/1340) -48.4 W15-W24 N14-N23
      9  992 1041 (N04E02 L= 10)  1031    312        1 ( 897/ 971-1127/1159) -33.7 E08-W06 N00-N11
     10  831  920 (S04E12 L=  1)   375    157        1 ( 775/ 904- 898/ 940)   7.1 E15-E07 S05-S02
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y