Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  453 1332 (S01E42 L= 60)  1225    320        2 ( 366/1266- 570/1413) -11.3 E52-E32 S04-N02
      2 1031 1189 (N02E03 L= 99)   569    197        1 ( 990/1097-1108/1260) -36.9 E06-W02 S01-N06
      3  743 1154 (S06E19 L= 83)   651    209        2 ( 671/1053- 782/1222) -41.6 E24-E15 S11-S03
      4  636  983 (S18E22 L= 80)   344    141        0 ( 582/ 934- 703/1012)  49.4 E25-E18 S21-S14
      5 1675  919 (N04W42 L=144)   291    113        1 (1643/ 872-1720/ 963)  66.8 W40-W44 N00-N07
      6 1943  929 (N14W70 L=172)  1153    169        2 (1922/ 903-1966/ 963) -76.0 W66-W75 N12-N16
      7 1147  875 (S13W10 L=112)   241    121        0 (1084/ 854-1200/ 900)   8.1 W06-W13 S13-S11
      8 1152  692 (S23W16 L=118)   192    113        0 (1091/ 685-1207/ 703)  29.7 W12-W19 S24-S20
      9 1364  591 (S23W33 L=135)   275    128        0 (1342/ 553-1398/ 633) -31.0 W30-W35 S24-S19
     10 1467  469 (S26W45 L=147)  2350    437        9 (1293/ 409-1653/ 570)  43.5 W35-W55 S34-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