Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  821 1691 (N29E30 L= 16)  2044    405        3 ( 660/1576- 977/1798) -14.8 E49-E15 N24-N35
      2 1260 1753 (N46E03 L= 43)   557     80        1 (1225/1719-1293/1786)  56.3 E04-W00 N42-N49
      3  462 1668 (N22E58 L=349)  2402    244        7 ( 359/1619- 535/1736)  21.8 E69-E50 N18-N27
      4  753 1588 (N21E30 L= 16)   285    166        0 ( 689/1557- 809/1611)   7.1 E33-E25 N18-N22
      5 1532 1459 (N31W23 L= 70)  2139    490        2 (1338/1299-1663/1583) -19.7 W07-W35 N24-N36
      6 1119 1425 (N18E03 L= 43)   555    231        1 (1049/1352-1187/1509) -31.0 E08-W02 N15-N22
      7  684 1234 (S02E24 L= 22)   309     72        0 ( 653/1211- 710/1267)  71.6 E25-E22 S03-N00
      8 1808 1150 (N19W50 L= 96)  2366    288       10 (1741/1054-1900/1263) -20.2 W42-W62 N14-N24
      9 1841  965 (N10W57 L=103)  1229    176        1 (1814/ 878-1865/1044)  90.0 W53-W59 N04-N14
     10  516  927 (S22E30 L= 16)   304     83        1 ( 504/ 880- 529/ 971) -80.5 E31-E28 S25-S19
     11  163  827 (S33E68 L=338)  1055    184        3 ( 134/ 744- 213/ 900)  27.8 E79-E58 S38-S28
     12 1717  649 (S11W56 L=102)   410     33        1 (1701/ 620-1733/ 674) -56.3 W53-W58 S12-S09
     13  534  624 (S41E27 L= 19)   294     79        0 ( 495/ 604- 576/ 639)  26.6 E30-E22 S43-S39
     14 1401  519 (S27W38 L= 84)   667    163        0 (1342/ 483-1480/ 569)  -4.1 W31-W45 S28-S25
     15 1489  314 (S36W61 L=107)  1237    114        1 (1441/ 271-1539/ 348)   9.5 W53-W66 S38-S33
     16 1390  274 (S41W57 L=103)   510     52        1 (1366/ 246-1414/ 304) -11.3 W52-W62 S42-S39
  

    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