Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1378 1594 (N38W05 L=279)   474    126        1 (1326/1553-1433/1619)   0.0 W00-W09 N36-N39
      2  406 1488 (N05E53 L=221)   409     64        2 ( 378/1456- 430/1513) -63.4 E55-E51 N02-N06
      3  474 1465 (N05E46 L=228)   352     55        1 ( 454/1427- 490/1493) -63.4 E48-E43 N02-N06
      4  929 1381 (N11E15 L=259)   902    242        1 ( 853/1294-1007/1480) -26.6 E22-E08 N07-N15
      5  316 1306 (S08E53 L=221)   691    202        0 ( 289/1219- 378/1371)  74.7 E55-E49 S13-S02
      6  194 1075 (S23E61 L=214)  1137    172        1 ( 163/ 994- 248/1143)  55.0 E64-E55 S28-S18
      7 1476  899 (S00W29 L=303)   236     71        0 (1443/ 883-1499/ 929)  45.0 W27-W30 S01-N01
      8 1028  870 (S15W04 L=279)   311     76        1 ( 990/ 844-1059/ 902)  45.0 W02-W06 S16-S12
      9 1364  727 (S13W28 L=302)   487    165        1 (1311/ 659-1395/ 797) -26.6 W22-W32 S15-S10
     10 1272  642 (S20W26 L=300)   453    149        0 (1204/ 618-1342/ 675)  11.3 W20-W30 S20-S17
     11  358  630 (S48E41 L=233)   534     83        0 ( 316/ 613- 404/ 650)  18.4 E47-E35 S49-S45
     12 1171  577 (S26W23 L=297)   382    148        0 (1102/ 545-1236/ 596)  12.5 W18-W27 S28-S25
     13 1170  495 (S31W27 L=301)   235     45        2 (1146/ 470-1192/ 521) -76.0 W25-W27 S33-S29
     14 1642  454 (S17W62 L=337)  1155    157        1 (1580/ 420-1706/ 490)  33.7 W57-W69 S21-S14
     15  875  399 (S47W13 L=287)   803    210        3 ( 792/ 360- 966/ 444)   7.1 W05-W21 S48-S44
     16 1087  310 (S45W34 L=308)   414     81        0 (1045/ 286-1125/ 332)   7.1 W29-W37 S46-S44
  

    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