Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  771 1412 (N25E12 L= 50)   855    192        1 ( 669/1364- 845/1438) -22.6 E19-E07 N21-N28
      2  228 1245 (N21E57 L=  5)  1713    276        2 ( 150/1160- 343/1332) -24.0 E72-E44 N14-N27
      3  682 1207 (N14E19 L= 43)   193     59        1 ( 655/1183- 711/1234) -33.7 E20-E17 N12-N15
      4  642 1021 (N03E23 L= 39)   238     55        0 ( 613/1001- 673/1038) -18.4 E24-E21 N01-N04
      5 1701  891 (S16W44 L=106)   957    288        3 (1628/ 770-1762/1017)  51.8 W38-W49 S22-S08
      6  941  853 (S10E07 L= 55)   300     62        2 ( 905/ 826- 972/ 876)   0.0 E08-E04 S11-S08
      7 1657  663 (S30W42 L=104)   622    122        1 (1634/ 618-1682/ 724) -74.1 W40-W45 S32-S25
      8 1744  659 (S31W52 L=114)   700     87        0 (1699/ 636-1787/ 680)   0.0 W46-W56 S32-S29
      9  361  495 (S25E59 L=  2)  1784    210        4 ( 327/ 382- 399/ 581)  30.6 E76-E52 S32-S19
  

    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