Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1645 1501 (N21W57 L=136)   582     57        1 (1619/1473-1676/1531)   7.1 W53-W61 N19-N22
      2 1541 1438 (N20W44 L=124)   611    163        1 (1475/1404-1612/1465) -15.3 W38-W49 N17-N21
      3 1272 1200 (N12W18 L= 98)   286     75        0 (1229/1186-1313/1216) -21.8 W15-W20 N11-N13
      4  200 1160 (N26E63 L= 17)   796     91        3 ( 168/1118- 233/1197) -68.2 E68-E58 N23-N28
      5 1824 1144 (S04W59 L=138)   654    141        1 (1771/1088-1871/1187) -45.0 W53-W63 S08-N00
      6 1386  942 (S05W20 L=100)   646    248        1 (1321/ 868-1450/1054) -71.6 W17-W23 S10-N02
      7 1721  982 (S10W45 L=125)   304     94        0 (1708/ 936-1733/1029) -90.0 W43-W46 S12-S07
      8 1775  933 (S15W50 L=129)   767     82        1 (1750/ 901-1811/ 979) -56.3 W47-W53 S17-S11
      9 1904  912 (S20W67 L=146)   853     81        0 (1886/ 880-1934/ 957)  18.4 W63-W73 S22-S17
     10 1498  740 (S20W25 L=104)   449    140        1 (1482/ 676-1518/ 801)  82.9 W23-W26 S24-S15
     11  439  604 (S10E47 L= 32)  1126    235        4 ( 379/ 538- 511/ 682) -68.2 E51-E43 S14-S04
     12  866  591 (S17E18 L= 62)   446    149        0 ( 850/ 524- 884/ 653)  71.6 E19-E15 S21-S12
     13  596  525 (S17E38 L= 42)  1021    332        3 ( 554/ 402- 623/ 644)  53.1 E45-E33 S25-S09
     14  306  575 (S11E62 L= 17)  1049    105        0 ( 282/ 529- 332/ 632)  66.8 E65-E59 S13-S06
     15 1574  574 (S33W30 L=110)   958    126        2 (1511/ 545-1616/ 617)   0.0 W24-W34 S35-S30
     16  381  521 (S15E57 L= 23)   744     88        1 ( 362/ 481- 402/ 565)  71.6 E60-E54 S17-S11
     17 1154  497 (S28E02 L= 78)   436    114        2 (1132/ 443-1181/ 551)  63.4 E04-W00 S32-S24
     18  295  502 (S16E70 L= 10)  1754     94        4 ( 264/ 462- 324/ 542) -90.0 E78-E65 S18-S12
     19  689  355 (S30E39 L= 40)   336    126        0 ( 630/ 339- 745/ 367)   0.0 E44-E34 S30-S29
  

    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