Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1453 1733 (N52W15 L=171)   441     54        0 (1422/1716-1478/1753)  45.0 W11-W17 N48-N54
      2  551 1626 (N20E47 L=110)   349     93        2 ( 527/1589- 573/1670) -35.5 E51-E43 N17-N22
      3  820 1635 (N25E28 L=128)   247     47        0 ( 794/1619- 846/1651)  45.0 E29-E26 N23-N26
      4 1430 1457 (N28W15 L=172)   292    111        0 (1390/1414-1467/1499) -26.6 W12-W18 N26-N30
      5 1230  689 (S21W20 L=177)   298     84        0 (1203/ 649-1255/ 730)  90.0 W19-W21 S23-S18
      6  907  645 (S32W02 L=158)   542    230        0 ( 860/ 600- 993/ 729) -21.8 E02-W08 S34-S26
      7 1254  565 (S28W27 L=183)   376     90        1 (1217/ 530-1287/ 598) -18.4 W22-W29 S29-S25
      8  816  318 (S57W06 L=163)   328     30        0 ( 797/ 297- 834/ 342)  90.0 W04-W08 S59-S55
  

    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