Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1019 1861 (N51E24 L=293)   671     79        1 ( 983/1837-1067/1886)  -5.7 E29-E19 N48-N53
      2  748 1695 (N30E34 L=283)  1396    341        2 ( 697/1551- 813/1821) -65.6 E40-E29 N19-N41
      3  964 1689 (N34E18 L=299)   362    101        0 ( 945/1640- 979/1739) -49.4 E21-E15 N30-N37
      4 1503 1511 (N33W23 L=340)   323     80        1 (1460/1496-1544/1524)  23.2 W19-W26 N31-N34
      5 1423 1387 (N23W18 L=335)   219     66        0 (1405/1354-1448/1406) -56.3 W16-W19 N20-N23
      6  615  939 (S18E24 L=293)   699    177        2 ( 560/ 867- 676/1000)  63.4 E26-E20 S23-S13
      7  425  937 (S21E38 L=279)   233     50        0 ( 396/ 916- 452/ 957)  36.9 E40-E35 S22-S19
      8  417  578 (S44E41 L=276)  1247    165        3 ( 371/ 531- 470/ 635) -50.2 E45-E35 S47-S40
      9 1763  495 (S20W73 L= 30)  1818    173        2 (1711/ 443-1819/ 547)  63.4 W65-W79 S23-S15
     10 1098  216 (S57W36 L=353)   890    134        1 (1032/ 183-1150/ 247)  -6.3 W26-W45 S59-S54
  

    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