Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  818 1768 (N32E41 L=199)   353     89        0 ( 801/1728- 838/1807) -50.2 E43-E37 N29-N35
      2  819 1519 (N16E27 L=213)   496    115        0 ( 768/1480- 859/1553) -14.0 E30-E22 N14-N17
      3 1747 1129 (N22W43 L=282)   818    180        1 (1677/1058-1797/1182) -14.0 W35-W48 N18-N23
      4  378 1009 (S24E41 L=198)   347    189        1 ( 317/ 939- 446/1068)  48.0 E45-E36 S29-S18
      5  577  921 (S24E24 L=216)   309     99        2 ( 529/ 898- 629/ 948)  39.8 E27-E20 S26-S21
      6  761  741 (S30E07 L=232)  1171    344        4 ( 619/ 683- 918/ 793)   5.2 E18-W04 S31-S27
      7  449  439 (S62E33 L=206)   393     67        1 ( 428/ 411- 474/ 470) -31.0 E36-E29 S64-S58
      8 1435  384 (S28W52 L=291)   499    102        0 (1393/ 357-1481/ 419)  71.6 W49-W53 S30-S24
  

    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