Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1012 1726 (N36E19 L= 89)   799    167        3 ( 935/1682-1086/1761)  -8.7 E25-E12 N34-N37
      2  373 1556 (N13E58 L= 50)   952    119        4 ( 330/1502- 415/1600)  51.3 E61-E55 N09-N16
      3 1200 1510 (N25E00 L=108)   257     87        0 (1170/1472-1235/1551)  80.5 E01-W01 N22-N28
      4 1549 1515 (N35W24 L=132)   842    170        1 (1470/1465-1620/1547)   0.0 W17-W30 N32-N37
      5 1033 1267 (N07E04 L=103)   217     52        1 (1005/1250-1066/1284)  45.0 E05-E02 N05-N08
      6  185 1221 (S09E64 L= 44)   897    129        3 ( 155/1160- 228/1273)  49.4 E67-E60 S12-S05
      7  304 1144 (S13E50 L= 58)   384     99        2 ( 262/1122- 353/1172)  -7.1 E53-E45 S13-S11
      8 1808 1103 (N16W50 L=158)   714    109        0 (1763/1059-1841/1144) -18.4 W45-W54 N14-N17
      9 1410 1052 (N03W21 L=129)   304    104        0 (1398/ 998-1434/1102) -59.0 W19-W23 N00-N05
     10  868  982 (S13E08 L=100)   993    515        1 ( 665/ 904-1064/1061)  31.8 E19-W02 S20-S05
     11  285  961 (S24E50 L= 58)   414     64        1 ( 253/ 937- 319/ 989)  26.6 E52-E46 S25-S21
     12 1827  836 (N02W58 L=166)  1599    323        2 (1771/ 684-1860/ 977) -80.0 W54-W61 S07-N10
     13 1273  739 (S18W21 L=129)   232     51        1 (1254/ 715-1300/ 760)  56.3 W19-W22 S19-S16
     14 1111  693 (S24W13 L=121)   246     79        0 (1078/ 677-1156/ 715)   0.0 W10-W15 S24-S22
     15 1776  587 (S12W65 L=173)  2830    259       11 (1709/ 516-1854/ 676)  42.3 W56-W78 S17-S07
     16 1144  615 (S28W17 L=125)   212     43        1 (1118/ 596-1175/ 631)  56.3 W15-W18 S29-S26
     17 1664  402 (S25W69 L=177)  1044    102        1 (1625/ 380-1727/ 419)  14.0 W63-W79 S27-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