Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  899 1691 (N50W17 L= 97)   426    100        0 ( 874/1650- 945/1719)  12.5 W12-W21 N47-N51
      2 1660 1628 (N19W69 L=149)  1701    176        1 (1613/1545-1710/1684)  68.2 W64-W75 N13-N23
      3 1037 1428 (N29W13 L= 93)   342    132        0 ( 980/1388-1087/1460) -60.3 W10-W14 N24-N32
      4  481 1357 (N40E29 L= 52)   596    156        0 ( 454/1294- 513/1435) -68.2 E31-E26 N35-N45
      5  790 1268 (N26E07 L= 73)   534    171        1 ( 694/1245- 865/1288) -28.6 E13-E02 N23-N29
      6  413 1190 (N31E36 L= 44)   252     78        0 ( 385/1163- 429/1213)  56.3 E38-E34 N29-N32
      7 1874 1152 (S12W63 L=143)  1574    164        2 (1808/1097-1946/1188)  -6.6 W53-W79 S15-S10
      8  455 1135 (N27E34 L= 47)   490    103        3 ( 404/1101- 497/1170) -14.0 E37-E29 N24-N28
      9  545 1079 (N21E28 L= 53)   319     53        1 ( 525/1042- 564/1112)  76.0 E29-E26 N19-N23
     10  317  982 (N21E48 L= 32)   551    144        1 ( 276/ 925- 378/1023) -41.2 E51-E43 N16-N23
     11  618  941 (N12E26 L= 55)   257     92        0 ( 587/ 901- 657/ 977) -71.6 E27-E24 N08-N14
     12  496  878 (N11E36 L= 45)   672     99        4 ( 467/ 827- 527/ 929)  71.6 E37-E33 N08-N14
     13 1536  766 (S23W24 L=104)  1047    404        0 (1388/ 655-1647/ 902)  14.6 W11-W34 S25-S18
     14 1376  805 (S15W14 L= 95)   277     52        2 (1349/ 780-1400/ 825)  18.4 W12-W15 S16-S13
     15 1793  754 (S33W47 L=127)   812    102        0 (1748/ 717-1836/ 787)   5.2 W41-W52 S34-S31
     16 1223  615 (S22W00 L= 81)   394    162        0 (1171/ 562-1275/ 684)  24.0 E04-W05 S24-S19
     17  719  368 (S24E40 L= 40)  1443    344        0 ( 618/ 248- 804/ 461)  19.7 E58-E30 S31-S19
  

    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