Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1262 1853 (N65W45 L=349)  1043    116        4 (1207/1830-1314/1878)   7.1 W32-W64 N62-N66
      2 1554 1576 (N37W49 L=353)   440     60        1 (1520/1560-1587/1591)   0.0 W44-W52 N36-N38
      3  970 1455 (N32E02 L=302)   322    110        0 ( 925/1428-1010/1472)   0.0 E05-W00 N30-N33
      4  922 1378 (N27E06 L=299)   430    122        2 ( 899/1319- 949/1435) -74.1 E07-E03 N23-N30
      5 1557 1382 (N24W40 L=345)   331     74        0 (1531/1346-1578/1418) -78.7 W38-W41 N21-N26
      6 1005 1304 (N22E00 L=304)   282    117        0 ( 952/1291-1063/1317)  -9.5 E03-W03 N20-N22
      7  235 1119 (N11E59 L=246)   448     47        0 ( 221/1090- 253/1150) -76.0 E60-E56 N08-N12
      8 1748  852 (S10W51 L=356)   643     63        1 (1725/ 815-1778/ 880) -45.0 W48-W54 S12-S08
      9  800  816 (S08E15 L=290)   781    217        2 ( 705/ 772- 861/ 859)  26.6 E21-E10 S10-S05
     10  338  651 (S19E53 L=252)   606     77        2 ( 309/ 621- 378/ 674)  20.6 E57-E48 S20-S17
     11 1559  590 (S26W38 L=342)   279     79        1 (1532/ 574-1595/ 617)  -9.5 W35-W41 S27-S24
     12 1174  153 (S67W17 L=321)   786     63        1 (1148/ 119-1190/ 184) -36.4 W10-W29 S75-S61
  

    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