Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1166 1646 (N36W30 L=103)   447     82        2 (1133/1622-1211/1667) -53.1 W26-W32 N33-N38
      2 1683 1603 (N17W65 L=137)   582    181        0 (1630/1586-1738/1625) -24.4 W60-W71 N14-N19
      3 1462 1401 (N13W37 L=109)   510    129        1 (1401/1351-1512/1431) -60.3 W33-W38 N09-N16
      4  624 1106 (N19E21 L= 51)   354     73        2 ( 582/1081- 664/1124) -36.9 E23-E18 N17-N21
      5 1903 1078 (S18W62 L=135)  1516    205        4 (1860/1036-1935/1124) -53.1 W58-W66 S21-S13
      6  186  901 (N17E61 L= 11)   523    198        1 ( 139/ 868- 235/ 935) -21.8 E66-E56 N14-N19
      7 1529  821 (S20W24 L= 97)   442    168       -1 (1499/ 789-1562/ 860)  26.6 W21-W26 S22-S17
      8 1636  643 (S35W30 L=102)  1917    286        2 (1537/ 522-1687/ 779)  42.9 W20-W35 S40-S27
      9 1096  537 (S25E09 L= 63)  1033    197        4 (1057/ 421-1138/ 643)  65.0 E12-E04 S33-S18
     10  484  450 (S15E52 L= 21)  2326    406        5 ( 406/ 335- 550/ 628)  74.1 E56-E47 S23-S02
     11  773  415 (S24E33 L= 39)  1027    255        0 ( 712/ 342- 805/ 504)  31.6 E41-E28 S27-S18
  

    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