Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1372 1625 (N34W41 L=172)   289     45        0 (1343/1608-1396/1639) -33.7 W39-W43 N32-N35
      2  711 1446 (N39E12 L=119)   604     75        1 ( 664/1414- 755/1474) -39.8 E15-E08 N36-N41
      3  460 1236 (N30E34 L= 97)  1248    252        2 ( 391/1146- 527/1334) -52.4 E39-E29 N23-N37
      4  961 1252 (N22W01 L=132)   984    266        4 ( 850/1190-1087/1321) -47.5 E04-W07 N15-N28
      5 1276 1032 (N02W15 L=146)  1242    272        2 (1181/ 900-1342/1148) -82.4 W11-W17 S05-N10
      6  485 1024 (N17E34 L= 97)   363    111        0 ( 457/ 989- 525/1072) -56.3 E35-E31 N14-N20
      7 1827 1008 (S14W55 L=185)  1087    211        2 (1792/ 934-1850/1068)  49.4 W51-W58 S17-S10
      8  250  979 (N17E55 L= 76)   427     63        1 ( 232/ 948- 269/1015) -63.4 E56-E53 N15-N19
      9 1274  706 (S17W08 L=139)   884    170        2 (1199/ 628-1311/ 761) -63.4 W03-W10 S22-S13
     10  290  520 (S12E68 L= 63)  1356    278        1 ( 273/ 468- 326/ 575)  36.0 E74-E63 S15-S07
     11 1020  361 (S35E17 L=114)   855    170        1 ( 989/ 304-1050/ 430)  41.6 E21-E12 S38-S30
     12  584  327 (S29E50 L= 81)  1031    145        2 ( 523/ 294- 624/ 356) -14.0 E54-E45 S31-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