Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  974 1700 (N41E22 L= 89)   355     47        0 ( 949/1679-1003/1716)   0.0 E24-E19 N39-N42
      2  373 1641 (N24E72 L= 39)  2848    146        6 ( 319/1584- 416/1698) -82.9 E81-E66 N19-N28
      3  966 1409 (N20E12 L= 99)   517    108        1 ( 929/1355-1000/1464) -50.2 E14-E08 N17-N23
      4  367 1280 (N01E49 L= 62)   538    143        0 ( 323/1229- 428/1325)  45.0 E51-E45 S02-N04
      5 1014 1304 (N15E06 L=105)   287     54        1 ( 987/1278-1047/1325)  33.7 E07-E04 N13-N16
      6  643 1107 (S04E25 L= 86)   652    133        1 ( 587/1073- 698/1129)  33.7 E28-E21 S06-S01
      7 1854  964 (N12W62 L=173)   609    116        0 (1842/ 905-1871/1016)  81.9 W60-W64 N08-N15
      8 1039  832 (S12W05 L=116)   933    361        2 ( 910/ 717-1162/ 929)  62.1 W00-W10 S21-S04
      9  315  891 (S23E47 L= 64)   324     77        0 ( 282/ 867- 352/ 909)  23.2 E50-E43 S25-S21
     10 1403  696 (S13W31 L=142)   464    111        1 (1371/ 642-1447/ 744)  90.0 W29-W32 S16-S09
     11 1644  605 (S13W53 L=164)   346    119        0 (1614/ 561-1664/ 666) -81.9 W51-W53 S16-S09
  

    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