Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1004 1602 (N37W11 L= 52)   326    160        0 ( 935/1583-1077/1619) -18.4 W06-W15 N34-N38
      2  605 1564 (N43E20 L= 21)   283     49        0 ( 578/1545- 631/1579)   0.0 E22-E17 N42-N44
      3 1737  844 (S22W45 L= 86)   333    113        0 (1712/ 799-1768/ 896)  39.8 W42-W48 S24-S19
      4 1657  597 (S37W39 L= 79)   398     48        0 (1618/ 582-1685/ 615)  -8.1 W34-W41 S38-S36
      5  372  469 (S20E61 L=340)  2184    199        2 ( 317/ 388- 431/ 546) -90.0 E67-E56 S26-S14
      6  568  392 (S28E48 L=353)  1788    365        5 ( 434/ 284- 703/ 493) -63.4 E54-E43 S38-S19
      7 1224  391 (S42W00 L= 41)   549     86        1 (1192/ 351-1257/ 438) -90.0 E01-W02 S45-S38
  

    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