Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  772 1676 (N37E33 L=155)   354     98        2 ( 740/1649- 803/1701)  14.0 E34-E29 N35-N38
      2  396 1447 (N15E52 L=136)   368     59        2 ( 380/1417- 416/1470) -71.6 E53-E50 N13-N16
      3  311 1369 (N09E57 L=131)  1448    367        7 ( 252/1249- 356/1443) -45.0 E63-E50 N02-N13
      4  644 1160 (N02E25 L=162)   269     66        1 ( 621/1127- 671/1189) -36.9 E27-E23 N00-N03
      5  259 1041 (S12E54 L=134)  1652    263        3 ( 228/ 953- 325/1150) -45.0 E58-E46 S16-S05
      6  867 1074 (S00E10 L=178)   249     98        0 ( 808/1059- 912/1091)  18.4 E13-E07 S01-N01
      7  816  986 (S06E12 L=176)   191     75        3 ( 796/ 952- 836/1024)  90.0 E12-E10 S08-S03
      8 1688  743 (S06W50 L=238)  1030    242        5 (1656/ 650-1709/ 837) -67.4 W45-W53 S12-N00
      9 1285  434 (S33W31 L=218)   397    145        0 (1220/ 419-1351/ 456)   5.7 W25-W35 S34-S31
  

    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