Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1705 1458 (N10W57 L=257)   976    222        2 (1672/1375-1750/1513)  45.0 W53-W64 N05-N12
      2 1018 1218 (N15W04 L=204)   632    134        2 ( 962/1189-1095/1263) -45.0 W01-W08 N12-N18
      3  484 1162 (N25E30 L=170)   842    198        0 ( 408/1109- 569/1203) -39.3 E35-E24 N19-N28
      4 1022  602 (S20E11 L=189)  1407    389        2 ( 930/ 466-1076/ 752)  62.1 E19-E06 S26-S09
      5 1702  513 (S48W38 L=238)  1296    205        2 (1674/ 438-1720/ 583) -60.3 W34-W44 S55-S41
      6  476  404 (S19E55 L=145)   621    121        2 ( 452/ 371- 507/ 468)  74.1 E57-E52 S21-S14
      7 1339  272 (S54E01 L=198)   759    121        1 (1290/ 250-1387/ 302) -40.6 E05-W03 S58-S50
  

    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