Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  433 1529 (N19E49 L=137)  1026    162        3 ( 392/1466- 468/1601)  90.0 E53-E45 N14-N23
      2 1283 1464 (N25W10 L=197)   484     68        2 (1250/1431-1311/1503) -78.7 W08-W12 N22-N27
      3  452 1369 (N08E42 L=144)   305     71        2 ( 414/1344- 488/1395)   0.0 E45-E38 N06-N10
      4 1344 1312 (N16W16 L=202)   286     88        0 (1303/1291-1391/1334)   0.0 W13-W18 N14-N16
      5 1884 1255 (N23W64 L=250)  1974    258        2 (1786/1162-1937/1333) -16.7 W51-W72 N17-N26
      6  658 1230 (N02E25 L=161)   300    131        0 ( 607/1200- 707/1253)  -8.1 E28-E21 N00-N03
      7 1134 1165 (N04W05 L=191)   292    124        1 (1093/1110-1168/1204) -33.7 W01-W07 N00-N05
      8  785 1042 (S08E14 L=172)   208     95        0 ( 761/ 991- 808/1074)  80.5 E15-E13 S11-S05
      9 1568  865 (S07W36 L=223)   419     84        0 (1533/ 840-1614/ 883)  31.0 W34-W39 S08-S05
     10  849  775 (S23E08 L=178)   544    153        0 ( 786/ 753- 918/ 809)   0.0 E11-E03 S24-S21
     11 1490  714 (S17W35 L=221)   483    212        0 (1395/ 680-1583/ 740)  26.6 W28-W41 S20-S14
  

    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