Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  802 1607 (N47E10 L=282)   259     72        1 ( 758/1595- 832/1620)   8.1 E13-E06 N45-N47
      2 1312 1509 (N32W28 L=320)  1149    320        2 (1173/1451-1449/1552) -26.6 W18-W36 N26-N36
      3  222 1379 (N35E69 L=222)  1236     92        1 ( 177/1357- 270/1399)  -2.4 E84-E60 N33-N35
      4  748 1308 (N26E15 L=277)   463    151        2 ( 669/1286- 818/1323) -11.3 E20-E10 N24-N27
      5 1375 1228 (N13W25 L=316)   262     67        0 (1361/1192-1396/1258)  76.0 W23-W26 N11-N15
      6 1679 1156 (N04W46 L=337)   324     78        0 (1658/1132-1711/1186) -53.1 W44-W48 N02-N06
      7 1856  853 (S18W61 L=353)  4632    433        7 (1783/ 654-1925/ 996)  51.7 W55-W73 S29-S09
      8 1539  870 (S10W31 L=323)   657    171        3 (1491/ 811-1591/ 926)  81.9 W27-W35 S14-S07
      9  236  697 (S08E65 L=227)  7254    617       13 ( 163/ 527- 311/ 886)  49.0 E81-E52 S20-N03
     10 1605  747 (S19W36 L=328)   894    269        4 (1516/ 686-1698/ 811)  22.6 W28-W43 S22-S16
     11 1035  711 (S14E03 L=289)   326     94        1 ( 989/ 687-1076/ 730) -31.0 E05-W00 S15-S12
     12  978  564 (S22E09 L=283)  1441    365        4 ( 916/ 477-1096/ 699) -82.4 E12-E01 S29-S13
     13 1134  594 (S22W02 L=293)   383    120        0 (1111/ 544-1155/ 651)  66.8 W00-W03 S25-S18
     14  770  482 (S26E25 L=267)   327     77        1 ( 735/ 459- 813/ 505)  18.4 E27-E21 S27-S24
     15  762  405 (S31E28 L=264)   307     64        0 ( 733/ 380- 787/ 431)  31.0 E30-E25 S33-S29
     16 1201  345 (S43W04 L=295)   763    103        1 (1148/ 322-1256/ 380)   0.0 W00-W08 S44-S39
  

    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