Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  922 1456 (N15E17 L=238)  1061    293        1 ( 877/1337- 966/1594) -49.4 E24-E12 N08-N22
      2  370 1430 (N02E53 L=201)  1636    375        2 ( 245/1362- 450/1472)   2.7 E65-E44 S01-N05
      3 1527 1314 (N23W23 L=277)   278    102        2 (1511/1265-1542/1352) -78.7 W21-W24 N19-N25
      4 1067 1242 (N06E03 L=252)   477    137        1 ( 997/1194-1121/1289) -14.0 E07-W01 N04-N07
      5  250 1104 (S19E54 L=201)  1095    265        6 ( 213/1031- 276/1220) -57.5 E58-E50 S23-S11
      6  205 1132 (S18E59 L=195)   402     73        1 ( 190/1098- 221/1159) -53.1 E60-E57 S19-S15
      7 1525 1050 (N07W28 L=283)   884    186        0 (1489/ 953-1554/1147) -45.0 W23-W32 N02-N11
      8  675  850 (S25E17 L=238)   704    259        1 ( 560/ 820- 795/ 902)  35.5 E24-E10 S29-S19
      9  277  778 (S38E48 L=206)  1498    471        1 ( 179/ 713- 406/ 813)  23.5 E59-E36 S43-S33
     10 1656  405 (S22W65 L=319)  1612    203        2 (1593/ 360-1702/ 470)  69.4 W61-W69 S25-S17
  

    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