Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  581 1699 (N59E20 L=265)  1449    194        4 ( 495/1641- 661/1754)  14.7 E32-E09 N55-N61
      2  400 1616 (N57E45 L=240)   670     78        0 ( 374/1578- 427/1646)  78.7 E49-E41 N53-N59
      3  320 1593 (N57E60 L=225)   805     73        0 ( 274/1576- 359/1609)   0.0 E75-E50 N55-N58
      4 1315 1559 (N27W33 L=318)   313     56        1 (1285/1541-1344/1579) -26.6 W30-W35 N25-N28
      5 1397 1490 (N21W35 L=320)   601    150        1 (1347/1451-1429/1530)  78.7 W32-W37 N18-N23
      6  817 1434 (N32E03 L=282)   244     56        1 ( 798/1409- 840/1466) -78.7 E04-E02 N29-N34
      7  250 1350 (N39E55 L=230)   394     77        1 ( 217/1329- 288/1371)   0.0 E60-E50 N37-N40
      8 1170 1214 (N11W13 L=297)   254     93        0 (1137/1181-1198/1247)  21.8 W09-W14 N09-N11
      9  239  759 (N04E60 L=225)   565     83        0 ( 222/ 716- 259/ 804)  80.5 E62-E58 N00-N06
     10 1636  687 (S31W34 L=318)   294     69        0 (1598/ 677-1673/ 698) -18.4 W30-W36 S32-S30
     11 1703  648 (S36W40 L=325)   326     48        1 (1685/ 627-1721/ 672)  45.0 W38-W42 S37-S34
     12  512  614 (S10E42 L=243)   341    101        0 ( 487/ 563- 539/ 657)  35.5 E45-E38 S12-S07
     13 1021  408 (S34E16 L=268)  1290    401        0 ( 871/ 325-1150/ 544)  13.0 E29-E03 S37-S27
  

    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