Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1207 1576 (N39W10 L=304)  7087   1369        9 ( 643/1398-1741/1743)  -9.8 E43-W56 N29-N47
      2  423 1341 (N17E44 L=250)   755    158        2 ( 354/1300- 458/1395)  45.0 E49-E41 N14-N20
      3 1570 1299 (N22W36 L=330)   593    189        0 (1529/1232-1590/1391)  70.0 W32-W38 N17-N28
      4  739 1192 (N10E19 L=275)  1758    423        4 ( 714/ 999- 774/1387) -82.6 E22-E16 S01-N22
      5  185 1191 (N06E66 L=229)  2206    358        2 ( 155/1037- 210/1350) -80.5 E72-E62 S03-N15
      6 1492 1300 (N22W30 L=324)   288     72        0 (1473/1261-1514/1336) -63.4 W28-W31 N19-N24
      7  287 1220 (N09E54 L=240)   933    155        3 ( 263/1161- 317/1294) -76.0 E56-E50 N05-N13
      8  340  767 (S19E47 L=247)  1718    421        2 ( 302/ 576- 366/ 959)  76.5 E52-E43 S31-S06
      9  744  645 (S23E16 L=278)   570    190        0 ( 668/ 594- 823/ 693)  36.0 E21-E10 S27-S19
     10 1576  558 (S24W44 L=338)  3688    551        8 (1338/ 449-1786/ 673)  27.3 W27-W61 S32-S15
     11 1201  414 (S37W18 L=313)   428     85        0 (1160/ 394-1247/ 433)  26.6 W15-W21 S38-S34
  

    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