Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1189 1600 (N39E01 L=195)   290     61        0 (1169/1581-1217/1629) -33.7 E02-W01 N37-N40
      2 1486 1397 (N31W25 L=220)   728    153        4 (1427/1340-1534/1454) -26.6 W19-W28 N27-N34
      3 1326 1061 (N07W17 L=213)   238     80        0 (1298/1029-1347/1096)  63.4 W16-W18 N04-N09
      4 1546 1044 (N10W32 L=228)   693    150        1 (1486/1016-1628/1072)  21.8 W28-W38 N07-N11
      5  525  854 (S20E29 L=167)  2144    499        4 ( 353/ 722- 707/1023)  45.0 E44-E18 S31-S06
      6 1704  656 (S10W56 L=251)   370     40        0 (1680/ 634-1725/ 677)  71.6 W53-W57 S11-S08
      7 1419  629 (S17W34 L=229)   316     24        0 (1395/ 612-1447/ 650)   0.0 W31-W35 S18-S16
      8 1781  582 (S13W70 L=266)  1043    102        0 (1738/ 562-1834/ 596)  10.6 W64-W81 S14-S11
      9 1355  542 (S24W32 L=228)   667    305        0 (1226/ 494-1480/ 592)  33.7 W24-W39 S28-S18
     10 1067  238 (S53W29 L=225)  1282    214        2 ( 982/ 216-1178/ 261)  11.3 W19-W39 S57-S48
  

    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