Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1019 1862 (N48E37 L=220)   715     79        1 ( 988/1824-1054/1885)  90.0 E40-E32 N44-N51
      2  933 1789 (N38E35 L=221)   755    135        1 ( 917/1714- 958/1845) -45.0 E40-E29 N32-N43
      3 1279 1538 (N32E00 L=256)   612    112        1 (1260/1468-1294/1592) -63.4 E02-W02 N27-N35
      4 1516 1462 (N35W18 L=275)  1376    357        1 (1432/1302-1605/1592) -37.9 W09-W27 N26-N42
      5 1667 1083 (N17W37 L=294)  1253    251        6 (1625/1002-1750/1192) -21.0 W32-W45 N12-N21
      6  915  547 (S36W09 L=265)   774    257        1 ( 867/ 453- 947/ 657) -71.6 W05-W12 S43-S28
      7  731  607 (S38E06 L=251)   385     98        1 ( 698/ 575- 767/ 643)  66.8 E07-E03 S41-S34
      8 1163  334 (S41W37 L=294)   347     74        0 (1130/ 312-1197/ 360)   0.0 W33-W41 S41-S39
      9 1420  325 (S32W56 L=313)   901    100        1 (1369/ 299-1465/ 350)  11.3 W50-W61 S32-S29
  

    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