Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  967 1704 (N41E23 L= 89)   381     52        1 ( 935/1681- 999/1719)   9.5 E25-E19 N39-N42
      2  369 1641 (N24E73 L= 38)  2603    126        3 ( 319/1587- 406/1696) -69.4 E82-E66 N19-N27
      3  668 1638 (N29E41 L= 71)   351    110        1 ( 630/1614- 721/1675) -11.3 E45-E35 N27-N31
      4  979 1368 (N18E10 L=102)   884    225        2 ( 925/1280-1040/1464) -48.0 E15-E04 N13-N23
      5  362 1282 (N01E50 L= 62)   579    142        0 ( 319/1225- 421/1328)  50.2 E52-E45 S02-N04
      6  642 1111 (S04E25 L= 87)   677    104        1 ( 584/1085- 694/1136)  29.7 E28-E21 S05-S01
      7 1849  956 (N11W62 L=173)   496     87        0 (1838/ 905-1862/1000) -80.5 W60-W63 N08-N14
      8 1031  834 (S12W04 L=116)   966    354        1 ( 905/ 717-1150/ 929)  59.5 E01-W09 S21-S04
      9  312  895 (S23E48 L= 64)   316     76        0 ( 278/ 873- 345/ 913)  23.2 E51-E44 S24-S21
     10 1517  796 (S05W35 L=147)   253     83        0 (1498/ 765-1529/ 828) -45.0 W33-W36 S06-S03
     11 1396  699 (S13W30 L=142)   426    117        1 (1365/ 645-1441/ 745)  81.9 W28-W32 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