Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1228 1831 (N54W11 L=282)   507     72        0 (1208/1804-1245/1861)  66.8 W08-W12 N51-N58
      2  653 1636 (N32E33 L=238)   460     99        2 ( 620/1604- 690/1662)  -9.5 E36-E29 N29-N33
      3 1633 1536 (N34W43 L=314)   431    105        0 (1617/1483-1655/1586) -90.0 W41-W45 N30-N37
      4 1369 1247 (N11W20 L=291)   436    117        0 (1346/1200-1398/1306) -56.3 W17-W21 N08-N14
      5 1341  799 (S16W22 L=293)   236     60        0 (1316/ 772-1366/ 831)  63.4 W20-W22 S17-S13
      6 1461  704 (S20W32 L=303)   424    193        0 (1418/ 655-1510/ 756)  66.8 W29-W34 S23-S16
      7 1574  684 (S20W41 L=312)   245     63        1 (1551/ 659-1593/ 711) -45.0 W38-W41 S21-S18
      8 1319  425 (S40W30 L=301)   774     96        2 (1273/ 392-1355/ 475)  68.2 W26-W33 S42-S36
      9  991  446 (S42W03 L=274)   542     90        0 ( 937/ 431-1044/ 464)   7.1 E01-W07 S43-S40
  

    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