Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  582 1718 (N23E53 L= 75)  1574    214        3 ( 535/1636- 627/1806) -19.8 E69-E44 N18-N29
      2  796 1726 (N29E38 L= 90)   828    139        0 ( 774/1649- 817/1789) -42.0 E44-E33 N24-N33
      3 1075 1603 (N29E15 L=114)  1218    394        1 ( 957/1531-1210/1697)  67.6 E22-E09 N23-N40
      4 1303 1636 (N39E02 L=127)  1277    149        3 (1246/1575-1353/1692) -12.5 E06-W03 N35-N43
      5  273 1355 (S05E60 L= 69)  2750    289        2 ( 228/1237- 327/1447) -50.7 E69-E52 S11-N00
      6  808 1169 (S03E16 L=113)  1509    431        2 ( 738/ 946- 908/1312) -41.4 E21-E04 S12-N03
      7 1628 1276 (N26W30 L=159)   754     91        2 (1579/1238-1666/1308)   8.1 W26-W33 N24-N28
      8 1329 1246 (N15W11 L=139)   254     82        0 (1310/1199-1342/1288) -59.0 W08-W12 N12-N17
      9 1179 1135 (N05W05 L=134)   451    138        1 (1159/1061-1211/1192) -69.4 W03-W07 N00-N08
     10 1237 1022 (N00W12 L=140)   529    194        3 (1204/ 934-1266/1116) -48.8 W07-W15 S04-N04
     11 1570  465 (S19W55 L=183)   949    150        2 (1496/ 445-1648/ 483)  32.5 W49-W61 S22-S15
  

    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