Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1721 (N54E54 L=319)   738     99        1 ( 410/1681- 485/1747)  68.2 E60-E49 N50-N56
      2 1535 1606 (N30W46 L= 60)   470    100        1 (1483/1568-1579/1637) -39.8 W42-W48 N27-N32
      3 1064 1474 (N26W08 L= 21)   291     92        1 (1022/1450-1109/1492)   9.5 W04-W10 N25-N27
      4 1772 1458 (N18W63 L= 76)   411     37        1 (1748/1445-1796/1475)   7.1 W58-W66 N17-N19
      5  783 1334 (N20E12 L=  2)   193     43        0 ( 762/1311- 807/1359) -63.4 E13-E10 N18-N22
      6 1713  643 (S31W49 L= 62)   388     47        0 (1693/ 615-1735/ 667)  56.3 W46-W50 S33-S29
      7  341  611 (S18E55 L=319)   465    101        2 ( 310/ 586- 388/ 635) -15.9 E58-E50 S19-S16
      8  393  518 (S25E54 L=319)   726    139        3 ( 352/ 466- 423/ 575)  26.6 E61-E49 S27-S21
      9 1695  411 (S49W62 L= 76)   664     71        0 (1672/ 385-1721/ 443) -14.0 W53-W73 S51-S45
     10 1635  388 (S50W53 L= 67)   512     60        1 (1600/ 366-1665/ 405) -11.3 W46-W61 S52-S48
     11 1526  349 (S52W40 L= 54)  1232    136        2 (1476/ 317-1580/ 395) -21.8 W32-W52 S55-S47
  

    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