Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  756 1803 (N63W09 L= 67)   340     38        0 ( 739/1778- 773/1824)  68.2 W06-W11 N60-N65
      2  660 1752 (N62E05 L= 52)   410     34        1 ( 636/1726- 687/1772) -18.4 E07-E01 N59-N64
      3  461 1668 (N62E34 L= 23)  1654    198        5 ( 390/1611- 554/1705)   0.0 E48-E19 N58-N64
      4 1554 1274 (N05W38 L= 95)   365    167        0 (1491/1250-1599/1320) -29.7 W33-W41 N02-N07
      5 1792 1237 (S04W55 L=113)  1133    131        2 (1743/1176-1846/1271)   8.7 W49-W62 S06-S02
      6 1861 1148 (S11W60 L=117)   507     77        0 (1833/1121-1896/1184) -51.3 W57-W64 S13-S08
      7  861 1131 (N14E07 L= 51)   348    100        1 ( 806/1109- 893/1159) -11.3 E09-E04 N12-N15
      8  335 1113 (N25E43 L= 14)   564     81        1 ( 305/1074- 368/1158)  63.4 E46-E40 N23-N27
      9 1120 1073 (N05W06 L= 64)   233     26        0 (1097/1051-1139/1096) -45.0 W05-W07 N03-N06
     10  907  824 (S04E11 L= 46)   990    313        2 ( 791/ 740-1011/ 933)  16.4 E19-E02 S05-N00
     11 1695  579 (S41W37 L= 94)   607     89        0 (1660/ 546-1729/ 621)  14.0 W32-W40 S43-S38
     12  494  494 (S13E48 L=  9)   793    161        1 ( 464/ 426- 547/ 561)  78.7 E50-E45 S18-S08
     13  342  441 (S14E67 L=350)   680     71        0 ( 327/ 405- 359/ 487)  50.2 E71-E63 S16-S10
     14 1391  392 (S44W07 L= 64)   537    123        1 (1325/ 376-1441/ 414) -16.7 W01-W11 S46-S41
     15 1153  299 (S43E15 L= 43)   616    129        1 (1090/ 280-1224/ 324) -15.3 E19-E08 S45-S41
     16  928  226 (S42E35 L= 22)   499     78        0 ( 893/ 197- 972/ 252)  11.3 E40-E30 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