Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  582 1688 (N26E50 L=247)  1462    429        6 ( 529/1570- 626/1789) -31.8 E62-E41 N18-N32
      2  807 1740 (N35E38 L=259)   592     97        1 ( 765/1708- 853/1765)   0.0 E41-E33 N33-N37
      3  263 1488 (N07E71 L=226)   626    130        0 ( 251/1451- 276/1521) -51.3 E75-E67 N04-N09
      4  519 1418 (N08E42 L=255)   556     90        2 ( 481/1387- 553/1445)  -8.1 E45-E38 N05-N09
      5  225 1192 (S11E60 L=237)  1160    131        3 ( 199/1117- 261/1271) -84.3 E62-E56 S15-S05
      6  167 1125 (S16E65 L=232)  1146    148        4 ( 137/1064- 188/1166) -71.6 E67-E61 S19-S13
      7  625 1001 (S14E23 L=274)   451    106        2 ( 594/ 941- 664/1041)  90.0 E24-E20 S17-S11
      8  598  767 (S29E21 L=276)   339     66        0 ( 571/ 733- 628/ 792)  51.3 E22-E18 S31-S26
      9 1213  668 (S19W20 L=317)   235     59        0 (1181/ 654-1256/ 681)  26.6 W18-W22 S20-S17
     10  454  629 (S41E31 L=266)   705    140        2 ( 411/ 594- 503/ 667)  33.7 E35-E26 S44-S38
     11 1551  371 (S27W58 L=355)   416    113        0 (1518/ 357-1580/ 389)   9.5 W54-W61 S27-S25
  

    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