Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1058 1425 (N24E05 L=272)   717    163        2 ( 996/1346-1094/1489) -41.2 E10-E02 N19-N26
      2  697 1353 (N13E26 L=251)  1570    260        8 ( 600/1258- 747/1466) -85.9 E33-E21 N06-N20
      3  340 1365 (N07E54 L=224)   771    142        0 ( 324/1292- 358/1429) -69.4 E56-E50 N03-N11
      4  536 1379 (N12E38 L=239)   607    116        3 ( 499/1333- 572/1420) -51.3 E41-E35 N08-N13
      5  591 1309 (N08E32 L=245)   677    141        1 ( 537/1257- 642/1354)  45.0 E34-E28 N04-N11
      6  505  975 (S13E32 L=246)  1122    268        2 ( 459/ 849- 569/1082)  78.7 E34-E28 S20-S05
      7  657  923 (S13E21 L=257)   232     37        1 ( 633/ 907- 682/ 940)   0.0 E22-E18 S14-S12
      8  975  807 (S15W01 L=278)   558    226        0 ( 877/ 734-1060/ 861) -14.0 E05-W07 S17-S12
      9 1629  795 (S04W43 L=320)   491    148        0 (1595/ 731-1652/ 850) -50.2 W39-W45 S07-N00
     10  731  773 (S21E14 L=263)   264    111        0 ( 681/ 758- 781/ 789)  23.2 E17-E10 S22-S19
     11  803  587 (S32E06 L=271)   256     65        1 ( 763/ 567- 832/ 607)  31.0 E09-E04 S32-S29
     12 1715  497 (S20W66 L=343)   507     55        0 (1685/ 484-1748/ 513)   5.7 W61-W71 S20-S18
  

    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