Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1364 1764 (N55W42 L=163)   639     88        2 (1325/1730-1398/1794)  14.0 W35-W51 N51-N57
      2 1327 1565 (N38W26 L=147)   338     60        1 (1295/1538-1352/1581) -14.0 W23-W28 N36-N39
      3  588 1426 (N30E32 L= 90)  1875    440        2 ( 436/1292- 738/1539) -26.6 E50-E18 N21-N37
      4 1578 1496 (N32W46 L=167)   307     52        1 (1551/1473-1610/1517)   8.1 W42-W49 N30-N33
      5  226 1288 (N20E64 L= 57)   921    174        4 ( 181/1234- 271/1345)  54.5 E69-E58 N16-N23
      6 1090 1251 (N18W05 L=126)   510    252        0 (1049/1169-1165/1314)  48.8 W01-W09 N12-N21
      7  715  820 (S08E20 L=101)  2793    808        4 ( 561/ 503- 900/1102)  53.6 E36-E07 S29-N09
      8 1502  919 (S04W31 L=152)   301     72        2 (1480/ 888-1530/ 961) -63.4 W29-W32 S05-S01
      9 1432  514 (S31W29 L=151)  3261    624        6 (1200/ 378-1624/ 701)  37.2 W12-W42 S40-S18
     10  972  671 (S18E04 L=117)   225     51        2 ( 953/ 650- 995/ 694)  33.7 E05-E02 S19-S16
     11  407  631 (S21E46 L= 75)   408     74        2 ( 384/ 596- 431/ 664)  39.8 E49-E43 S23-S18
     12 1015  484 (S31E02 L=119)  1283    454        0 ( 878/ 390-1214/ 616)  34.8 E12-W11 S38-S22
  

    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