Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1024 1598 (N40W19 L= 33)   277     68        1 ( 994/1569-1052/1621) -63.4 W17-W20 N38-N42
      2  470 1308 (N37E31 L=343)   260     93        0 ( 441/1272- 512/1341)  23.2 E33-E26 N35-N38
      3  706 1244 (N27E14 L=  0)   359    106        1 ( 646/1230- 766/1260) -26.6 E17-E09 N25-N29
      4 1218 1130 (N08W13 L= 27)   253     79        0 (1194/1090-1249/1159)  26.6 W11-W15 N06-N09
      5  209 1039 (N25E59 L=315)  1995    217        3 ( 125/ 969- 299/1111) -21.4 E72-E49 N20-N30
      6 1837  961 (S20W53 L= 67)   372     32        1 (1819/ 936-1852/ 986) -90.0 W50-W54 S22-S18
      7 1363  846 (S12W14 L= 28)   318    129        1 (1341/ 804-1379/ 910)  59.0 W12-W16 S14-S08
      8 1179  799 (S09W03 L= 17)   228     97        0 (1139/ 764-1221/ 825)  11.3 W00-W05 S10-S08
      9 1142  727 (S12E01 L= 13)   272     87        1 (1100/ 708-1183/ 743) -36.9 E03-W01 S13-S10
     10 1130  650 (S17E04 L= 10)   676    203        2 (1063/ 613-1203/ 710)  20.6 E08-W01 S18-S14
     11  807  622 (S10E23 L=351)   605    164        0 ( 758/ 555- 853/ 698)  33.7 E27-E18 S13-S06
     12 1676  580 (S39W34 L= 47)   361     52        0 (1656/ 556-1697/ 609)  45.0 W31-W35 S40-S37
     13  615  450 (S17E42 L=332)   483    152        1 ( 593/ 393- 641/ 527)  41.6 E45-E36 S20-S12
     14  541  476 (S13E46 L=328)   683    118        1 ( 487/ 444- 582/ 511) -45.0 E49-E43 S16-S10
     15 1468  467 (S40W14 L= 27)   292     37        1 (1443/ 446-1490/ 483) -14.0 W11-W15 S41-S37
     16  746  349 (S26E39 L=335)   508     80        1 ( 706/ 319- 784/ 378)   0.0 E42-E34 S27-S24
     17  723  249 (S33E48 L=325)  2374    368        3 ( 632/ 125- 807/ 353)  58.0 E61-E43 S47-S23
     18 1094  135 (S56E35 L=338)   928    119        1 (1037/ 118-1148/ 154) -40.6 E38-E29 S60-S51
  

    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