Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  511 1751 (N24E65 L= 66)   774     63        2 ( 481/1729- 541/1776)  63.4 E70-E60 N21-N26
      2  827 1612 (N22E30 L=101)  1129    227        2 ( 793/1513- 871/1711) -55.3 E34-E24 N15-N28
      3 1448 1629 (N44W09 L=140)   851    118        1 (1386/1599-1512/1652)  10.3 W03-W14 N41-N46
      4 1273 1540 (N31E01 L=130)   456     70        1 (1241/1511-1312/1566)   0.0 E03-W02 N28-N32
      5 1327 1321 (N19W09 L=140)   883    348        2 (1225/1271-1467/1401)  14.0 W01-W17 N15-N23
      6 1815 1316 (N37W49 L=180)   776     80        1 (1790/1275-1839/1351) -26.6 W45-W52 N34-N38
      7 1773 1170 (N25W45 L=176)   321     79        0 (1730/1141-1806/1188)  40.6 W40-W48 N21-N27
      8 1639 1046 (N13W36 L=167)   823    225        0 (1554/ 995-1709/1136) -11.3 W27-W42 N11-N15
      9  201 1051 (S25E61 L= 70)   399     85        0 ( 181/1006- 212/1089) -51.3 E63-E58 S27-S22
     10 1775 1039 (N17W47 L=178)   437     78        0 (1739/1017-1801/1069)  38.7 W44-W49 N15-N19
     11  785  805 (S25E08 L=123)  1321    211        5 ( 699/ 747- 878/ 867)  54.0 E12-E02 S30-S19
     12 1873  828 (N10W66 L=197)   466    107        0 (1827/ 813-1908/ 862)  31.0 W59-W70 N07-N13
     13 1693  640 (S06W55 L=186)   435     70        1 (1655/ 627-1731/ 653)  33.7 W52-W58 S08-S04
     14  741  615 (S38E05 L=126)   554    132        1 ( 698/ 590- 813/ 647)  26.6 E08-W00 S39-S35
     15  567  565 (S47E18 L=113)  1452    225        6 ( 507/ 514- 644/ 638)  -8.1 E24-E10 S50-S41
     16 1487  539 (S19W45 L=176)   639    114        2 (1452/ 508-1544/ 582)   6.3 W40-W49 S20-S16
     17  847  518 (S41W06 L=137)   772    210        1 ( 749/ 480- 933/ 576)  50.7 W00-W09 S46-S34
     18 1712  460 (S14W72 L=202)  3619    196        5 (1653/ 386-1781/ 522)  70.0 W63-W83 S19-S08
     19 1581  438 (S20W59 L=190)   432     87        1 (1540/ 406-1614/ 472)  90.0 W57-W60 S23-S17
  

    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