Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1624 1258 (N16W43 L=335)   391    105        0 (1575/1225-1662/1281) -33.7 W39-W45 N13-N17
      2 1340 1148 (N11W20 L=313)   349     65        1 (1318/1119-1368/1173) -45.0 W18-W22 N08-N12
      3 1428 1101 (N07W26 L=319)   295     69        3 (1397/1075-1467/1123)  26.6 W23-W28 N05-N08
      4  397 1065 (N07E42 L=250)   282    100        0 ( 363/1027- 433/1102) -33.7 E45-E39 N05-N09
      5 1844  884 (S09W62 L=355)   803    118        3 (1813/ 833-1867/ 924)  90.0 W57-W65 S12-S06
      6  468  853 (S06E37 L=255)   805    177        1 ( 381/ 804- 531/ 884) -18.4 E44-E32 S08-S03
      7 1897  824 (S14W72 L=  5)   823     84        1 (1886/ 781-1909/ 864) -68.2 W69-W75 S16-S11
      8  274  623 (S21E62 L=231)   897    138        2 ( 256/ 564- 300/ 686)  83.7 E66-E58 S25-S16
      9  543  604 (S22E35 L=257)   257    101        0 ( 522/ 560- 561/ 646)  63.4 E37-E33 S24-S18
     10  791  534 (S27E18 L=275)   275     61        1 ( 774/ 503- 813/ 565)  51.3 E19-E15 S29-S24
     11 1641  503 (S33W50 L=342)   673     97        0 (1606/ 475-1685/ 540) -22.6 W44-W56 S35-S30
  

    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