Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1646 1351 (N33W33 L=124)   566     87        2 (1621/1309-1666/1397) -51.3 W30-W34 N30-N35
      2 1821 1094 (N23W53 L=145)  2532    374        1 (1731/ 958-1914/1177)  46.6 W45-W68 N12-N30
      3  939 1101 (S02E07 L= 85)   242     60        1 ( 907/1076- 965/1128) -11.3 E09-E04 S03-S01
      4  417  762 (S35E35 L= 57)  2265    314        7 ( 332/ 663- 493/ 894)  53.6 E44-E29 S44-S25
      5 1769  638 (S04W64 L=155)   706     89        1 (1750/ 592-1787/ 672) -63.4 W59-W66 S06-S02
      6 1399  589 (S18W36 L=128)   596    184        1 (1333/ 551-1468/ 657)  68.2 W32-W38 S22-S12
      7 1191  589 (S24W23 L=115)   351     90        2 (1155/ 572-1239/ 612)   9.5 W19-W25 S25-S22
      8 1482  553 (S18W43 L=135)   450     72        0 (1444/ 534-1517/ 577)  53.1 W40-W45 S19-S15
      9 1040  512 (S33W17 L=109)   809    244        2 ( 938/ 451-1144/ 575)  63.4 W13-W21 S40-S26
     10 1650  498 (S15W61 L=152)   476     56        0 (1618/ 477-1682/ 523)  68.2 W58-W62 S17-S12
     11  700  480 (S46E06 L= 86)   264     74        0 ( 681/ 448- 717/ 516) -78.7 E06-E04 S48-S43
  

    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