Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1592 1380 (N11W43 L=314)   448    116        0 (1563/1316-1635/1439) -83.7 W42-W44 N06-N15
      2  372 1246 (N27E39 L=231)   838    195        1 ( 295/1204- 437/1290)  13.0 E46-E33 N24-N28
      3  409 1021 (N13E39 L=232)   294     71        0 ( 372/1001- 444/1044) -31.0 E41-E36 N10-N14
      4  826  954 (N01E13 L=257)   332     68        2 ( 796/ 920- 860/ 991) -45.0 E14-E10 N00-N03
      5  109  852 (N07E75 L=195)   978    140        0 (  95/ 812- 123/ 885)  45.0 E80-E71 N05-N09
      6  704  823 (S04E23 L=248)   341     94        0 ( 657/ 792- 742/ 864)   9.5 E25-E19 S05-S01
      7 1885  799 (S29W59 L=329)   499     97        2 (1859/ 781-1911/ 818)   9.5 W55-W62 S30-S28
      8  899  364 (S36E24 L=246)   803    167        1 ( 815/ 342- 981/ 392)  -9.5 E30-E18 S38-S33
  

    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