Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  480 1417 (N37E37 L=127)   457    103        1 ( 434/1398- 529/1439)   0.0 E41-E32 N35-N38
      2  898 1370 (N29E04 L=161)   277     83        0 ( 862/1355- 943/1387)  -9.5 E06-W00 N27-N30
      3  553 1134 (N18E30 L=135)   269     79        1 ( 524/1094- 573/1165) -68.2 E31-E28 N15-N20
      4  465 1080 (N15E37 L=128)   549    185        0 ( 434/1004- 520/1155) -59.0 E38-E32 N10-N20
      5  315  904 (N05E51 L=114)   909    125        2 ( 274/ 841- 359/ 950)  33.7 E55-E46 N01-N08
      6  324  514 (S20E64 L=100)   681     58        1 ( 309/ 481- 345/ 540)  29.7 E68-E59 S22-S18
      7  528  380 (S31E51 L=114)   671    164        1 ( 504/ 318- 563/ 456)  34.5 E59-E42 S36-S25
      8 1600  397 (S48W44 L=208)  1363    141        3 (1530/ 363-1662/ 434)   0.0 W34-W52 S49-S44
      9  624  288 (S39E49 L=116)  1207    144        2 ( 569/ 257- 683/ 331)  14.0 E59-E39 S42-S36
  

    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