Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1306 1291 (N14W24 L=279)  1363    369        2 (1187/1220-1471/1320) -11.3 W14-W34 N09-N17
      2 1839 1033 (S17W54 L=309)  1303    188        1 (1794/ 990-1894/1097)  45.0 W49-W60 S20-S12
      3  642  995 (N15E23 L=232)  1999    450        5 ( 453/ 957- 839/1041) -21.8 E36-E11 N08-N19
      4 1714  832 (S24W39 L=294)   780    159        0 (1670/ 754-1743/ 916)  41.2 W33-W42 S28-S19
      5  327  848 (N13E50 L=205)   285     57        0 ( 312/ 823- 349/ 871)  45.0 E51-E47 N11-N14
      6 1187  801 (S10W03 L=258)   607    172        0 (1138/ 715-1251/ 861)  32.0 E01-W08 S14-S07
      7 1278  784 (S13W08 L=263)   407    149        0 (1248/ 711-1318/ 858)  36.9 W04-W12 S17-S10
      8  332  693 (N04E55 L=200)  1246    188        1 ( 293/ 625- 384/ 797) -84.3 E57-E51 N00-N10
      9 1693  649 (S36W35 L=290)   693    160        1 (1641/ 611-1729/ 671) -45.0 W30-W38 S40-S32
     10 1292  560 (S28W03 L=258)   358     57        0 (1257/ 534-1322/ 592)  21.8 W00-W05 S28-S26
     11  707  414 (S21E38 L=217)  2228    362        3 ( 611/ 263- 787/ 540)  77.7 E43-E35 S33-S10
     12  578  347 (S23E51 L=204)   709    208        0 ( 558/ 283- 602/ 405)  41.6 E55-E45 S27-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