Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1413 1358 (N19W22 L= 13)   390    110        1 (1359/1336-1476/1385)  -7.1 W17-W25 N17-N20
      2 1108 1286 (N11W03 L=354)  1087    290        1 ( 983/1187-1241/1362) -26.6 E05-W11 N06-N15
      3 1657 1290 (N18W39 L= 31)   876    171        0 (1621/1186-1682/1357) -78.7 W36-W41 N11-N21
      4 1813 1137 (N10W53 L= 44)   497    106        3 (1785/1107-1848/1174) -31.0 W50-W57 N08-N12
      5 1709 1113 (N07W44 L= 35)   265     78        0 (1674/1099-1740/1126)  11.3 W41-W46 N06-N08
      6  159 1003 (S11E64 L=288)   635    127        1 ( 137/ 957- 194/1042)  51.3 E66-E59 S13-S08
      7  925  954 (S10E05 L=346)   351    133        0 ( 859/ 931-1001/ 984)  -6.3 E09-W00 S11-S08
      8  606  730 (S26E25 L=326)  1025    352        0 ( 485/ 605- 721/ 853)  43.4 E35-E17 S34-S17
      9  420  595 (S36E42 L=309)   340     61        0 ( 392/ 580- 454/ 614)   0.0 E44-E38 S36-S34
     10 1548  511 (S30W44 L= 35)   771    218        0 (1497/ 454-1610/ 577) -24.8 W38-W51 S33-S26
  

    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