Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1698 1464 (N27W55 L= 20)   697    204        3 (1636/1448-1758/1485) -15.3 W48-W61 N25-N28
      2 1155 1386 (N26W11 L=336)  1161    339        3 ( 967/1333-1285/1430) -21.8 E01-W19 N22-N30
      3 1425 1269 (N17W28 L=353)  1031    195        4 (1348/1212-1485/1331) -48.8 W23-W31 N13-N21
      4 1752 1266 (N14W54 L= 19)   743    161        3 (1721/1200-1786/1314) -60.3 W52-W57 N10-N17
      5  687 1191 (N16E21 L=304)   628    212        3 ( 652/1090- 727/1292) -77.0 E23-E18 N09-N22
      6 1783  888 (S09W54 L= 19)  2048    383        4 (1692/ 783-1863/ 962) -33.7 W45-W64 S16-S03
      7 1331  711 (S17W18 L=343)  1266    371        2 (1168/ 576-1458/ 868)  41.8 W07-W27 S24-S07
      8 1527  578 (S27W34 L=359)  2640    478        3 (1340/ 431-1672/ 757)  42.3 W20-W44 S37-S16
      9  333  542 (S25E56 L=269)   772    192        1 ( 280/ 511- 389/ 579) -16.7 E61-E50 S27-S22
  

    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