Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1767 1359 (N18W55 L=276)  1256    289        0 (1678/1246-1851/1470) -51.8 W49-W62 N10-N24
      2  326 1273 (N14E48 L=173)   556    141        1 ( 275/1238- 397/1290)   5.2 E52-E41 N12-N15
      3  184 1080 (N03E61 L=160)  1655    235        5 ( 142/ 983- 235/1188)  65.6 E66-E55 S02-N09
      4  122 1094 (N04E71 L=151)  1026    104        4 (  93/1047- 154/1145)  90.0 E76-E65 N01-N07
      5  508  818 (S14E34 L=187)   833    246        1 ( 487/ 702- 540/ 934)  77.9 E37-E32 S21-S07
      6 1063  863 (S13W02 L=223)   227     86        1 (1020/ 838-1088/ 893) -53.1 W00-W03 S14-S10
      7 1851  840 (S14W62 L=283)   390     57        1 (1829/ 810-1868/ 859) -36.9 W58-W64 S15-S12
      8 1034  794 (S17W00 L=221)   279    142        0 ( 974/ 777-1094/ 835) -23.2 E03-W04 S17-S14
      9 1693  552 (S32W54 L=275)   460     93        0 (1668/ 511-1725/ 591)  71.6 W52-W56 S35-S29
     10 1557  456 (S39W44 L=266)  1017    224        2 (1457/ 435-1655/ 482)   2.9 W34-W54 S40-S37
  

    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