Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1086 1743 (N40E21 L= 95)  2198    388        8 ( 949/1626-1186/1830) -19.8 E34-E09 N34-N44
      2  487 1602 (N14E53 L= 64)  4360    855       15 ( 282/1433- 686/1740)  41.2 E66-E40 N00-N24
      3  672 1476 (N10E34 L= 83)  1030    315        3 ( 595/1359- 725/1556) -36.5 E41-E26 N05-N15
      4  742 1467 (N12E29 L= 87)   385    123        1 ( 716/1405- 770/1516) -76.0 E30-E27 N07-N15
      5 1520 1242 (N21W23 L=139)   748    117        1 (1480/1200-1565/1302) -26.6 W19-W27 N17-N23
      6  282 1201 (S14E54 L= 63)   307     60        0 ( 263/1167- 301/1236)  90.0 E54-E52 S15-S11
      7 1611  438 (S19W60 L=176)   460     92        1 (1573/ 419-1652/ 459)  45.0 W57-W62 S21-S16
      8 1241  239 (S44W48 L=165)  1200    180        2 (1179/ 217-1322/ 277)  21.8 W42-W53 S46-S40
  

    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